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ABSTRACT 


Much  work  has  been  done  on  outliers  in  normal  populations  and 
recently  in  exponential  populations.  Here  we  extend  the  study  to  the 
examination  of  outliers  in  three  competing  life-testing  distributions. 
Assuming  the  exchangeable  model  of  random  variables,  we  examine  the 
concepts  of  outlier-prone  and  outlier-resistant  families  as  they  apply 
to  the  Gamma,  Lognormal,  and  Weibull  families  of  density  functions.  We 
examine  also  the  detection  of  outliers  for  these  distributions  and  the 
estimation  of  parameters  in  the  presence  of  one  or  more  spurious 
observations . 
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CHAPTER  I 

OUTLIERS  AND  OUTLIER-PRONENESS 
What  is  an  outlier  and  what  is  the  outlier  problem? 

To  quote  Ferguson  (1961): 

In  a  sample  of  moderate  size  taken  from  a  certain  population 
it  appears  that  one  or  two  values  are  surprisingly  far  way 
from  the  main  group.  The  experimenter  is  tempted  to  throw 
away  the  apparently  erroneous  values,  and  not  because  he  is 
certain  that  the  values  are  spurious.  ...  It  is  rather 
because  he  feels  that  other  explanations  are  more  plausible, 
and  that  the  loss  in  the  accuracy  of  the  experiment  caused  by 
throwing  away  a  couple  of  good  values  is  small  compared  to 
the  loss  caused  by  keeping  even  one  bad  value. 

Barnett  and  Lewis  (1978)  indicate  that  in  the  light  of  developments  in 
outlier  methodology  over  the  last  15  years,  Ferguson's  formulation  may 
be  unduly  restrictive.  Outlying  values  need  not  be  "bad"  or 
"erroneous";  in  fact  they  may  be  welcomed  as  indicating  a  useful 
treatment  or  a  strain  that  was  unexpectedly  good. 
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1.1  What  is  an  outlier? 

What  do  we  mean  by  an  "outlier"?  We  shall  use  Barnett  and  Lewis' 
definition  (1978):  "an  outlier  in  a  set  of  data  (is)  ...  an 
observation  (or  subset  of  observations)  which  appears  to  be 
inconsistent  with  the  remainder  of  that  set  of  data."  We  use  the  term 
"outlier"  to  characterize  an  observation  that  stands  out  in  contrast  to 
the  rest  of  the  data  and  is  not  consistent  with  what  our  mind  feels 
constitutes  a  reasonable  data  set  nor  with  our  initial  view  of  an 
appropriate  probability  model  to  describe  the  generation  of  our  data. 

An  observation  is  "spurious"  if  it  is  statistically  unreasonable  on  the 
basis  of  some  prescribed  probability  model.  This  would  include  an 
observation  known  to  be  generated  by  a  different  probability  model; 
however  a  spurious  observation  will  not  necessarily  show  up  as  an 


outlier 
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1.2  The  outlier  problem 

Experimental  scientists  and  others  who  deal  with  data  are  forced 
to  make  decisions  about  outliers  -  whether  or  not  to  include  them  in 
analysis,  whether  to  make  allowances  for  them  on  some  compromise  basis, 
etc.  What  worries  an  experimenter  is  whether  or  not  some  observations 
are  genuine  members  of  the  main  population.  If  these  observations 
appear  in  the  midst  of  the  data  they  may  not  be  conspicuous  and  are 
unlikely  to  distort  inferences  seriously.  However,  should  the  outlier 
appear  extreme,  it  could  create  problems  in  attempts  to  represent  the 
population  and  it  could  contaminate  estimates  and  tests  of  population 
parameters . 

Once  we  decide  that  outliers  exist  in  a  data  set,  we  must  decide 

on  how  to  react  to  them.  Methods  to  support  outright  rejection  or  to 

adjust  values,  prior  to  processing  the  principal  mass  of  data  will  be 
related  to  any  postulated  model  for  that  population.  We  are  concerned 
with  whether  the  extreme  values  are  so  extreme  as  to  be  unreasonable 
under  our  original  model.  This  may  indicate  the  data  contains  a 
"mistake"  (that  should  be  rejected  or  corrected)  or  it  may  indicate  an 
alternative  model  for  which  the  complete  data  set  appears  as  a 
homogeneous  sample.  The  outlier  may  be  a  'foreign*  random  influence  in 
an  otherwise  homogeneous  data  set  -  interesting  to  study  in  itself  or 

only  serving  to  obstruct  analysis  of  the  main  data  mass. 

An  assessment  of  the  discordancy  of  some  outliers  is  just  the 
first  stage  of  the  study  of  outliers.  Once  an  outlier  is  judged 
discordant  we  may: 

1)  decide  to  reject  (or  correct)  it  and  analyse  the  remaining 
(modified)  data  on  the  original  model, 


. 
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2)  decide  to  modify  the  model  to  incorporate  the  outliers  in  a 
non-discordant  manner, 

3)  refine  the  analysis  of  the  entire  data  set  to  accommodate 
outliers  (i.e.  render  the  analysis  relatively  impervious  to 
their  presence), 

4)  focus  attention  on  the  discordant  outliers  because  they 
identify  factors  of  unexpected  value. 

As  examples  of  each  of  these  situations,  consider  the  following: 

Fifteen  residuals  (about  a  simple  model)  of  observations  of 
the  vertical  semi-diameter  of  Venus,  in  seconds,  made  by  Lt.  Herndon  in 
1846  were 


-0.30 

+0.06 

in 

o 

• 

o 

l 

<r 

• 

o 

1 

+0.63 

+0.20 

+1.01 

CO 

1 — 1 

. 

o 

1 

+0.18 

+0.48 

-1.40 

+0.39 

-0.24 

-0.22 

+0.10 

Chauvenet  (1863)  declared  -1.40  and  +1.01  outliers.  We  may  choose  to 
reject  them  as  "gross  errors";  incorporate  them  by  changing  our  model 
to  a  non-normal  one  or;  accommodate  them  by  using  a  robust  estimation 
technique  that  employs  Winsorization  or  an  a-trimmed  sample. 

Consider  also  the  data  described  by  Pearson  and  Pearson  (1931)  on 


the  capacity  (in  c.c.)  of  a  sample  of  17  male  Moriori  skulls. 
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1230 

1318 

1380 

1420 

1630 

1378 

1348 

1380 

1470 

1445 

1360 

1410 

1540 

1260 

1364 

1410 

1545 

observation 

1630  is 

suspicious 

and  as 

such  may 

be  rejected  as  a 

measurement /recording  error;  or  we  may  Incorporate  it  in  a  non- 
discordant  way  by  assuming  a  non-normal  model  (since  biological  data 
often  require  skew  distributions);  or  identification  of  the  outlier  may 
reflect  the  presence  of  a  small  number  of  another  species  in  the 
population  being  studied. 
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1.3  Models  for  discordancy 

There  are  several  possible  models  for  discordancy. 

i)  Deterministic:  this  covers  the  cases  of  outliers  caused  by  obvious 

,  n 


identifiable 


gross  measurement  errors,  etc.  If  the  data  set  {x^}_^_^ 


contains  one  observation,  ,  clearly  resulting  from 

measurement /recording  error,  we  reject  H:  all  x^  e  F  in  favor  of 
H^:  all  Xj  e  F  (j^i)  and  x^  is  different  (requiring  rejection  or 

correction) . 

ii)  Inherent:  here  we  reject  H:  all  observations  are  from  F  in 
favor  of  H^:  all  observations  are  from  G,  where  G  has  more  or 
different  inherent  variability  than  F. 

iii)  Mixture:  this  model  allows  contamination  of  the  sample  by  a  few 
members  of  a  population  other  than  F.  We  reject  H:  all  x^  e  F  in 
favor  of  H^:  all  x^  £  pF+(l-p)G.  Thus  outliers  reflect  probability 
1-p  that  observations  arise  from  G. 

iv)  Slippage:  here  states  all  observations  apart  from  some  small 

number  m  arise  independently  from  the  initial  model  F  indexed  by 
location  and  scale  parameters  x  and  9,  while  the  remaining  m  are 
independent  observations  from  a  modified  version  of  F  in  which  x  or 
9  is  shifted.  Which  observation(s)  come(s)  from  the  shifted 
distribution  is  not  specified.  In  most  published  work,  F  represents 
the  normal  distribution. 


v)  Exchangeable:  this  is  an  extension  of  the  slippage  alternative,  i 
set  of  n  observations  =  {x^,...,x  }  ideally  comes  from  a  target 

population  P  with  probability  density  function  (p.d.f.)  f(x;9). 

However  the  suspicion  is  that  one  of  the  observations,  say  x 


,  is  not 
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from  the  target  population  P  but  from  a  different  population  Q  with 
p.d.f.  f(x;£).  Prior  to  the  experiment,  there  is  no  way  to  identify 
the  possible  (at  most  one  spurious  observation).  Hence,  a  priori,  each 
observation  is  equally  likely  to  be  the  discordant  one.  Thus  the 
random  variables  X.  ,...,X  are  not  independent  but  are  exchangeable. 
For  one  outlier,  the  likelihood  of  the  sample  is  given  by 

1  n 

L(x|9_,0  =  ±  l  f(x  n  f(x  ;_9)  . 

i=l  j*i  J 

Definition  1.3.1:  Let  f  be  a  family  of  absolutely  continuous 

univariate  distribution  functions  (F(x;_9)}  indexed  by  a  parameter 

9  e  0.  Then  the  random  variables  X, , . . . ,X  are  said  to  be 
—  I  n 

exchangeable  random  variables  based  on  f  if  their  joint  p.d.f.  is  of 

the  form 


l 

n 


n  F’(xi;91)F’(xr;_92), 

i  Ar 


e  0 


where  F’(x;_9)  =  ^  F(x;_9). 

2 

Anscombe  (1960)  used  this  model  for  f(x;9)  as  N(t,9  )  and  f(x;j|) 

2 

as  N(t+a9,9  ).  Guttman  and  Smith  (1969,1971)  and  Guttman  (1973a)  used 

2  2 

it  for  f(x;_9)  as  N(t,9  )  and  f(x;£)  as  N(r+a,9  ).  Kale  and 

Sinha  (1971),  Joshi  (1972b),  Kale  (1974),  Chikkagoudar  and  Kunchur 

(1980)  and  Rauhut  (1982)  have  used  it  for  f(x;_9)  as  EXP(9)  and 

f(x;p  as  EXP(  9/k*) ,  0  <  k*  <_  1. 

For  m  outliers,  we  assume  x.  ,x  ,...,x,  come  from  a  target 

Xl *  1 22  1n-m 
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population  P  with  p.d.f.  f(x;_9),  while  x.  , ...,x,  come  from 

~  1n-m+l  1n 

populations  ^i»^2,***,<^m  w^t:^  P*d.f.  f  (x;_|^ ),...,  f  (x;J^)  .  Some  or 

all  of  the  ,  i  =  l,...,m  may  be  identical  (i.e.  =  _|) .  The 

association  of  the  different  observations  with  the  different 
distributions  is  assumed  to  occur  at  random.  For  the  case  of  all 
i  =  l,...,m,  the  likelihood  of  the  sample  is 

L(x|e,c)=-i-  I  n  f(x.;o  n  f(x.;e) 

-  (")  leT  lei  1  _  3  ~ 

m 

where  I  =  (i^ , ^ • ,i  )  is  a  selection  of  m  integers  from 
{l,2,...,n}  and  T  is  the  set  of  all  fnl  possible  such  choices. 
Certain  relationships  must  exist  between  f(x;_9)  and  f(x;£)  for  it  to 
be  reasonable  that  the  suitability  of  the  model  will  be  reflected  in 
outliers.  For  the  case  of  one  outlier,  we  need  the  discordant 
observation  from  Q  to  show  up  at  one  of  the  sample  extremes.  The 
exchangeable  model  also  assumes  that  the  maximum  number  of  possible 


outliers  is  known. 
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1.4  The  concept  of  outlier-proneness 

Related  to  model  development,  though  not  actually  a  model  to 
describe  the  occurrence  of  outliers,  is  a  concept  introduced  by  Neyman 
and  Scott  (1971)  and  furthered  by  Green  (1974,1976)  and  Kale  (1975b, 
1975c,  1976).  These  papers  considered  a  method  of  distinguishing 
between  families  of  distributions  by  examining  the  extent  to  which  they 
are  liable  to  exhibit  outliers. 

Let  S  be  a  sample  of  n  >  3  observations  and  let  x/1N,...,x/  N 
n  —  (1)’  (n) 

be  the  order  statistics  for  this  sample. 


Definition  1.4.1:  For  a  positive  number  k,  x.  N  e  S  is  a  k-outlier 

(n)  n 

on  the  right  if  x,  N  -  x,  , N  >  k{x,  1N-x/1N}  .  The  definition  of  a 

(n)  (n-1 )  (n-1)  (1) 

k-outlier  on  the  left  is  analogous  (i.e  x,,  N  e  S  is  a  k-outlier  on 

(1)  n 

the  left  if  x/ON  -  x,.  .  >  k{x,  .-x/0.}. 

(2)  (1)  (n)  (2) 

Let  P(k,n|F)  denote  the  probability  that  a  sample  of  n  observations 
from  a  distribution  F  will  contain  a  krmtlier  on  the  right. 

Then,  for  jointly  distributed  random  variables 


x. 


°°  (n-1) 

P(k,n|F)  =  |  )  j 


—  oo  —  oo 


(k+1)x(n-l)'kX(l) 


g(x(l)’X(n-l)’X(n)) 


dx  .  .dx,  ,dx,  .  x 
(n)  (1)  (n-1) 


where  g(x/1N,x,  N)  is  the  marginal  joint  p.d.f.  of 

(1)’  (n-1)’  (n)'  &  j  ^ 


X(l)>  X(n-1)>  X(n)  glve"  by 


10 


X(n-1)  X(3) 

«(x(l)’x(n-l)’x(n))  =  I  •••  I  f(yl’”-’5'n)dy2,"dyn-2 

x(l)  x(l) 

where  f  is  the  joint  p.d.f.  of  X^,...,X  . 

For  the  special  case  of  i.i.d.  random  variables  we  would  have 

oo  ao  n-  2 

P(k,n|F)  =  n(n-l)  J  dF(x)j  {f  )  “  F(x)}  dF(y). 

-CO  X 

Let  J"  be  the  family  of  distributions  and  let  IL^(k,n|J”)  be  the 
least  upper  bound  of  probabilities  P(k,n|F)  for  F  e f . 

Definition  1.4.2:  A  family/  of  distributions  is  (k, n)-outlier- 

prone  on  the  right  if  II^(k,n|J*)  =  1. 

Definition  1.4.3:  A  family  f  of  distributions  is  (k,n)-outlier- 
resistant  on  the  right  if  11^  (kjnljf )  <  1. 

Definition  1.4.4:  If  a  family  J"  of  distributions  is  (k,n)-outlier- 
prone  on  the  right  Vk  >  0,  Vn  >  2  it  is  outlier-prone  completely 
on  the  right  (o.p.c.r.). 

Definition  1.4.5:  If  a  family  /  of  distributions  is  (k,n)-outlier- 
resistant  on  the  right  Vk  >  0,  Vn  >2  it  is  outlier-resistant 
completely  on  the  right.  (o.r.c.r.). 
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Theorem  1.4.6:  (Green  (1974)) 

If  J  is  a  family  of  distributions  and  is  a  random  sample  of 

n  i.i.d.  observations  from  F  e  f ,  then  f  is  outlier-prone 
completely  on  the  right  (o.p.c.r.)  iff  f  is  (k,n)-outlier-prone  on 
the  right  for  some  k  >  0,  n  >  2. 

For  i.i.d.  observations,  Theorem  1.4.6  shows  that  it  is  not 
necessary  to  distinguish  between  the  concepts  of  "(k,n)-outlier-prone 
on  the  right"  and  "outlier-prone  completely  on  the  right".  This  leads 
to  the  following  definitions  and  theorem  for  cases  in  which  the 
observations  are  i.i.d. 

Definition  1.4.7:  With  respect  to  a  random  sample  of  n  i.i.d. 
observations  a  family  f  is  outlier-prone  on  the  right  (o.p.r.)  iff  it 
is  (k,n)-outlier-prone  on  the  right  for  some  k  >  0  and  n  >  2,  or, 
equivalently,  iff  it  is  outlier-prone  completely  on  the  right. 

Definition  1.4.8:  With  respect  to  a  random  sample  of  n  i.i.d. 
observations  a  family  f  is  outlier-resistant  on  the  right  (o.r.r.) 
iff  it  is  not  o.p.r. 

Theorem  1.4.9:  With  respect  to  a  random  sample  of  n  i.i.d. 
observations  a  family  f  is  o.r.r.  iff  n^(k,n|jO  <  1  for  some 
k  >  0,  n  >  2,  F  e  f . 
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Proof:  i)  Assume  that  family  T  is  o.r.r.  .  Then  J  is  not  o.p.r.  . 
Therefore  n^(k,n|J)  +  1  Vk  >  0,  yn  >  2. 

i.e.  There  exists  at  least  one  k*  >  0,  n*  >  2,  F  e  f  ^  . 

sup  P(k* ,n’ I F)  <  1  . 

F# 

ii)  Assume  11^ (k , n |^")  <  1  for  some  k  >  0,  n>2,  F  e  J" .  Then 
n  ( k,n|J)  *  1  Vk  >  0,  Vn  >  2  . 

Therefore  f  is  not  o.p.r. 

From  Definition  1.4.8  it  follows  that  f  is  o.r.r.  . 

We  will  use  the  same  definition  of  a  (k,n)-outlier  on  the  right 
when  represents  n  observations  from  the  exchangeable  model. 

The  implication  is  that  for  an  outlier  resistant  family  we  are 
justified  in  seeking  out  and  eliminating  outliers.  On  the  other  hand, 
outlier-prone  families  should  be  used  to  model  cases  where  outliers  are 
common  and  in  these  cases  we  should  seek  ways  of  accommodating 
outliers . 

Neyman  and  Scott  (1971)  showed  that  in  general  families  differing 

only  in  location  or  scale  parameters  are  outlier-resistant.  Thus  we 

may  limit  studies  to  subfamilies  with  standard  values  for  location  and 

2 

scale.  Both  the  family  N(t,9  )  or  normal  distributions  (with  mean  % 
2 

and  variance  9  )  and  the  family  C(£,9)  of  Cauchy  distributions 
(centered  at  £  and  having  scale  parameter  9)  are  outlier- 


resistant  . 


■ 


CHAPTER  II 

SURVEY  OF  THE  LITERATURE 

Most  studies  of  the  outlier-problem  assume  an  initial  distribution 
that  is  either  exponential  or  normal. 


-13- 
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2.1  Detection  of  outliers  in  exponential  and  normal  families  using  the 
exchangeable  model. 

A  semi-Bayesian  approach  was  used  by  Kale  (1969),  Kale  and  Sinha 

(1971),  Chikkagoudar  and  Kunchur  (1980),  and  Rauhut  (1982)  with  respect 

to  the  one-parameter  exponential  family.  It  was  assumed  that 

X, ,...,X  were  such  that  n-1  observations  were  from 
1  n 

1  r  — X  l 

f(x;9.)  =  -s—  exp  {  f,  x  >  0,  9.  >0  while  one  of  the  X's  was 
1  —  l  i 

distributed  as  f(x;92),  9^  >_  9^  (i.e.  9^  =  9^/k*  ,  0<k*<l).  A 

priori,  each  X^  was  equally  likely  to  be  distributed  as  f(x;92). 

Kale  and  Sinha  (1971)  calculated  u(r;n,k*),  the  probability  that  ^(r)> 
til 

the  r  order  statistic,  corresponds  to  the  spurious  observation 
distributed  as  f(x;92).  It  was  shown  that 


u(r;n,k*) 


k*r(n)T(n-r+k*) 
F(n+k*) T(n-r+l ) 


was  monotone  increasing  in  r  and  hence  X,  N  had  maximum  posterior 

(n) 

probability  of  being  an  outlier.  Mount  and  Kale  (1973)  generalized 

this  result  for  the  case  of  n-1  observations  with  distribution 

function  F  and  one  with  distribution  function  G  where  F  and  G 

are  stochastically  ordered  (G<F),  where  a  priori  each  observation  is 

dG 

equally  likely  to  be  the  spurious  one,  and  where  ¥(x)  =  —  is 

dF 

monotone  increasing  in  r.  Kale  (1974a)  generalized  these  results  for 
the  case  of  m  (>1)  possible  outliers,  where  a  priori  each  group  of  m 


observations  (X.  ,X.  ,...,X.  )  is  equally  likely  to  come  from 

h  12  m 

distribution  function  G.  Then  (x,  11N,...,X/  N)  has  maximum 

(n-m+1)  (n) 


. 
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posterior  probability  of  corresponding  to  the  set  of  spurious 

observations,  provided  Y  is  monotone  increasing.  Kale  restricted  the 

distributions  to  the  single-parameter  exponential  family.  In  another 

paper,  Kale  (1974b)  gave  a  completely  Bayesian  approach  where  n-m 

observations  were  distributed  as  f(x;9),  m^  observations  were  from 

f(x;9_.),  j  =  l,2,...,m^,  9^  <_  9,  and  m^  observations  were  from 

f  (x;  9^) ,  1=1,2,...  9^  >_  9,  m  =  m^-t-m^  and  f  belonged  to  the 

single-parameter  exponential  family.  Under  the  exchangeable  model,  it 

was  shown  that  { (x,,  N  , . . .  ,x,  N),(x/  .)}  had  maximum 

(1)  (m)  (n-m+1)  (n) 

posterior  probability  of  being  the  set  of  spurious  observations. 

For  the  case  of  the  normal  family  of  distributions,  a  completely 

Bayesian  approach  has  been  used  by  Box  and  Tiao  (1968)  where  f(x;9)  is 

N(x,9^)  and  f(x;_9^)  is  N(r,k*9^),  k*  >_  1  (Model  A)  and  by 

2 

Dempster  and  Rosner  (1971)  where  f(x;9)  is  N(t,9  )  and  f(x;_9^)  is 
2 

N(x  ,9  ),  j  =  1,2, ...,m,  _>  x  (Model  B).  For  m  =  1,  Model  B  has 

been  handled  as  a  slippage  test  by  Ferguson  (1961),  Kudo  (1956)  and 


Paulson  (1952). 
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2.2  Estimation  in  the  presence  of  outliers. 

Anscombe  (1960)  suggested  a  premium-protection  approach  (see 
Appendix  1)  to  estimation  that  has  subsequently  been  used  by  Kale  and 
Sinha  (1971),  Joshi  (1972b),  Chikkagoudar  and  Kunchur  (1980),  and 
Rauhut  (1982)  for  estimation  of  the  mean  in  the  single-parameter 
exponential  distribution  and  by  Guttman  and  Smith  (1969,  1971,  1973a) 
and  Veale  and  Huntsberger  (1969)  for  estimation  in  the  normal 
distribution. 


2.2.1  Estimation  for  the  single-parameter  exponential  distribution. 

n 


Under  homogeneity,  the  optimal  estimator  of  9  is  T 


o,n 


y  x* 

i-i  1 
n+l  * 


MSE (T  Ik*)  =  02{4-+  -  -  as  k*  -  0. 


However,  under  the  exchangeable  model  with  one  outlier, 

r  l  .  2  fi-k* 

Ln+1  n+l  ^  k* 

Among  restricted  L-type  estimators  T(_£)  =  £  i.x...  which  ignore 

j=l  J 

the  largest  m  observations,  Kale  and  Sinha  (1971)  and  Veale  and  Kale 


n-m 


(1972)  advocated  the  one-sided  Winsorized  mean 


n-m 

T  X/j.+mx,  v 
(i)  (n-m) 

T  =  - 

m ,  n  n-m+1 

for  which  MSE (T( 1) I k*=l )  is  minimum  and  MSE(T  |k*<l)  shows  a  gain 

—  '  m, n 1 

in  efficiency  relative  to  T  for  k*  sufficiently  small.  Joshi 

o,n 

(1972(b))  considered  choice  of  m  and  showed  for  k*  small, 
substantial  gains  in  relative  efficiency  are  possible.  Samples  up  to 
size  n  =  20  were  considered. 
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Kale  (1974a)  considered  maximum  likelihood  estimation  (M.L.E.)  for 


the  exchangeable  model  with  m  possible  upper  outliers  and  obtained 
n-m 

l  xd) 

-  i=l 

9  =  -  ,  the  trimmed  mean,  as  MLE  of  9.  In  comparing  the 

m,n  n-m  r  ° 

n-1 

l  x(i> 

-  i=l 

trimmed  mean  9.  =  - = -  with  the  Winsorized  estimator 

l,n  n-1 

n-1 

X(i)+  X(n-1 ) 

T-  =  -  for  the  case  of  a  single  upper  outlier,  it  was 

l,n  n  o  rr  > 

shown  that 


MSE( 9,  |k*=l)  >  MSE(T.  |k*=l) 

1 ,  n  '  1 ,  n  1 


but  lim  MSE(9,  Ik*)  =  — t-  <  lim  MSE(T,  Ik*).  Thus  9,  provides 
k*+0  1*n  n_1  k**0  1>n  X>n 


more 


protection  than  T  but  for  a  higher  premium  (See  Appendix  I) 

l ,  n 

Rauhut  (1982)  suggested  two  ' testimators " : 


Dk*(x)  = 


o,n 


>  lf  x(n)~x(n-l)  <  k(x(n-l)'x(l)) 


n-1 

I  xd) 

i=l  K  ' 


n 


=  T  ,  otherwise 


V©  - 


T  ,  if  x.  .  <  Cx,  C  >  0 
o,n  (n) 

T  ,  otherwise 


and  showed  T  .(x)  is  preferred  over  T  ,  T  and  D  since  the 
tc  o  y  n  u  tc 

premium  is  small  compared  to  the  protection  it  provides. 
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Chikkagoudar  and  Kunchur  (1980)  proposed  using 
n 


where 


j-l 

0.  = 


x  . 
J 


2r. 

i 


*  Wi  n  ( n+1 ) 


r.  =  rank  of  x. 
l  i 


i  n 

sample.  This  estimator  is  more  efficient  than: 


n 


I 

i=l 


w .  9 
l  i 


in  the  complete 


TABLE  2.2.1 

Alternative  Estimator 

T 

T 

A 

9 

o  ,n 

m,n 

1  ,n 

When 

n>4,k*<.75 

n>6,  .  40<k*<_.  75 

all  n,  k*>.45 

to 

n=3 ,k*<. 70 

n=4 , 5 ,  . 35<k*<. 75 

use 

““ 

-  - 

qck 

n=2  ,k*<_.  65 

n=3 ,  .  30<k*<.  70 

n=2 ,  . 35<k*<. 65 

and  it  is  independent  of  k*. 

2.2.2  Estimation  for  the  normal  distribution 

Kale  (1974a)  has  shown  that  the  method  of  maximum  likelihood 

applied  to  the  exchangeable  model  involving  normal  distributions  with 

change  in  location  leads  to  estimators  that  are  trimmed  means. 

2 

For  the  case  of  n-1  observations  from  N(x,9  )  and  one 

2 

observation  from  N(x+k*9,9  )  ,  k*  >  0  ,  x  is  UMVUE,  MLE  for  x  if 

2  2 

-  -  9  k* 

k*  =  0  but  x  is  biased  and  MSE(x)  =  -  if  k*  >  0.  Thus  an 

n 

alternative  estimator  might  be  T(x)  where  T(x)  is  one  of  the 


following: 
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T(x)  = 


n-1 


i=l 


x 


(i) 


n 


n-1 


=  T 


—  {  y  x.  N  +  x.  i n  } =  T- 
"  1  L  (i)  (n-l)J  l,n 


n 


i=l 


A(x)  = 


x  ,  if  x>  .  -  x  <  C  S 
(n)  -  a 


n-1 

I  X 

i=l 


(i) 


n-1 


,  otherwise 


The  premium  paid  for  using  T  instead  of  x  is  given  by 


^y  MSE(T|k*=0)  -  -y  MSE(x|k*=0). 


The  protection  obtained  by  using  T  instead  of  x  is  given  by 


y-  MSE(x|k*>0)  -  \r  MSE(T|k*>0). 


A  detailed  study  of  the  accommodation  of  outliers  in  slippage  models 
appears  in  Guttman  and  Smith  (1969,  1971)  and  Guttman  (1973a),  using 
premium-protection.  They  consider  three  methods  for  estimating  the 


mean  p 

=  t: 

i) 

modified 

trimming 

A 

^A 

(A-rule) 

ID 

modified 

winsorization 

My  (W-rule) 

ill) 

semi-winsorization 

A 

(S-rule ) 

(see  Appendix  I). 
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Table  2.2.2.  Comparison  of  p^,  p^  and  p^ 


Form  of  (one  or  two  outliers) 


Location-slippage 

Scale-slippage 

N ( x+k* , 92 ) , k*>0 

N(x,92k*),k*>l 

Best  for  large  k* 

Not  a  contender 

A 

Best  for  intermediate  k* 

Best  for  large  k* 

A 

Best  for  small  k* 

Best  for  small  k* 

2  2 

The  results  are  based  on  sample  sizes  up  to  n  =  20,  9  -a  known, 


and  only  for  n  =  3,  when  9  is  unknown. 

-2  -2 

Guttman  and  Smith  (1971)  defined  dispersion  estimators  cr^,  and 

*22  *2. 

Og  of  a  analogously  to  p^,  p^  and  jig  (see  Appendix  I);  is 

-v  2  *2 

not  worth  considering  when  p  =  x  is  known,  cr^  and  cr^  are  not 


worth  considering  when  x  is  unknown. 


CHAPTER  III 

The  Gamma  Distribution 

One  of  the  most  common  life-testing  distributions  is  the  gamma 
distribution.  We  shall  examine  the  exchangeable  model  based  on  the 
gamma  distribution  and  show  that,  in  the  case  of  a  shape  change,  it  is 
outlier-prone  completely.  For  changes  in  shape  or  scale  parameter,  we 
shall  determine  which  observation  is  most  likely  to  be  the  spurious 
one.  We  shall  also  consider  estimation  in  the  presence  of  an  outlier. 
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3.1  Characteristics  of  the  Gamma  Distribution 

Consider  the  three-parameter  gamma  distribution  given  by 

x-t 

(  _  9 

f(x;9,p,-c)  =  -  -X  - - — -  ,  x  >  t,  9,  p  >  0,  — °°  <  t  <  °°. 

e^rcr,) 

We  may  denote  this  p.d.f.  by  GAM(  9,  p,  t)  .  The  shape  parameter  is  p, 
the  scale  parameter  9,  and  the  location  parameter  t.  For  i  =  0  or 
known,  this  reduces  to  the  two-parameter  gamma  distribution  where 

_x 

9  p-1 

f  (x;  9,  p)  =  - — — -  ,  x  >  0,  p,  9  >  0. 

9  T(p) 

For  p  =  1  we  have  the  exponential  distribution  with  parameter  9  and 
for  integer  p  we  have  the  Erlang  distribution.  For  p  <_  1  and  fixed 
9,  the  p.d.f.  is  decreasing  in  x  and  unbounded  (p<l)  near  the 

T 

2  r  ©  tv  tH“1t  ") 

origin.  The  mean  is  9p,  variance  is  9  p,  E(X  )  =  — f(70 an<^ 

— p  1 

moment-generating  function  is  M  (t)  =  (l-9t)  ,  t  <  -^  .  The  shape 

parameter  p  is  the  reciprocal  of  the  squared  coefficient  of 


variation. 
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o 

ii 


<o 


to 

toC 

•H 

t-u 


X 

4-1 


Gamma  p.d.f.  for  various  values  of 
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If  T  denotes  the  waiting  time  until  the  q  occurrence  in  a 
Poisson  process  with  intensity  then  T  ~  GAM( 9  =  -7  ,  9,0).  The 

waiting  time  until  the  first  occurrence  may  be  denoted  as 
GAM( 9  =  ,1,0)  which  is  EXP (9  =  y  ) .  Thus  the  gamma  distribution  is 

a  natural  extension  of  the  exponential  distribution.  The  two-parameter 

exponential  may  be  denoted  as  GAM( 9  =  \  ,  q  =  1 ,  x)  where  x  is  a 

A. 

location  parameter. 

n 

If  X1,...,Xn  i.i.d.  GAM(9,mi,0)  then  Y  =  l  X± 
n 

~  GAM( 9,  y  m.,0).  Consider  the  case  of  a  component  with  n-1  spare 
i=l  1 

parts.  If  X_^  ,  i  =  l,2,...,n  denote  the  lifetimes  of  the  component 

and  the  spares  and  if  each  is  distributed  GAM( 9,1,0)  =  EXP (9  =  — )  , 

A. 

then  the  lifetime  of  the  system,  assuming  use  of  the  n-1  spares,  is 

v  2Y  2 

Y  =  }  X.  ~  GAM(9,n,0)  or  -5—  ~  .  Note  also  that  if 

.  1  9  A,2ndf 

i=l 

2 

Y  ~  GAM(2,q,0)  then  ^  ~  ^2q  df  *  hazard  function  (HF) 

h(x)  =  X  =  ^  for  9=1;  h(x)  ■>  \  as  x  -►  ®  for  9  <  1  and 

h(0)  =  0;  h(x)  ->•  x"1"  as  x  ■*  and  h(0)  =  ®  for  9  >  1. 

Consequently  the  gamma  distribution  can  model  systems  in  a  regular 
maintenance  program  where  the  failure  rate  is  likely  to  increase 
initially  but  then  stabilize. 

Outliers  in  gamma  samples  arise  in  any  context  where  Poisson 

processes  are  appropriate  basic  models,  e.g.  traffic  flow,  biological 

aggregation,  failure  of  electronic  equipment.  They  also  occur  in  the 

context  of  a  shifted  exponential  or  gamma  distribution.  Outliers  in 
2 

X  samples  arise  in  ANOVA;  outliers  in  gamma  samples  of  arbitrary 
shape  parameter  arise  with  skew-distributed  data,  for  which  the  gamma 


distribution  is  often  a  useful  model. 
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3.2  Outlier-proneness  of  the  exchangeable  model  with  the  Gamma 
distribution. 

Neyman  and  Scott  (1971)  showed  that  for  i.i.d.r.v.’s  the  family  of 
gamma  distributions  indexed  by  shape  parameter  p  is  outlier-prone 
completely  on  the  right.  Let  us  consider  now  the  exchangeable  model. 

CASE  I:  Scale  change 

Kale  (1975b)  showed  that  the  exchangeable  model  with  (at  most)  one 
spurious  observation  for  scale  parameter  families  for  non-negative 
random  variables  involving  a  possible  change  in  scale  is  outlier-prone 
completely.  This  would  apply  to  the  exchangeable  model  involving  the 
gamma  distribution  with  possible  change  in  scale  parameter  9. 

CASE  II:  Shape  change 

Now  consider  the  exchangeable  model  with  (at  most)  one  spurious 
observation  based  on  the  gamma  distribution  with  possible  change  in  the 
shape  parameter  p.  Without  loss  of  generality  (w.l.o.g)  we  may  take 
9=1,  t  =  0. 

Then  n-1  observations  have  p.d.f.  f(x;_9)  which  is  GAM(l,p,0) 
and  one  observation  has  p.d.f.  f(x;£)  which  is  GAM(l,k*p,0)  where 
k*  >  1.  We  shall  show  that  this  model  is  outlier-prone  completely  on 
the  right. 

Now  we  may  write  the  likelihood  as 


x .  >0,  k*  >  1 ,  p>0 

l  —  — 


(k*>l  since  we  are  considering  x 


(n) 


as  a  possible  outlier).  Letting 
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E  denote  the  event  that  x(n)  1 s  a  (k,n)-outlier  on  the  right  and 

P(k,n|L)  =  P(E|n,k*),  we  must  show  that  sup  P  (E  |  r| , k* )  =  1  in  order 

r|>0 

k*>l 

for  this  family  to  be  (k,n)-outlier-prone . 

The  joint  p.d.f.  of  the  order  statistics  x,,  ...... x,  .  is  given 

(1)  (n)  6 

by 


8(xq)»  •  •  •  ,x(n) ’  h>k*) 


n!L(X(l)» *  *  *  »x(n) ; h»k* ) 


=  (n-1 ) !  I  n  f(x.  .;p)f(x,  .;k*p), 
r=l  i^r  K  J  K  } 


0  <  x 


(1) 


< 


<  X.  v  < 

(n) 


Thus 


P (E |  p,k*)  =  J  gCx^y  •>x(n)^»k*)dx(i)‘-*dx(n) 

E 


n 


I  i  (n-1)!  IT  f  (x^)  ’  r,)f  (X(r)  ;k*p)dx(.1.  ...  dx,  . 
r=l  E  i^r  K  K  ^  K 


n 


-  l  i  (ti.k*) 

r=l 


kx  .+X 

and  E  is  such  that  0  <  x  >  -  .  <  x,~N  <  ...  <  x,  <  - r~n - 

(1)  (2)  (n-1)  k+1 

<  X  x  <  °°. 

(n) 
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Lemma  3.2.1:  0  _<  Ir("n,k*)  <_  u(r;n,r),k*)  where  u(r;n,r|,k*)  =  P{X^^ 

is  the  "spurious"  observation  whose  p.d.f.  is  f(x,k*Ti)}. 


Proof : 

u(r;n,p,k*)  =  (n— 1 )  !  J  II  f  (x  ;  p)f  (x,-,k*p)dx dx  where 

s  Ci)  (r)  (1)  (n) 

S  is  such  that  0  <  x,..  N  <  x,_N  <  . . .  <  x,  ,  N  <  x,  N  <  °°. 

(1)  (2)  (n-1)  (n) 

Now  I  (ruk*)  =  ( n— 1 )  !  [  II  f  (x, ,  p)f  (x,  N  ;k*ri)dx/,  N ,  .  .  .dx,  N  where 

t  El±r  ^  (r)  C1)  <n) 

kx^.+x 

E  is  such  that  0  <  tl,.  n  <  x/0.  <  ...  <  x,  <  - m - 

(1)  (2)  (n-1)  k+1 


<  x  v  <  °°* 

(n) 


Since  E  c  S,  0  <_  I  ( Ti,k*)  <^u(r;n,p,k*),  r  =  1,2,  . . .  ,n. 


Theorem  3.2.2:  lim  u(r;n,p,k*)  = 


1  , 

r  = 

n 

0  , 

r  = 

n-1 

Proof:  Now  u(r;n,p,k*)  =  ("  J)]  (F(y;p)}r  1  {1-F(y;  r))  }n  rf(y;k*p)dy 

0 

y 

where  F(y;r|)  =  j  f(x;p)dx. 

—03 

Consider  now 


u 


00  y  -x  n-1  ,  -y  k*n-l 

/  ,j.n  f  ft  ex  ,i  n-1  e  y 

(n;n,p,k*)  =  J  {}  —  dxJ  - Z - 

y=0  x=0  ^ 


T(k*p) 


dy,  k*  >  1,  r|  >  0 


(without  loss  of  generality  0=1,  x  =  0) 


Let 
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00  y  -x  rr-1  .  _ 

I'(Tl,k*)  =  j  {J  6  -3_v-  dx}n  e  7yVdy  for  v  =  k*T)-l 
y=0  0  ^  \  M/ 


Now  setting  y  =  vt,  dy  =  vdt 


-  I  (J 


00  vt  -x  n-1  .. 

ex  n-1  -  vt ,  N  v 


t=0  0 


r(T)) 


.  -I  n-i  -  vr ,  s  v  , 
dx }  e  (  vt )  vd  t 


v+1  <  Q,  .  v(Xn  t-t). 

=  v  J  iHv,t)e  dt 

t=0 


where 


tf(v,t) 


-x  n-1 
e  x 


r(n) 


and  where  there  is  no  loss  in  generality  in  assuming  v  >  0  (i.e. 

k*ri  >  1)  since  we  are  interested  in  lim  I f  (  r),k* )  . 

k*+~  n 

For  0  <_  t  <_  a,  i?(v,0)  <_^(v,t)  _<i>(v,a)  and  0  _<  <?(  v,t)  <_  1  for  all 
v,t  and  $(v,t)  ■»  1  a.e.  in  t  as  v  ■*  00  (i.e.  as  k*  ®).  We  may 

00  \ 

.  1  0/  N  v(in  t-t),^  ,  .  v(  In  t-t), 

write  J  i?(v,t)e  dt  =  J  $(v,t)e  dt 

0  0 

CO 

I  n  /  N  v(  Jin  t-t )  , 

+  J  tf(v,t)e  dt. 

1 

For  each  of  the  integrals  on  the  right  hand  side,  the  dominant  part 
occurs  in  the  neighborhood  of  the  point  when  Jin  t-t  is  maximum,  i.e. 
t  =  1.  Following  Copson  (1967  ,  p.  36  ff),  since  0  <_  t?(v,t)  _<  1  for 
all  v  and  t  and  $(v,t)  is  continuous  and  increasing  in  t  there 
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is  a  number  tQ(0<tQ<l)  such  that 


/1«v>t)ev(to  ^dt  -  Wv.t.)  | 1  ev(,!n  C'C)dt 
0  0 


0  <  co  <  1 


We  cannot  have  t^  =  0  since  $(v,0)  =  0  while 


j  i?(v,t)eV(An  C  t)dt  >  j  l?(v,t)eV(int  t}dt 
0  .5 


,  -v(-y  -in  2) 

2)e  >  ° 


and 


j  $(v,t)eV^n  C  t'^dt  =  lim  J  i)(v,t)eV^n  t'^dt 
1  1 


a 


=  lim  ?Kv,t1)  J  e  V(ln  C  t)dt, 
a-*30  1 


1  <  h 


<  a. 


Thus 


j“«,,t)eV(ln  t_t)dt  ~  j"  ev(ln  t_t)dt 
0  0 


1 


r(vhi) 

v+l 

V 


as 


V  ->  °o 


/■N/ 
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(see  Appendix  II).  As  a  result,  we  may  write 


I ’  (r|,k*)  ~  v^-l  =  r(vt-l)  as 

n  v+1 

v 


V  <=° 


and 


r(p,k*) 

u(n;n,ri,k*)  =  ■  where  v  =  k*p-l 


•>  1  as  v  -*  00  (i.e.  as  k*  ■>  co) . 


For  u(r;n,p,k*),  r  =  I,2,...,n-1  we  have 


0 


_<  u(r;n,p,k*)  =  )j  {J  ~p7'“T —  <*x}r  l 

y=0  0 


y  -x  ri — 1  -y  k*n-l 

fi  i  e  x  j  i  n-r  e  J  y 
{1  "  1  dx}  ~~~~  dy 


o 


r(k*p) 


.  i 1 ( n>k*) 

n-l>  r  * 

T(k*p) 


=  r  ) 

ir-i  > 


where 


i;(n,k*) 


J  (J 

t=0  0 


Vt  -X  r|-l 


e  x 


dx}r  1  {l  -  J 


r(n) 


Vt  -x  r|— 1 


e  x 


0 


T(p) 


,  in-r  -vt  \H-1  ' 
dx  [  e  v  t 


dt , 


v  =  k*p-l,  y  =  vt,  dy  =  vdt  and  again  we  assume  v  >  0.  Therefore 
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i;(n,k*) 


V+l  / 

v  j 


tK  v,t)e 


v(ln  t-t)dt 


v+l  ,a  .  N  v(ln  t-t), 

=  v  lim  J  !?(v,t)e  dt  . 

a-x»  0 

vt  -x  r|-l  vt  -x  r| —  1 

Now  i?(v,t)  =  {J^  -YO)j —  dX^  f1  "  i  ~e~f(>j) —  dx}n  r  1  1  and  as 

k*  -+  00 ,  v  ->  00  and  $(v,t)  ->>0  almost  everywhere  on  [0,°°].  Since 
i?(v,t)  is  continuous  and  bounded  for  all  t 


ia  Q/  .  N  v(ln  t-t),.  Q/  i  f ^  v(ln  t-t)  , 

J  tf(v,t)e  dt  =^(v,t  )J  e  'dt 

0  U 


a 

0 


0  <_  t^  _<  a 


<*(v,t0)j  evan  t-t)dt ,  (,(V)tn)  . 


0 


O'  v+l 
v 


("  !)  «v,t 0)r(v+i) 

Therefore  0  <  u(r;n,r|,k*)  <  -  v  lim  r 


r(k*ri) 


(X+co  vfl 

V 


=  fn~L ) 
lr-l  J 


lim  ??(v,tQ).  But  (^_d)  lim  $(v,tQ)  >0  as  v  -*  °°  (i.e.  as  k*  -*  <=) 

a-*°  a-xn 

and  hence  u(r;n,-n,k*)  ->  0  as  v  -*  °°  for  r  =  l,2,...,n-l. 


Theorem  3.2.3:  The  exchangeable  model  with  (at  most)  one  spurious 
observation  based  on  the  gamma  distribution  with  a  change  in  shape 
parameter  is  outlier-prone  completely  on  the  right. 

Proof:  We  need  to  show  sup  P(E|p,k*)  =  1. 

T)>0 

k*>l 
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n 

Since  P (E  |  r),k*)  =  l  I  ( t),k*)  , 

r=l 


n 

lim  P  (E  |  r|,k*)  =  l  litn  I  ( r),k*)  . 
k*-Hu  i=l  k*-*° 

From  Lemma  3.2.1,  0  <_  I  (r|,k*)  _<u(r;n,r|,k*),  r  =  1,2, ...  ,n. 

0  ,  r  1,2,..., n- 1 

By  Theorem  3.2.2,  lim  u(r;n,r|,k*)  = 

k*-Ko 

1  ,  r  =  n 

Thus  lim  I  (ri,k*)  =0,  r  =  l,2,...,n-l  . 
k*+~  r 

and  lim  P(E|p,k*)  =  lim  I  (p,k*). 

k* ->-ao  k*  ~^co  ^ 

Now  we  need  only  show  that  lim  I  (r|,k*)  =1  in  order  for  this  model 


k*-Ko 


n 


to  be  outlier-prone  completely  on  the  right. 


But 


I  (T|,k*)  =  (n-1 ) !  ]  IT  f(x.  ,;r|)f(x,  ,;k*r))dxn.  . 
n  E  i#n  (1)  (n)  (1) 


.  dx 


(n) 


X(i)  r)— 1  X(n)  k*ri-l 

=  (n’1)!^E  irt  r<Ti)(  >  r(k*<n))  dx(l)  •••  dx(n) 


where'  E  is  such  that  0  <  x^,^  <  X(2)  ^  * 


kx,  .+x,  N 

<  x  <  (X) 

*  (n-1)  k+1 

<  x  <  °°. 

(n) 


Let 


n 


denote  the  order  statistics  x 


(1)’ 


,X(n)’ 


respectively.  Then  we  may  write 


In(T1,k*)  =  (n_1) 


JJ  IJ 

0<y  <y  <®  0 

I  n 


k+1 


-x  ri-l 
e  x 


r(n) 


d  x  -  J 


0 


e  XxTri  ,  i n-2 
dx} 


71  T) —  1  yn  k*rrl 

e  yl  e  yn 

r(  -n)  r(k*r))  dyidyn 


00  yn  ky  +y 

=  J  [J  (n-l)  |F(--k--j'-n  ;Tl)-F(y1;ri)}n  f(y1;ri)dy]L] 

yn=o  yx=o 

f (y  ;k*T])dy 
n  n 


-y  r(— 1 

where  f  (y ;  ri)  =  e-  ,  y  >  0.  Then 

oo  — y 

I  (ri,k*)  =  -  ,  v  J  i?(y  )e  nyk  11-1  dy  where 

n  ’  r(k*ri)  J  _q  n  'n  yn 

yn 

y 

n  ky,  +y  _o 

i>(yn)  =  j  ( n-l )  {f  (  k+1  n  ;  Ti)-F(y1;  Ti)  }n  £(y1;n)dy1.  Setting 

y1=0 

y  =  vt  where  v  =  k*r|-l,  we  obtain 


(3.2.1) 


In(T1,k*) 


v+1 

V 

r(  v+i) 


/”  ,?<v,t)ev(to  t_t)dt 
t=0 


v+1 

V 

r(  v+i) 


lim  j  l?(v,t)eV(ln  C  C)dt 

(X^oo  t  =  0 


where  t?(v,t)  =  t?(vt) 


vt  ky.  +  vt 

j  (n"l)lF("^i —  ;  h)“F(y1;  ti)  }n  f(y1;ri)dy 
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As  k*  ->  °°,  v  ■*  00  and  the  major  contribution  to  this  integral  in 
3.2.1  occurs  in  the  neighborhood  of  t  =  1.  Also  0  ^$(v,t)  _<  1  for 
all  v,t  and  $(v,t)  ->  1  a.e.  on  [0,®)  as  v  ■*  00  (i.e.  as  k*  -*  °°), 

since 

ky^+vt 

lim  [F(-g^- —  ;ri)-F(y1;ri)]  =  1-F(y1;  n)  a.e.  on  [0,°=). 
v-*=° 

vt  _  ky..  +  vt 

If  t  =  0,  $(v,t)  =  (n-l)J  [F("kTl —  5  rl)-F(y1;ri)  ]n  f(y1;Ti)dy1  and 

$(v,0)  =  0.  But  #(v,t)  is  nonnegative,  bounded,  increasing  and 

a  ,  . 

absolutely  continuous  in  t.  Thus  2t~  e  [0,°°)3  .  J  i?(v,t)e  'dt 

0 

a 

=  t?(v,t„)  j  ev^  n  t  dt.  On  the  other  hand  for  any  0  <  E,  <  1 
U  0 

,a  Q,  N  v(£n  t-t).  .  Q,  N  v(Jln  t-t), 

J  tf(v,t)e  dt  >  J  i>(v,t)e  v  'dt 

0  5 

>  r?(v,e)eV(^1  ^“^(l-£) 

>  0. 


Therefore  t^  ±  0. 


_ _  ,a  n.  .  v(ln  t-t),  (  v(£n  t-t),.  r(v+l) 

Thus  J  tf(v,t)e  dt  ~  J  e  'dt  =  --  ■— —■ 


0 


0 


as  v  -*  °° 
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and 


in(n,k*) 


v+1 

v 

r(  vfi) 


a  ,  „  . 

,  .  I  q/  v(  Jin  t-t ) 
lim  J  tf(v,t)e  dt 

a+®  0 


V+1 

V 

r(  v+i) 


r(  vfi) 
v+1 

V 


as 


v  ->  1:0  (i.e.  as  k*  ->  co)  . 


Since  lim  I  (p,k*)  =  1  and  lim  I  (ri,k*)  =  0,  r  =  l,2,...,n-l  and 
k*-*=°  n  k*-**>  r 

n 


P(E|n,k*)  =  l  I  ( ri,k*) , 


r=l 


sup  P(E  |  r|,k*)  =  1 

p>0 

k*>l 


and  the  family  is  outlier-prone  completely  on  the  right. 


Thus  for  either  a  scale  or  shape  change,  the  exchangeable  model 
with  at  most  one  spurious  observation  based  on  the  gamma  distribution 
is  outlier-prone  completely  on  the  right. 
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3.3  Detection  of  Outliers 

Mount  and  Kale  (1973)  considered  a  general  model,  assuming 

X. , . . . ,X  are  such  that  n-1  of  them  are  distributed  with 
1  n 

distribution  function  (d.f)  F(x)  and  one  is  distributed  with  d.f. 

G(x),  where  F  and  G  are  stochastically  ordered,  i.e.  G  <  F.  A 

priori,  each  X^  has  probability  -i-  of  being  the  spurious  observation 

dG 

distributed  as  G.  Let  Y(x)  =  —  •  If  Y(x)  is  monotone  increasing, 
then  u(l;n,k*)  <  u(2;n,k*)  <  ...  <  u(n;n,k*)  where  u(i;n,k*)  =  P(X^^ 
is  the  spurious  observation  in  a  sample  of  size  n  with  k*  the 
coefficient  of  spuriosity). 


CASE  I:  Scale  change:  Consider  a  situation  where  n-1  observations 
are  from  GAM(9,"n,0)  and  one  is  from  GAM(k* 9,  r),0) ,  k*  >0.  If  k*  =  1 
we  have  homogeneous  data.  Now 


Y(x) 


dG(x) 

dF(x) 


x  /-1-k*  \ 
9^  k*  j 


k* 


h 


and 


(x) 


X/'l-k*'! 


k*11 


k*-l 

k* 


>  0  if  k*  >  1 

<  0  if  k*  <  r 


Thus,  if  k*  >  1 ,  Y(x)  is  monotone  increasing  and 
probability  of  being  spurious;  if  k*  <  1 ,  H  is  monotone 


and  X.-j.  has  maximum  probability  of  being  spurious. 


has  maximum 
decreasing 
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CASE  II:  Shape  change:  If  n-1  observations  come  from  GAM(9,p,0) 
and  one  is  from  GAM(  9,k*r),  0) ,  k*  >  0,  assuming  the  exchangeable 
model , 


nx) 


dG 

dF 


0k*'rr|r(k*Tl) 


and 


^t(x)  =  r(ri)Ti(k*-l)x 


k*r|-r|-l  l  >  0  if  k*  >  1 


T(k*n)9 


k*TT-r) 


<  0  if  k*  <  1  . 


Thus,  if  k*  >  1 ,  Y(x)  is  monotone  increasing  and  X.  N  has  maximum 

(n) 

probability  of  being  spurious;  for  0  <  k*  <  1,  Y(x)  and  u(r;n,k*) 
are  monotone  decreasing  and  has  niaximum  probabilty  of  being 

spurious . 

This  now  shows  that  the  spurious  observation  resulting  from  a 
scale  or  shape  change  is  most  likely  to  occur  at  the  sample  extremes 


i.e.  it  tends  to  show  up  as  an  outlier. 
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3.4  Estimation  for  Gamma  parameters 
3.4.1  Standard  Estimators 

The  gamma  distribution  is  a  member  of  the  exponential  class, 
n  n 


consequently  (  £  X  ,  £  Jin  X  )  are  complete  sufficient  statistics 

i=l  1  i=l  1 

for  (9,q). 

~  In 

For  known  shape  q,  the  MLE  of  X  =  -q  is  3  which  is 

x 


biased  but  consistent  and  asymptotically 

9  (q  known)  is  —  . 

For  unknown  q,  we  obtain 


N(>,— ). 
v  nq J 


Thus  the  UMVUE  for 


MLE 


x 


h 


W  =  In  W  '  +  to  ;  =  0 


n 


where  x  =  [  II  x.  )  and  cJj(z)  is  Euler’s  psi  functi 


ion  i.e. 


i=l 


e  t  -zt 

-  e 


(Kz)  =  -j—  Jin  T(z)  =  J  — - — -  dt  =  -y  +  —  +  z  £  [i(i+z)]  ^  where 

az  A  *  z  ,  4 


0  l-e 


i=l 


y  =  .5772157  (Euler's  constant).  These  equations  must  be  solved 
iteratively  (see  Choi  and  Wette  (1969)).  For  large  q^^ ,  Linhart 
(1965)  suggested  approximating  Jin  q^g  -  (Rq^g)  by  (2hMLE“l/3) 


thus 
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=  {(Jin  x- Jin  x)  1+l/3}l/2. 

Moment  estimators  give 


n  9 

(  l  xt)2 

i=l 

^  =  — - : 

l  n(X.-X) 


These  have  lower  asymptotic  efficiency  than  MLE’s  but  Lilliefors  (1971) 
showed  that  for  n  _<  20,  p  >_  2,  the  MSE’s  of  the  moment  estimators  are 
close  to  those  of  the  MLE’s.  Both  the  MLE’s  and  the  moment  estimators 
are  biased.  To  approximate  zero  bias  and  smaller  MSE,  Lilliefors 
suggested  the  following  corrections: 


TABLE  3.4.1  Lilliefors’  improved  estimators  of  r\  and  \  =  1/9 


Based  on 

M.L.E.  Moment  Estimators 


A 

A 

^MLE 

f) 

1+3 /n 

c 

1+2 /n  5/3 

/\ 

A 

1 

n,1MLE 

(  -5/3}  1 

e 

7  11 

l(l+2/n)  '  }  n 

(1+  l 

i=l 

l  Xi 

i=l 
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Thom  (1968)  has  given  estimators  very  similar  to  the  MLE’s  for 


T)  >  1  • 


1  + 


= 


VI  -L  4M 
1  +  ~ 


4M 


0*  = 


where  M  =  ln{x/x)  =  -In  and  is  sufficient  for  p  and 

independent  of  9.  Bain  and  Englehardt  (1975)  have  a  chi-square 
approximation  for  M  valid  for  all  p  ( 9  acts  as  a  nuisance 
parameter) . 


2nr)Mc  ~  y 


vdf 


where  W  =  2npM,  c  =  2 


E(W)  . 


Var(W)  nw2(q)-w2(np) 


and  v  =  2 


[E(W)]  = 

Var(W) 


[nw^(  p)-w^(np)  ]  c  and  w^(z)  =  2z  {in  z—  cl>( z )  } ,  w2(z)  =  2z  (zcl*  (z)-l  } 

d  2 

where  <Kz)  is  the  psi  function.  As  p  -*  0,  W  -*  X2(n-l)df* 


W1(T1)  1  +  1+6  p 


w2(p)  = 


1  + 


1+2 . 5p 


1  +  x— 

3p 


,  0  <  p  <  2 


,  p  >  2 


v 


=  1  + 


( 1+4.3  p)' 


n-1 
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3.4.2  Estimators  suggested  for  use 


Case  I:  Scale  Change 

Assuming  the  exchangeable  model  where  n-1  observations  are 

distributed  as  GAM(9,r|,0)  and  one  is  distributed  as  GAM(  9k*,  rj, 0 )  ,  r) 

known,  k*  _>  1,  we  have  for  k*  =  1,  homogeneous  data  and  the  best 

n 


l  Xi 

i=l 

linear  unbiased  estimator  (BLUE)  for  9  is  9,„_  =  - -  and 

MLE  nr) 

E(W  -  0  and  MSE<W  "T"<W  +  {Blas(W)2  -  hi-  FOr 

^  1^1  1 

the  heterogeneous  case  (k*>l),  E,  (9.__)  =  —  Y  f— —  9rH - k*9n] 

het  MLE  nr)  v  n  n  ; 

9  - 

=  — (n-l+k*)  and  Bias  ..(Own,)  =  -  •  As  k*  >  this  bias  ■>  ®. 

n  het  MLE  n  ’ 

Now  MSEhet(0MLE>  ■  Varhat<W  +  {Biashet(  W}2  and 


Varhet(9MLE) 


(nri)  i=l 


l  02 ^  +  |  (k*9)Zri} 


-  -y-  {k*2+n-l}. 
n  T| 


therefore 


MEEhet  ^  9MLE  ^  * 


9 


2 

n  r) 


(k*  +n-l )  + 


2  2 
9  (k*-l) 


n 


=  ~y~  {k*2+n-l+ri(k*-l)2} 
n  T) 
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n 

l  xi 

-  *  i=l 

As  k*  ->  MSE,  (9.__)  00  and  hence  0.  _  _  =  -  proves  to  be  a 

het  MLE  MLE  np 

poor  estimator  of  9. 

If  we  consider  a  random  sample  of  size  n  where  n-m 
observations  are  from  GAM(9,p,0)  and  m  are  from  GAM(k* 9, p, 0) ,  k*  >  1 
(i.e.  same  shape  but  a  change  in  scale  parameter)  and  a  priori  each 
subset  of  m  observations  is  equally  likely  to  be  the  "outlier 
subset",  then  the  likelihood 


l(x|  e,p,k*,n  = 


H 


-  i  !i 

xil  9 

1  rr  h“l 

e  II  x. 

x.  l\  1 

_ 1 _ 

(r(p)eT1}n"in 


-  I 


.k*  9 


x.  el  . 

e  1  n  xT1 

V1  1 

{r(n)(k*e)r|}ln 


X.  X. 

_r  y  -A+  T  — 1 

L  9  %  91 

x.el  x. el  1 

n  ,  1  1 

H  xj  e 
1  i=l  1 


n 

vm  ' 


{r(n)}neTl(n-m)0j|m 


for  k*9  =  9  where  I  =  (x.  ,x.  ,...,x.  )  and  I  e  f,  the  collection 

X1  X2  Xm 

of  all  possible  combinations  of  m  observations  out  of  n.  (I 

represents  the  subset  of  m  spurious  observations).  To  maximize 

L(xJ  9,p,k*,I)  for  9,p  >  0,  I  e  f,  k*  >  1  we  use  the  fact  that 

T(x)  =  is  monotone  increasing  for  k*  >  1  and  hence 

dr  ix; 

max  L(x| 9, p,k*,I)  occurs  at  1=1=  (x.  11N,...,x/  .)•  Therefore 

Vjr-  — 1  *  ”  ’  (n-m+1)  (n) 
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n-m  x 


n 


-(  I 


(i) 


n  x 


(i) 


L(x|  9,  ri,k*,I)  = 


1 

(!) 


T|—  1  i=l 
n  x  e 

i=l  1 


i=n-m+l 


vm 


{ro1)}neTl(n~m)01rp 


and 


K(xJ  9,  r|,k*,i}  =  In  L(xJ  9,  r|,k*  ,1) 


n 

=  C  +  ( r|-l )  l  Jbi  x 
i=l 


n-m  x 


(i) 


i=l 


I 


i=n-m+l 


-  nAn  r(  r|)  -  q(n-m)Jln  9  -  qmJhi  9^  . 


Assuming 


i)  9,  9^  and  r)  are  unknown,  we  obtain: 


5K 

5r| 


n 


=  l  ^  x  - 

i=l 


n^- —  {(n-m)  In  9  +  mJbi  9.  } 
r(T])  1 


5K  _  nym  X(i)  _  r|(n-m) 
59  "  L 


i=l  9 


9 


5K 

59, 


n  x 


(i) 


i=n-m+l  ^1 


r|m 

97 


Setting  the  above  three  equations  equal  to  zero,  we  obtain 
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n-m 


9  = 


1  x/  •  \ 
i-i  (1) 


ri(n-m) 


n 

I  X 

i  _  i=n-m+l 

1 ' — ; — r 

ri(n-m) 


(i) 


n 


h(ri)  =  l  in  x(i)  -  n-F-^-vy  +  nAn  T) 


-  (n-m)  An 


n-m 

l  xa) 

i=l  V  ' 


n-m 


-mAn 


n 

I  X 

i=n-m+l 


(i) 


m 


=  0 


ii)  for  the  case  of  known  r|,  we  obtain 


n-m 

I  x 

i=l 


(i) 


T)(n-m) 


A 

and  9^ 


n 


I 

i=n-m+l 


x 


(i) 


mr) 


which  are  trimmed  means. 


Xi 

iii)  for  the  case  of  known  9,  we  may  use  Y.  =  —  . 

1  U 


n 


L(y;  r),k*  ,1)  = 


n  h-l  i A 
n  y.  e 

1  i=l  L 


-{  l  y±  +  l  yJ^*} 

J  i  61 


(")  {r(r,)}"k*rp 

m 
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and 


n-m 


n 


L(y;  ri,k*,i)  = 


"  ,-l  '(Jl  y^>+  i=n-m+l  y^>/k<1 

n  y '  e 

1  i=l 


n 


(r(  ri)  >\*Tlm 


and 


K(y;ri,k*,I)  =  JinL(y;  ri,k*,I ) 


n 


n-m 


n 


=  C'  +  ( r)-l )  l  hi  y  -  {  l  y  +  T  y 

i=l  i=l  k  }  i=n-m+l  K  ‘ 


-  n  Jin  TCp)  -  rp  In  k*. 


Then 


5K 

9ri 


n 


i=i 


In  y. 
1 


nr’Cn) 

r(h) 


m  Jin  k* 


and  thus 


and 


5K 

5k* 


l 


i=n-m+l 


mr| 

k* 


k* 


n 


i=n-m+l 


(i) 


mt] 


nr’(p) 

r(h) 


n 

-  m  Jin  T)  =  1 

i=l 


yn  y± 


m  Jdi 


( i=n-m+l ^ ^ 1 ^ 

l  ” 


/k*} 
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iv)  for  the  case  of  known  k* 


5K 


I 


i=i 


In  x. 
1 


nrf(ri) 
r(  r,) 


(n-m)  in  9  -  m  J?n(k*9) 


,  n-m  X(i)  n  X(i) 

=  y  _ +  y  _ 

50  4  2  L  2 

i=i  9  i=n-m+l  k*9 


Tl(n-m) 

9 


pm 

9~ 


Setting 


5K 

5r) 


0 


and 


5K 

59 


0,  we  obtain 


n-m 


n9  = 


i=l  k*  i=n-m+l  (l) 


n 

l 


n 


In  9  + 


rctj) 

r(fi) 


n 

y  An  x.  -  m  Jin  k* 
i=l  1 
n 


Case  II.  Shape  change 

Assuming  n-1  observations  distributed  as  GAM(9,p,0)  and  one 
distributed  as  GAM(  9,k*r|, 0) ,  k*  >  1  ,  we  have,  for  heterogeneous  data 


(k*>l), 


E,  ( 9  )  =  E  (— )  =  —  l  9r)  +  —  9k*  n) 

het  MLE  '“n '  nr|  .  ,  n  n 

1=1 


=  9  + 


9(k*-l ) 


n 
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As  k*  -*  00 ,  bias  >  ro.  Also 


Varhet('  9MLE 


)  =  Var(^)  = 

v  Tp 


2  2 
n  ri 


n 

I 

i=l 


Var(X.) 


1 


2  2 
n  r| 


n 


I  ( 

i=l 


n-i  q2  ,  1 
— —  9  ti  +  - 
n  n 


e2k*ri) 


Then 


=  (n+k*-l). 

n  ri 


MSEhet(W  ' 


Var,  (©„__) 
bet  MLEy 


+  {Bias 


(9MLE^  * 


02(n+k*-l)  02p(k*-l)2 

2  2 
n  r)  r|n 


e2  {n+k*-l+p(k*-l)2} 

2 

n  r) 


As  k*  -*  MSE^et^  °°  anc*  ^ence  ®  =  ~  proves  to  be  a  poor 

estimator  of  9  (We  have  assumed  p  known) . 

Consider  now  the  exchangeable  model  where  n-m  observations  are 
from  GAM(9,p,0)  and  m  observations  are  from  GAM(  9,k*p,  0) ,  k*  >  1 
(i.e.  possible  change  in  shape;  scale  constant).  Then 
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L(x|  9,  p,k*,I)  = 


1  e 


-  I 

x.  a 

i 


X. 


e 


n  x 


Tl-l 


-  y  — 

4  ft 


x_^  il 


X.  el 

l 


n  k*p-l 

n  x. 

x  el  1 
l 


(n)  (r( -n) }n  m(0Ti)n"m  (r(k*p)}n 


-T/e 


n  x 


rrl 


n  x 


k*p-l 


x.  il 
l 


x .  el 
l 


(n)  {r(r1)}n"m{r(k*r1)}m(9T1)n  m+k*m 

m 


n 


where  T  =  \  x.  ,  I  =  (x.  ,  x  , ...,x,  }  and  I  6  f  ,  the  collection 

i=l  1  1  2 


m 

of  all  possible  subsets  of  m  observations  out  of  n.  To  obtain 

dG 

MLE's,  note  that  Y(x)  =  -ttt  is  monotone  increasing  and,  by  Kale 

Qr 

(1975b),  we  know  max  L(x| 9, p,k*,I)  occurs  at  f  =  I(9,r|,k*  fixed) 

lej 


where  I  =  (x,  ,x,  N ) .  Thus 

(n-m+1)  (n) 


max  L(xJ  9, p,k*,I)  =  max  L(xJ  9,  p,k*,I) 
k*>l  “  k*>l 

Ief 


and 


L(x|  9,  p,k*,I)  = 


/  a  n*  m  -  tt  .  .  1 

-T/9  rrl  T,  k*p-l 

e  n  x  '  n  x, . N 

_ 1-1  (1)  l=n-m+l  (1) 

(n)  {r(,)}n-m{r(k^)}m(05n'm+k*m 


Using  p^  =  k*p  (  >  p  since  k*  >  1 ) 


K(x|  9,  p,  ti1,I)  =  In  L(x|  9,  p,  p1?f  ) 


n-m 

=  c  -  -g  +  (p-1)  I  An 


n 


1=1  +  ( v15±  Jm+1  ta  xd) 
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-  (n-m) In  r(  r|)  -  m  &i  T(  )  -  (n-m) p  Jin  9  -  mr^Jln  9 

and 

i)  assuming  9,  p,  and  p^  unknown,  we  obtain 


5K  _  T  (n-m)Ti  mT1l 
59  Q2  9  9 


n-m 


8K  =  “Y“  9n  (n-m)r’(p) 

5t|  isfL  1  (i)  r(n)  (n  m)  Jin  9 


5K 

5 


n  mr’Cr^) 

!  V  in  x,,.  -  —=r, - r - m  Jin  9 

P  i=n-m+l  (1)  F(V 


5K  -  T 

Now  =  0  implies  9  =  -  and 


A  A 


(n-m)  r)  +  mru 


n-m 


l  in  x 

*  =  0  implies  1^1  +  in  0  =  i=i - — 

an  r(p) 


and 


n 


5K  r'(V 

— —  =  0  implies  -  +  Jin  9  = 

5p, 


l  xm 

i=n-m+l  ^  ; 


m 


ii)  for  the  case  of  known  shape  parameter  p, 


5K  _  T  (n-m)p  mk*p 
59  “  02  9  9 


and 
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5K 

9k* 


n 

T]  l 

i=n-m+l 


In  x 


(i) 


mn 


rf  (k*n) 

r(k*ri) 


-  mr)  Jin  9 


and  setting  these  equal  to  zero,  we  obtain 


e  = 


q(n-m+mk*) 


and 


n 


I  In 

a  r’(k*r)  a  i=n-m+l 

«k*>  ■  +  9  -  m 


X 


(i) 


iii)  for  the  case  of  known  9  ,  we  may  use  Y  = 


Xi 

-r and 


n-m 

I  In  y 

r^h)  _  1=1 _ 


(i) 


r(h) 


n-m 


and 


rf(^) 

r(fi1) 


n 


i=n-m+l 


Jin  y 


(i) 


m 


iv)  for  the  case  of  known  k* 


9K  T  (n-m)q  mk*r) 
99  ~  q2  9  9 


5K 

9q 


n-m  n 

l  An  xm  +  k*  l 

i=l  ^  '  i=n-m+l 


Jin  x 


(i) 


(n-m)r'(n) 

r(  ri) 


=  0  . 


hence  obtain 


mk*r*  (k*r|) 

r(k*ri) 


; 
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and  thus 


T _ 

r){n-m(l-k*)  } 


and 


n-m  n 

l  x  +  k*  l 

i=l  ^  '  i=n-m+l 


In  x 


(i) 


(n-m)r(f) 

fmr 


mk*r,(k*f))  _  _ 
r(k*f|)  ” 


CHAPTER  IV 

The  Lognormal  Distribution 

We  now  examine  the  lognormal  distribution  as  a  life-testing 
distribution  that  is  a  competitor  to  the  gamma  distribution.  We  shall 
show  that  the  exchangeable  model  involving  the  lognormal  family  of 
distributions  indexed  by  the  shape  parameter  a  is  outlier-resistant 
completely  on  the  right.  We  shall  then  determine  that  as 
heterogeneity  increases  the  spurious  observation  tends  to  appear  as  an 
outlier  and  we  shall  determine  where  it  will  most  likely  appear.  We 
shall  also  consider  estimation  in  the  presence  of  this  outlier. 


-52- 


. 
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4.1  Characteristics  of  the  Lognormal  Distribution. 

If  a  random  variable  X  has  the  lognormal  distribution  with 
location  parameter  p  and  shape  parameter  a,  then  the  probability- 
density  function  of  X  is  given  by 

1  1  2 

- —  exp  {-  — j  (  Jin  x-  p)  }  ,  x  >  0  ,  a>0,-a°<p<c° 

ox/2u  2a 

0  ,  elsewhere 

We  shall  write  X  ~  A( p, a) .  The  following  table  summarizes  some  some 
of  its  properties: 


Table  4.1.1 


1  2 

M+  j  a 

mean 

e 

2  2  . 

variance 

0*0  ^  v  ^ 

e  (e  -l)e 

mode 

2 

p-a 

e 

median 

p 

e 

coefficient 
of  variation 

f  o2  ,1/2 
(e  -1) 

skewness 

2  2  i/o 

/  O  .  r)  n  ✓  O  -vl/Z 

(e  +2)(e  -1) 

kurtosis 

4 , 0  34.o  2  ,  a2 

co  +2  go  +3  go  -6,  oo  =  e 

The  lognormal  distribution  has  monotone  likelihood  ratio  (MLR)  in 

2 

T^(x)  =  Jinx,  a  known  and  in  T^Cx)  =  (Jin  x-p)  ,  p  known  but  not 
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n 


MLR  in  x  for  p  known  and  nonzero.  =  £  Jin  x_^  is  a  complete 


i=l 


n 


sufficient  statistic  for  p  (a  known)  and  (T~  =  £  (In  x,)' 

i=l  1 


T 

*  1 


n 


£  In  x  )  are  jointly  complete  sufficient  statistics  for  ( p, a  ). 


i=l 


n 


v  2  2 

T~  =  2  ( ^  x.-p)  is  sufficient  for  a  (p  known).  The  lognormal 

i=l  1 

distribution  can  assume  shapes  from  severely  right-skewed  to 


essentially  symmetrical. 
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b. 

•r*i 

u. 


Lognormal  p.d.f.  for  various  values  of  <r  (  =0) 
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If  X  ~  A(  |i,  a)  then  Y  =  Jin  X  ~  N(p,  a2). 

For  some  data,  the  lognormal  distribution  is  a  competitor  to  the 
gamma  distribution.  The  lognormal  p.d.f.  is  unimodal,  vanishes  at 

-a2  2 

x  =  0  and  its  mode  is  at  x  =  e  <1.  As  a  increases,  the  mode 
converges  to  zero.  From  the  above  graph  where  a  =  2,  the  lognormal 
p.d.f.  appears  monotonically  decreasing  for  almost  all  x  >  0  and  when 
0=5  the  mode  is  virtually  zero.  This  may  be  compared  to  the  gamma 
distribution  with  shape  parameter  q  <  1.  Here  the  density  is  infinite 
at  zero  and  monotonically  decreasing  thereafter.  The  lognormal 
distribution  has  a  non-monotonic  failure  rate. 

Neyman  and  Scott  (1971)  document  data  from  rain-making  experiments 
where  the  data  is  nonzero  rainfall  per  experimental  unit  (an 
experimental  day  or  storm).  Here  the  distributions  are  reverse 
J-shaped  with  long  "tails"  and  frequently  display  substantial 
"outliers".  They  suggest  that  an  outlier-prone  distribution  such  as 
the  gamma  or  lognormal  might  appropriately  model  this  data. 

Nelson  (1977)  points  out  the  usefulness  of  the  lognormal 
distribution  for  approximating  distributions  of  input  variables  such  as 
costs,  sales,  market  shares,  etc.  required  for  Monte  Carlo  simulations 
of  business  decisions.  This  distribution  has  been  used  by  Howard  and 
Dodson  (1961)  and  Peck  (1961)  to  study  semiconductor  devices  and  by 
Goldwaithe  (1961)  in  small-particle  statistical  economics  and  biology. 
Singpurwalla  and  Keubler  (1966)  used  it  to  study  the  lifetimes  of  high 
speed  steel  drills  since  the  failure  rate  first  increases  and  then 
decreases,  indicating  the  drill  could  resharpen  itself  and  prolong  its 
life.  It  has  wide  applicability  to  reliability,  especially 
maintainablity  and  fracture  problems.  If  X  <  X  <  ...  <  X  denote 


' 
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sizes  of  a  fatigue  crack  at  successive  stages  of  its  growth  and 

the  initial  size  of  the  crack,  assuming  a  "proportional  effect  model" 

for  growth  of  the  crack  (Kao  (1965)),  this  implies  crack  growth  at 

stage  i,  X^  -  X^_^,  is  randomly  proportional  to  the  size  of  the 

crack,  X.  ,  (i.e.  X.  -  X,  -  =  m.X.  ,,  i  =  1,2, ...,n  where  it. are 

l-l  l  i-1  i  l-l  ’  ’  ’  i 

independent),  and  the  item  fails  when  crack  size  reaches  X  .  It  can 

n 

be  shown  that  In  X^  is  asymptotically  normal  and  hence  X^  is 
lognormal.  This  model  is  also  true  for  the  distribution  of  oil  pool 
sizes  by  the  same  argument  on  how  the  pools  are  initially  formed. 
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4.2  Outlier-proneness  of  the  exchangeable  model  with  the  Lognormal 
distribution 

Neyman  and  Scott  (1971)  have  demonstrated  that  for  i.i.d.r.v.’s 
the  family  of  lognormal  distributions  indexed  by  a  is  outlier-prone 
completely  on  the  right.  If  we  now  consider  the  exchangeable  model 
with  at  most  one  outlier,  we  have  n-1  observations  from  A (  p,  a)  and 
one  observation  from  A(p^,o),  p^  _>  p  (Case  I)  or  n-1  observations 
from  A(  p,  a)  and  one  observation  from  A (p,  o^),  o^  >_  0  (Case  II). 

Case  I:  Scale  change 

We  first  consider  the  exchangeable  model  based  on  the  lognormal 
distribution  with  possible  change  in  p.  Then  f(x;j9)  is  A(  p,  a)  and 
f(x;£)  is  A( p^  ,  a)  where  p^  =  k*p,  k*  >_  1.  Kale  (1975b)  proved  that 
the  exchangeable  model  with  (at  most)  one  possible  outlier  observation 
for  scale  parameter  families  for  non-negative  random  variables 
involving  a  possible  change  in  scale  is  outlier-prone  completely  on  the 
right.  Thus  the  model  we  are  considering  would  be  outlier-prone. 

Case  II:  Shape  Change 

Consider  the  exchangeable  model  based  on  the  lognormal 

distribution  with  possible  change  in  shape  parameter  o.  Then  f(x;9) 

2  2 

is  A(0 ,  a)  and  f(x;_|)  is  A(0,o^)  where  =  k*a  ,  k*  _>  1 .  The 

likelihood  may  be  written  as 


, 
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L(x;cr,k*)  =  ^  I  n  f(x  ;9)f(x  ;_§), 
r=l  i^r 


x.  >  0  ,  k*  >  1, 
i  — 


a 


2 

1 


i  =  1,2 


.  ,n 


(k*  >_  1  since  we  are  considering  x(n)  as  a  possible  outlier).  Then 
the  joint  density  of  the  order  statistics  may  be  written  as 


£(X(1)(...,X(n))  nn!^if(Xr;a)  ^ 


n  f(x  jo.)  n 
I  —  —  n  f(x.,a) 


and 


s(x(l),x(n-l),x(n)) 


X(n-1)  X(n-2)  x(3 ) 

=  J  j  ...  j  f(x(1)>.-.-x(n))dx(2)-"dx(n-2) 

X(l)  X(l)  X(l) 


n 


*  (n-D!  I  hr(x(l)"x(n-U’x(n)) 


r=l 


where  Vx(l)  >x<n-l)’W  =  1 


f(xr;o1)  n 

-77 - c—  TI  f(x.;o)dx0  — dx  _ 

S  ,  f(xr;o)  1=1  1  2  n'2 

n-3 


n-3  is  the  reSion  X(l)  <  x(2)  <  •••  <  x(n-2 )  <  x(n-l)' 


and  S 

For  r  =  2 , . . . , n-2 


f(x(l):°)f(x(n-l);o)f(x(n);a) 
hr(x(l),x(n-l)’x(n))  (r-2) ! (n-r-2 ) ! 


(continued) 
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x( 


‘(n-1)  r-2 

I  [f (X(r ) 5 a)-F(X(1);  a) ]  [F(x(n_1 ^ ; a)-F(x(r) ;  a) ] 

X(l) 


n-r-2 


f(x(r);ai)dx(r)  , 


while 


hl(x(l),X(n-l)’X(n)) 


f(x(1);a1)f(x(n_1);a)f(x(n);a) 


fF(x(n-l);a)_F(x(l);a)] 


n-3 


(n-3 ) ! 


hn-l(x(l),X(n-l),X(n)) 


=  f(x(1);a)f(x(n_1);a1)  f(x(n);a) 


[F(x(n-l);a)"F(x(l);a)]n  3 

(^3)1 


y 


and 


hn(x(i),x(n_i),x(n)) 


=  f(x(1);a)f(x(n_1);a)f(x(n);a1) 


[F(x(n-1); a)_F(x(l); a) ^  3 

0^3)! 


Thus,  letting  t(y)  =  f(y;a)  and  s(y)  =  f(y;a^) 


(4.2.1)  .  . 


*(x(l)’x(n-l)>x(n)) 
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=  (n-l)!{t(X(1))t(x(n_1))t(x(n)) 


(n-1) '  iv"(r) 


(r) 


n-2  x(n-l)  [T(x(r))-T(x(1))]r  2  [T(x^_,  ,  )-T(x^)  ]n  r  2s(x,_,) 

X  I  I  - 

r=2  x(1) 


(r-2 ) ! (n-r-2 ) ! 


dx 


(r) 


+[s(x(1))t(x(n_1))t(x(n))+t(x(1))s(x(n_1))t(x(n)) 


+t(xa))t<X(n_1))s(X(n))] 


[T(x(n.i))-T(X(i))] 

(n-3)! 


n-3 


In  4.2.1,  we  may  replace  the  letter  of  integration  in  the  definite 
integrals  by  z  and  use  the  fact  that 


n-2 

(n-l )  !  I 
r=2 


r-2  n-r-2 

a  p 

(r-2) ! (n-r-2)! 


:  O'1)! 

(n-4 ) ! 


(a+P) 


n-4 


to  give 

g(X(l),X(n_i),X(n))  =  (n_3)  j  f(n“3)t(x(l))t(x(n-l))t(:x(n)  } 
[T(x(n-l))_T(x(n))1  [S(x(n-l))_S(x(l))] 

+  [^x(l)>t(x(n-l))t:(x(n))+t:(x(l))s(x(n-l))t(x(n)) 


(continued) 
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+  t(x(1))t(x(n_1))s(x(n))][T(x(n_1))-T(x(1))]n'3} 

=  tT(x(n-l))'T(x(l)))n"4l(n‘3)C<x(l))t(x(n-l))t(x(n))>< 

+  ts(x(l))t(x(n-l))t(x(n))+t(x(l))s(x(n-l))t(x(n)) 
+  t(x(1))C(x(n_1))s(x(n))][T(x(n_1))-T(x(1))]}. 

Let 


p  =  P{X(n)  >  (k+l)X(n_1)  -  kX  } 


®  x(n-l )  ® 

“  II  1 

0  0  (k+l)x(n_1)-kx(1) 


g(x(l)>x(n-l)>x(n))dx(n)dx(l)dx(n-l) 


and  let  x  =  x,..  .  and  y  =  x,  ,  N 
(1)  7  (n-1) 


Theorem  4.2.1:  If  p  =  P{X/  >  (k+l)X,  -  kX..,.}  then,  for  non- 

(n)  (n-1)  (1) 


negative  random  variables, 


n-3 


oo  y 

p  <_  l-k(n-l)(n-2)j  t(y)J  S(x)  (T(y)-T(x)  }“  ->t{(k+l)y-kx}dxdy 

0  0 


where  S  and  s  are,  respectively,  the  distribution  function  and 
probability  density  function  of  the  spurious  observation  and  T  and  t 
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are,  respectively,  the  distribution  function  and  probability  density 
function  of  the  non-spurious  observations. 


Proof : 

p  =  P{X(n)  >  (k+l)X(n_1^  -  kX(1)} 


co  x(n-l )  °° 

=  1  J  1  g(*m 

0  0  (k+l)X(n_1)-kx(1) 


,X(n-l ) 


,x(n))dx(n)dx(l)dx(n-l) 


(n-1 )  ! 
(n-3 )  ! 


00  X(n-1 ) 


I  J 
0  0 


{T(x(n.1))-T(Ka))} 


n-4 


{(n-3)t(x(,1))t(x(,n_1))[l-T{(k+l)x(,n_1)-kx(,1)}]  [S(x(.n_1))-S(x(1))] 
+  (s(x(1))t(x(n_1))[l-T{(k+l)x(n_1)-kx(1)}] 


+  t(x(1))s(x(n_1))[l-T{(k+l)x(n_1)-kx(1)}] 

+  t(x(1))t(x(n_1))[l-S{(k+l)x(n_1)-kx(1)}] ) 

*  lT(x(n-l))'T(x(l))lldXa)dX(n-l) 


(n-1)! 

(n-3)! 


(Il+I2+I3«4  -  (I5+I6+I7+I8)}. 
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Let  x  =  x(1)  ,  y  =  x(n_1} 


I  =  j  j  (n-3)t(x)t(y) (T(y)-T(x) }n  4{S(y)-S(x) }dxdy 
0  0 


(n-3) J  S(y)t(y) ]  t(x) (T(y)-T(x) }n  4dxdy 
0  0 


-(n-3)J  t (x)S (x) J  t(y){T(y)-T(x) }n  4dydx 
0  x 


00  n-3  CO 

(n-3)J  S(y)t(y)-T^^^ —  dy  -  (n-3)j  S  (x)  t  (x)— 


=  j  t(x)S(x)[{T(x)}n  3-{l-T(x)}n  3]dx 
0 


I  =  J  J  (T(y )-T(x) }n  3s(x)t(y)dxdy 
0  0 


CO  00 


n-3 . 


=  J  s (x) J  (T(y)-T(x) }  t(y)dydx 
0  x 


=  j  (1-T(x)}n  2dx. 


0 


(n-2) 


n.— 3 

Integrating  by  parts  where  u  =  (l-T(x)}  dv  =  s(x)dx 


du  =  -(n-2)t(x) (l-T(x) }n  3dx  v 


3 

—  dx 


=  S  (x) 
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then 

00 

I O  =  j  s(x)t(x) (l-T(x) }n  3dx. 
0 


go  y 

1 o  =  /  J  (T(y)-T(x)}n  3t(x)s(y)dxdy 
0  0 


co  y 

=  J  s(y)  J  {T(y)-T(x)}n  3t(x)dxdy 
0  0 


s(x) 


(T(x) }n  2 

(n-2) 


dx 


Integrating  by  parts  where  u 


{T(x)}n  2  dv  =  s(x)dx 


du  =  (n-2){T(x)}n  3t(x)dx 


I,  =  — ^  ~  ]  S(x)t  (x)  (T(x)  }n  3dx 
3  n"2  0 


00  y  _o 

I  =  j  J  t(x)t(y){T(y)-T(x)}n  Jdx 
0  0 


CO  oo 

=  i  t (x) J  t(y) (T(y)-T(x) }n  3dydx 
0  x 


00 


-  J  t(x) 

o 


(l-T(x) }"  2 

(n-2) 


dx 


=  S(x) 
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(n-1 ) (n-2 )  * 


Therefore  I,  +  I_  +  I_  +  I ,  =  ^  +  1 


1  2  3  4  (n-2)  (n-1 ) (n-2 ) 


n 


(n-1 ) (n-2 )  * 


n-4 


I  =  J  J  (n-3)t(x)t(y) {T(y)-T(x) }  T{(k+l)y-k(x) }{S(y)-S(x) }dxdy 


0  0 


>  J  J  (n-3)t(x)t(y){T(y)-T(x)}n  4T(y) (S(y)-S (x) }dxdy 
0  0 


J  (n-3)t(y)T(y)S(y) J  t (x) (T(y )-T(x) }n‘4dxdy 
0  0 


-  ]  (n-3)t(x)S(x) J  t(y) (T(y)-T(x) }n  3dydx 

0  x 


-  J  (n-3)t (x)T(x)S(x) J  t(y){T(y)-T(x)}n  4dydx 
0  x 


00  CO 

=  J  t(x)T(x)S(x)  (T(x)  }n  3dx  -  J  -|'^3 ^  t(x)S(x)  {l-T(x)  }n  2dx 


0 


0 


-  ]  t (x)T (x)S (x) (1-T (x) }n  3dx. 
0 


67 


Therefore 


IS  -  J  t (x)S (x) (T(x) }n  2dx 
0 


I  t(x)S(x)  (l-T(x)  }n  2dx 
(n"2)  0 


-  I  t(x)T(x)S(x) (l-T(x) }n  3dx 
0 


I  >  J  ]  t(x)s(y)T(y){T(y)-T(x)}n  3dxdy 

0  0 


I  s(y)T(y)J  t(x){T(y)-T(x)}n  3dxdy 
0  0 


j  s  (x)T(x)  -  dx 

o  (n'2) 


1:0  fr-,  .,n-l 

=  J  s (x)  -  dx  . 


0 


(n-2) 


n-1 


Integrating  by  parts  where  u  =  (T(x)}  dv  =  s(x)dx 


du  =  (n-l)t(x) (T(x) }n  2dx  v  =  S(x) 


I 


(n-1) 

(n-2) 


J  t (x)S (x) (T(x) }n  2  dx  . 
0 


00  y  _o 

I  >_  j  j  t(x)t(y)S(y){T(y)-T(x)}n  dxdy 
0  0 


t(y)S(y) 

(n-2) 


(T(y) }"  2dy 


00 


68 


=  -(-n_2 )  j  t(x)S(x)  {T(x)  }n  dx. 


I,  =  J  t(y)J  s(x) (T(y)-T(x) }n  3T { (k+1 )y-kx}dxdy  . 
0  0 


Integrating  by  parts  where 


u  =  {T(y)-T(x)}n  3T{(k+l)y-kx}  dv  =  s(x)dx 


du  =  {T(y)-T(x)}n  3t{(k+l)y-kx}{-kdx)  v  =  S(x) 


+  (n-3) {T(y)-T(x)}n  4T{(k+l)y-kx}{-t(x) }dx 


ao  y 

I  =  k  J  t(y )J  S(x) (T(y)-T(x) }n  3t {(k+l)y-kx}dxdy 
0  0 


+  (n-3) J  t(y) J  S(x) (T(y)-T(x) }n  4T{(k+l)y-kx}t(x)dxdy 
0  0 


®  y  _o 

>  k  j  t(y) j  S(x) (T(y)-T(x) }n  t {(k+l)y-kx}dxdy 
0  0 


co  y 

+  (n-3) j  t(y)J  S(x) (T(y)-T(x) }n  4T(y ) t (x)dxdy  . 
0  0 
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But 


°o  y 

(n-3)j  t(y)J  S(x){T(y)-T(x)}n'4T(y)t(x)dxdy 
0  0 


00  oo 

=  (n-3) ]  t(x)S(x)J  t(y) (T(y)-T(x) }n  3dydx 
0  x 


00  oo 

+  (n-3)]  t (x)T(x)S (x) j  t(y){T(y)-T(x)}n_4dydx 
0  x 


(n-3) 

(n-2) 


t(x)S(x) (l-T(x) }n  2dx 


+  J  t(x)T(x)S(x) (l-T(x) }n  3dx 
0 


00  oo 

1  t(x)s(x)  (l-T(x)  }n  2dx  -  J  t(x)S(x)  (l-T(x)  }n  2dx 
K  0  0 


+  J  t(x)S(x) (l-T(x) }n  3dx. 
0 


OO 

r\ 

Therefore  I-  +  I.  +  I,  +  IQ  >  j  t (x)S (x) (T(x) }n  dx 
5  6  7  8  -  Q 


(n-3) 

(n-2) 


OO  oo 

j  t (x)S (x) (l-T(x) }n  2dx  +  J  t (x)S (x) (l-T(x) }n  2dx 
0  0 


(continued) 
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-  J  t (x)S (x) (l-T(x) }n  3dx 
0 


CO  y 

+  kj  t(y)J  S(x){T(y)-T(x)}n_3t{(k+l)y-kx}dxdy 
0  0 


+  77TTT  J  t(x)S(x){l-T(x)}n  2dx  -  j  t (x)S (x)  (l-T(x)  }n  2dx 
^  ;  0  0 


+  /  t(x)S(x){l-T(x)}n  3dx  + 

n  \n  ^  / 


J  t(x)S(x){T(x)}n_2dx  +  -7-^-r  j  t(x)S(x){T(x)}n"2dx 

^  ;  0  v ;  0 


Therefore 


P 


<  (n-1 ) (n-2) 


r _ 2 _ 

L(n-l)(n-2) 


1 

(n-2) 


-  k J  t(y)J  S(x) (T(y)-T(x) }n  3t{(k+l)y-kx}dxdy ] 
0  0 


=  1  -  k(n-l ) (n-2 ) J  t(y)J  S (x) [T(y )-T(x) ] n  3t {(k+1 )y-kx}dxdy  . 

0  0 


Theorem  4.2.2:  The  exchangeable  model  with  (at  most)  one  spurious 
observation  based  on  the  lognormal  family  of  distributions  indexed  by 


shape  parameter  a  is  outlier-resistant  completely  on  the  right. 
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Proof:  From  Theorem  4.2.1  we  know  that  if 


p  -  P{X(n)  >  (k+l)X^_1)  -  kX(1)} 


then,  for  non-negative  random  variables, 


where  x 


distribution  function  and  p.d.f.  of  a  spurious  observation,  and  T  and 
t  are,  respectively,  the  distribution  function  and  p.d.f.  of  a  non- 
spurious  observation.  If,  for  all  k  >  0,  n  >  2  sup  p  <  1  then  the 
family  is  outlier-resistant  completely  on  the  right. 

It  is  then  sufficient  to  show  that 


y 

(4.2.2)  H  =  k(n-l)(n-2)J  t(y )j  S(x) (T(y)-T(x) }n  t { (k+1 )y-kx}dxdy 

0  0 


y 


>  0. 


Let 


in  x 


In  y 


v 


a 


a. 


1 


If  <})  and  <5  respectively  denote  the  standardized  normal  p.d.f.  and 
d.f.  then  expression  (4.2.2)  becomes 
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cn  S2  aS2  0S1 

k(n-l)(n-2)  j  <Ks2)J  $(vS;l)  {<S>( s 2 )-^>( s x )  }n~3({){-J^[-(k+1.^ - — - ]-} 

—  co  —co 


as. 


e  ds^ds^ 

as_  as. 
(k+l)e  -ke 


as . 


5  2 

>_  k(n-l ) (n-2 )  J  <t>(s?)  j 
3  0 


<3?(vs^)e 


as?  as 
(k+l)e  -ke 


(®(s2)-$(s1) } 


n-3 


as  9  as  - 

x  — L}dSlds2 


In  this  region,  s2  >_  3,  >_  0,  s2  >.  1  +  s^.  Therefore  3?(vs^)  >_  1/2 


as^ 

for  all  v(i.e.  for  all  a^),  e  1 , 


-as, 


$(s2)  -  $(s1)  >_  $(3)  -  $(2)  and  e  2  >  e  5  a. 


Then 


H  >  k-(~n2(k+iy2''  f^(3)  "  $(2)^ 


n-3  1  “5  a  „  ,  .  s  r  £n(k+l)>, 

j  e  2 J  <Ks2)<Hs2+  - - - Jds2 


Let  u  =  ^ "  .  Then  u  >  0  since  k  >  0  and 


5  1  5  1 

J  <|>(s9)<|>(s9+u)ds?  =  -  1  -  e 

3  /2n  3  >/ 2 tc 


s0+u/2  2 

-1/2  [-i— ]  -u2/4 

1//2 


ds. 
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-u2/4  /2(5+u/2) 

- - - -  j  _  <j)( z )dz 

/2n  /2  /2(3+u/2) 


-u2/4 

e-  -  [a{/T(5+u/2)>  -  ${/2  (3+u/2  )  }  ] 

2/ti: 


>  0. 

Then  for  all  k>0,  n>2,  v=a/a^>0,H>0  and  p  =  1-H  <  1  and 
therefore  sup  p  <  1.  Thus  the  exchangeable  model  with  (at  most)  one 
spurious  observation  involving  the  lognormal  family  of  distributions 
indexed  by  shape  parameter  a  is  outlier-resistant  completely  on  the 


right . 
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4.3  Detection  of  Outliers 


Case  I:  Scale  change. 


If  n-1  observations  are  from  A(  p ,  a)  and  one  observation  is 
from  A(p^,a),  p^  =  k*p,  k*  >_  1  then 


Jin  x-p-  2 

L) 


=  dG(x) 
nx;  dF(x) 


ox/2x 


J!n  x-li  2 
-1/2  - ) 


ox-/2iz 


1  2  2 
=  exp  [ - j  {(Jin  x-p  )  -(Jin  x-p)  }] 

2  a 


2 


=  [exp  { - ^ — }] 


x 


2  a 


P1-p 


-1 


M-i  “M- 

and  Y’(x)  =  ( — r— ) 


^  ^  a2  r  12  2 

2  ;x  exPl - f  ^ 

a  2  a 


Thus  V  (x) 


>  0  if  p^  >  p 


<  0  if  p^  <  p 


Consequently  ^(n)  has  maximum  probability  of  being  the  spurious 
observation  if  p^  >  p;  has  maximum  probability  of  being  the 
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spurious  observation  if 


^  <  p. 


Case  II:  Shape  change 


If  n-1  observations  are  from  A(  p,  cr)  and  one  observation  is 


2  2 

from  A(  p,  o^),  =  k*cr  ,  k*  >_  1  then 


'f(x)  =  frfy =  k expt"  r (ta  x-^)2(jT  -  ^ 

1  o^o 


is  not  monotone  in  x. 


Consider  u(r;n,o,k*)  =  P(X^^  is  the  spurious  observation)  (i.e. 
X^^  is  the  order  statistic  whose  distribution  has  parameter  o^). 

Then 


u(r;n,o,k*)  =  {F(y;o)}r  1{l-F(y;o)}n  rf(y;a1)dy 

0 


1  r  Jin  y-  p 

exp  -  T{ - 

n-r  2  ^  a. 


a^y/2-n; 


-dy 


where  $  is  the  standard  normal  distribution  function 
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01 


n-r 


1  fin  y-|i^ 
exp  '  2-1^-^ 

a^y/2-n; 


■dy 


and  setting  x  =  ■^n-  -  -  we  have 


,  00  _  ,  _  _  exp  {-  y~ } 

u(r;n,  cr,k*)  =  (n_-  )J  {$(x/k*)}r  (l-<J>(x/k*)  }n  r  - — -  dx 

r  -00  /2x 


Now 


j  0  if  x  <  0 

0  $(x/k*)  <_  1  and  lim  <$(x/k*)  =  J 

k*-**>  (  1  if  x  >  0 


Therefore 


2/0 

.  00  _  ,  _  -x  /2 

lim  u(r ; n,  o,k*)  =  (n_J)  lim{J  {$(x/k*)  }r  {l-$(x/k*)}n  r 
k*->-oo  r  k*+°°  0  •  /2ti: 


0  _  _  _  -x2/2 

+  J  {$(x/k*)}r  i{l-$(x/k*)  }n  r  - — — —  dx} 

/2n 


dx 


By  the  dominated  convergence  theorem 
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lim  u(n;n, o,k*) 
k*-*00 


00  -x2/2 

]  ^—zzr  dx  =  1/2 

0  /2ti 


0  -x2/2 

lim  u(l;n,o,k*)  =  j  - — -  dx  =  1/2 

k*>°°  -00  /2ti 


and  for  any  other  r  =  2,...,n-l 


{$(x/k*)}r  1  {l-<£(x/k*)  }n~r  +  0  as  k*  ->  <*>  , 


so  that,  by  the  dominated  convergence  theorem, 

<*»  -x2/2 

lim  u(r;n,o,k*)  =  0  — — -  dx  =  0,  r  =  2,...,n-l. 

k*-XD  -CO  /2  tc 

It  thus  appears  that  the  spurious  observation  tends  to  occur  at 
either  or  X^^  with  equal  probability  (approaching  1/2)  as  the 

shape  parameter  of  the  spurious  observation  increases. 
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4.4  Estimation  for  Lognormal  parameters 
4.4.1  Standard  Estimators 

If  X  ~  A(p,a),  M.L.E.'s  of  p  and  a  may  be  obtained  from 

? 

M.L.E.  of  W  =  In  X  since  W  ~  N(p,o  ).  Thus  the  M.L.E.’s  of  p  and 
a  are 


n 


n 


H  = 


l  wd)  £ 

i=l  v;  i=l 


Jin  X 


(i) 


n 


n 


n 


n 


-2  i=l 
a  =  - 


2  J  X(i)  2 

I  (to  X  )  -  n(LL -  ) 


n 


Since  p  and  o  are  based  on  complete  sufficient  statistics,  we  can 

2  ~2 

obtain  best  linear  estimators  p  and  s  =  --gT-  and  n , -  s  has 


n-1 


n+1 


2  2 

minimum  MSE  among  estimators  of  a  of  the  form  ks  .  The  BLIE 

2  ~2  n  -2 

estimators  of  p  and  a  are  p  and  a  =  -^j^-  a  .  Asymptotically, 

all  three  estimators  are  equivalent.  However,  using  the  property  of 
MLE ’ s  that  if  y  is  MLE  of  y,  g(y)  is  MLE  of  g(y),  we  may 

obtain  M.L.E.’s  of  the  mean  and  variance  of  the  lognormal  to  be 

q,  1  p.2 

P+  tt  <3  0  p2 

.  2  ,  2p+a  ,  a  -1N 

respectively,  e  and  e  (e  ). 


To  obtain  confidence  bounds  for  parameters  of  the  two-parameter 


lognormal  distributions,  we  use 
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p.  =  W  ~  N(  p,,  o'/n)  . 

2 

Thus,  for  a  known,  a  (1-a)  100%  confidence  interval  (C.I.)  for  p. 

is  given  by 


w  +  z ^^2  ~  |-i  _<  w  +  z 


v^n 


l-a/2  * 

/n 


n 


„  2  no 

For  a  unknown,  — - 


V  2  -2 

L  w  -  nw 

i=l  2 

- 2 - Vl  df 


and  ( 


o/  /n  /  f  o  (n-1 ) 


t  .  df 
n-1 


and  a  two-sided  (1-a)  100%  C.I.  for  p.  is 


w  +  t  -  <  p.  <  w  +  t.  /o 

a/2,n-ldf  > — = —  l-a/2, n- 

/n-1 


ldf 


/n-1 


4.4.2  Estimators  suggested  for  use 


Case  I:  Scale  change 


Consider  the  estimator  |i  = 


n 


l  h i  X. 


i=l 


n 


Under  the  homogeneous 
2 


model,  E((i)  =  p,  and  Var(|l)  =  .  Thus  its  MSE(p,)  =  . 


Under 
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the  exchangeable  model  with  (at  most)  one  possible  outlier, 


Ehet 


1_ 

n 


n-1 

n 


=  p{l  + 


k*-l 

n 


) 


and 


Varhet(^  “  T-  X  ^et  V 


n  i=l 


n 


1  v  21  2n 

=  ~T  L  hr  *  +  TT  <*  ) 


2  •  i 

n  i=l 


n 


n 


We  then  obtain 


MSEhet(^  =  Varhet(il)  +  {Blashet(W}' 


_  i_  +  m(k*-l)i2 
n  *•  n  ' 


n 

y  jin  x. 

i=l  1 

and,  as  k*  ->  00 ,  MSE^e{_(|i)  ->  00  .  Thus  p  =  - — -  appears  to  be  a 

poor  estimator  in  this  instance. 
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Assuming  the  exchangeable  model  with  m  outliers  (generally 

m  <_  10%  of  the  sample  size),  we  have  n-m  observations  with  p.d.f. 
112 

f(x;p,  a)  =  — — —  exp  -  — ^  (to  x-p)  and  m  observations  with  p.d.f. 
x  a/2  it  2  a 

f(x;p^,a),  p^  y_  p,  The  joint  likelihood  may  be  written  as 


L(x;  p,  px 


a,I)= 


l  (to 


2a  x^/I 


n  ax./2n: 
x^  H 


I  r-  2 

- 7  l  (to  x-p  ) 

2a  x  el 

II  ax./2-n; 
x^  el 


We  wish  to  obtain  max  L(x;  p,  p.  ,  a, I) 

I  ef  ~  1 

px>p 

a>0 

all  possible  subsets  of  m  outliers 
(1974b)  has  shown  for  Y(x)  = 

dF(x) 

A 

I  =  (x ,  , - N , . . . , x ,  v }  has  maximum 

(n-m-fl)  (n) 

spurious  observations.  For  p^  >  p, 
hence 


where  J  is  the  collection  of 

in  a  sample  of  size  n.  Kale 
monotone  increasing  in  x, 
probability  of  being  the  set  of 
Y(x)  is  monotone  increasing  and 


max  L (x;  p,  p.  ,  a, I 
lef  ~  1 

PX>P 

a>0 


) 


max  L(x;  p,  p^,  a,I)  . 

P;L>p 

a>0 


Now  K(x;  p,  p^,  a, I)  =  An  L(x;  p,  p^ ,  a, I ) 


=  C 


-  n.  m  o  ri  ^ 

17  {Ji (to  x^)^r+iJm+i(te  ] 
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-  n/2 


In  2na 


I 


i=l 


In  x . 
1 


n-m 


,  5K  1  r  .  .  N 

Then  af  "  J  (  (i)'1*5 


n 


fr  =  i2.  ‘t  <tox(ir,1i) 

1  a  i=n-m+l 


^2  =  \  {  I  (An  xm“^>Z  +  l  \ 

5  a  a  i=l  i=n-m+l  u;  a 


n-m 


n 


2 1  n 


n-m  In  x,,> 

and  hence  p,  =  Y  - - — 

Hiet  . ,  n-m 
i=l 


n 


l  An 


x 


Hi 


i=n-m+l 


(i) 


bet 


m 


n-m  2  n  2 

(to  x(ifiVet)  +.  I  (*>  x  ) 

£  =  ]_  v  y  -f  =n—  m-4-1  v  '  Kq+- 


i=n-m+l 


bet 


bet  n 

Thus  it  appears  that  maximum  likelihood  estimation  suggests 


trimmed  means  as  estimators  of  p^  and  p  (p^>p)  and  a  pooled 

2 

estimator,  for  a  ,  that  utilizes  these  trimmed  means. 
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Case  II.  Shape  change 


n 


I  x 


Under  the  homogeneous  model,  E(p)  =  e( 


i=l 


n 


-)  =  l_i  and 


Var(p)  =  — .  Under  the  exchangeable  model  with  at  most  one  possible 
outlier, 


Ehet(«  -  E( 


n 

l  *n  x 

i=l 

n 


)  -  hi  ^  - 


n  i=l 


and 


n 


£  In  x , 


Varhet^>  "  Var^ - +  “  h  X  Varhet(to  xi> 


n  i=l 


n 


1  v  rn_1  2  ,  1  2^ 

=  -S-  l  [ -  a  +  —  a.  j 

2  .  u.  ^  n  n  1 J 

n  1=1 


a  (n-l+k*) 
2 


n 


Thus  MSEhet((l)  =  Varhet(p)  +  (Bias^p)}' 


a  (n-l+k*) 
2 


n 


which  tends  to  infinity  as  k*  -»  °°.  Thus  p  appears  to  be  a  poor 


estimator  of  p. 


. 
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Since  in  this  instance  the  spurious  observation  tends  to  appear  at 


N  or  at  X,  an  estimator  of  li  based  on  the  A-rule,  W-rule. 
(1)  (n)’ 


or 


n 


y  jbi  x. 

-  i=l 

S-rule  (see  Appendix  I)  would  appear  preferable  to  p  =  - - - 

2  2 
Similar  estimators  for  o  appear  preferable  to  s 

xn  X 


CHAPTER  V 

The  Weibull  Distribution 

We  shall  show  that  for  i.i.d.r.v.'s  the  family  of  Weibull 
distributions,  indexed  by  the  shape  parameter  p,  is  outlier-prone 
completely  on  the  right.  On  the  other  hand,  we  shall  show  that  the 
exchangeable  model  based  on  the  Weibull  family  of  distributions  indexed 
by  the  shape  parameter  is  outlier-resistant  completely  on  the  right. 

For  this  exchangeable  model,  we  shall  determine  which  observation  is 
most  likely  to  be  the  spurious  one  and  find  asymptotic  limits  for  these 
probabilities  as  heterogeneity  increases.  Also  we  shall  examine  the 
problem  of  estimation  in  the  presence  of  outliers. 


-85- 


5.1  Characteristics  of  the  Weibull  Distribution 


The  random  variable  X  is  said  to  have  the  three-parameter 
Weibull  distribution  if  the  probability  density  function  of  X  is 

(  (~g~)T1  1  exp  >  x>x,  p>0,  9>0,  x> 


f  (x;  9,  p,  x)  =  \ 

0 

,  otherwise 

giving  a  d.f.  of  F(x;9,p,x)  =  1-  exp  {-  (— g— •)T1}  ,  and  hazard  rate 


(intensity)  of  h(x;9,p,x)  = 

p(x-x)T1"1 

=  — - - -  .  In  this  case  we  write 

9  71 

X  ~  WEI (  9,p,  x). 

Figure  5.1.1  illustrates  the  shape  of  the  probability  density 
function  f(x;l,p,0)  for  various  values  of  p. 
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bo 
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Weibull  p.d.f.  for  various  values  of 
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The  Weibull  distribution  has  many  applications  in  life-testing  or 
in  problems  where  a  skewed  distribution  is  required.  It  may  be 
considered  to  be  a  generalization  of  the  gamma  distribution  to  allow 
the  hazard  rate  to  depend  on  a  power  of  X.  In  contrast  to  the 
exponential  distribution  with  constant  hazard  (failure)  rate 
h(x; 9)  =  (i.e.  p  =  1),  the  Weibull  distribution  allows  decreasing 

hazard  rates  (i.e.  p  <  1  implies  work-hardened  materials)  or 
increasing  hazard  rates  (i.e.  p  >  1  implies  wearout).  Many  items, 
especially  nonelectronic  parts,  exhibit  increasing  failure  rates. 
Leiblein  and  Zelen  (1956)  used  it  to  model  ball-bearing  failures;  Kao 
(1959)  used  it  to  model  vacuum-tube  failure.  Weibull  (1951)  derived  it 
in  the  analysis  of  breaking  strengths.  Whatever  the  underlying 
distribution  for  positive  random  variables,  the  Weibull  distribution  is 
one  of  only  two  possible  limit  laws  for  min(X^ , . . . ,X^)  as  n  -*  00 .  As 
such  it  is  used  to  model  metal  fatigue  breaking  strength  where  strength 
is  that  at  the  weakest  flaw  (e.g.  breaking  strength  of  chain,  ceramics, 
lumber,  concrete,  aircraft  parts,  etc.).  In  many  applications  the 
location  parameter  t  is  known  and  without  loss  of  generality  we  may 
take  t  =  0. 

There  exists  a  relationship  between  WEI(9,p,0)  and  the  Extreme- 
Value  distribution  EV  (£,b)  given  by  F(y)  =  1-exp  {-  } , 


<  y  <  °°, 


—  oo 


<  l  <  «,  b  >  0.  If  X  ~  WEI( 9, p, 0)  then 


' 
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Y  =  An  X  ~  EV^.(£  =  In  9,  b  =  1/ p) .  Thus  the  Weibull  distribution 
competing  with  the  lognormal  distribution  is  similar  to  the  EV^. 
distribution  competing  with  the  normal  distribution.  For 

V 

Y  ~  EV_(5,b),  E(Y)  =  l  -  by  and  Var(Y)  =  where  y  =  .5772  . 

1  o 

(Euler’s  constant). 

For  known  shape  parameter  r),  we  may  use  the  transformation 

Y  =  X1"1  to  obtain  Y  ~  EXPCQ^)  =  GAM(  9^,  1 ,0) .  Outliers  here  may  be 
handled  using  methods  for  the  single  parameter  exponential  family  (see 
Chapter  II). 
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5.2.  Outlier-proneness  of  the  Weibull  and  related  models 

First,  let  us  restrict  ourselves  to  the  standard  subfamily  indexed 
only  by  shape  parameter  p.  The  probability  density  function  is 


f (x; p) 


p-1  -x 
px  e 


0 


x  >  0  ,  p  >  0 
otherwise 


The  corresponding  distribution  function  (d.f.)  will  be  denoted  as 

F(x;p)  and  the  family  of  such  distribution  functions  as  f  .  In  this 

model  we  consider  i.i.d.  random  variables  X,,...,X  with  order 

1  n 

statistics  X,, N  <  X/0.  <  ...  <  X,  .  . 

(1)  (2)  (n) 


Theorem  5.2.1:  For  i.i.d.  random  variables  the  family  of  Weibull 
distributions  is  outlier-prone  on  the  right. 


Proof;  In  order  to  prove  Theorem  5.2.1  we  notice  that 


P(X(n)  >  X(n-1)  +  k(X(n-l)  X(l))  n  X(l)  >  0} 


p{x(1)  >  o}.p{x(n)  >  x(n.:)  +  k(x(n.1)-x(1))|x(1)  >  0) 


P{X(n)  >  X(n-1)  +  k(X(n-l)-X(l))}P{X(l)  >  °lX(n)  >  X(n-1) 


+  k(X(n-l)-X(l))> 


Now  p(xn)  >  *  1  and 


P{X(1)  >  °lX(n)  >  X(n-1)  +  k(X(n-D-Xa))}  ‘  1 


' 
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so  that,  following  Neyman  and  Scott  (1971), 


P(k,nh)  '  P{X(„)  >  X(n-1) 


k(X(n-l)-X(l))}  "  P(X(n)  >  X(n-1 ) 


k(X(n-l)-X(l)) 


X 


(1) 


>  0} 


>  P (X(n)  >  (k+1  )X(n_1 } }  =  Q(k,n|ri)  . 

Thus  it  is  sufficient  to  prove  the  stronger  assertion  that,  as  r)  ->  0, 
Q(k, n |  ri)  ->  1 . 

00 

Q(k,n|p)  =  n  ]  Fn  1  ;  p)f  (x;  ri)dx 


Setting  u  =  1 


-f— l"1 

kk+l  ' 

e  we  obtain 


Q(k,n  p) 


1 

n(k+l ) 


0 


n-1,-  >> 

u  (1-u) 


^^du 


nr(n)r{(k+n)T1}(k+l)T1 

T {n+(k+l ) 


It  then  follows  that,  as  p  ->  0,  Q(k,n|  p)  ->  1  and  thus  P(k,n|  p)  -*  1 
and  II^(k,njJ")  =  1.  Since  this  result  holds  for  all  k  >  0,  n  >  2,  the 


family  of  Weibull  distributions  is  outlier-prone  on  the  right. 
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Consider  now  the  exchangeable  model  with  (at  most)  one  outlier 
observation.  Kale  (1975b)  has  shown  that  the  exchangeable  model 
based  on  scale  parameter  families  for  non-negative  random  variables 
with  a  possible  change  in  scale  is  outlier-prone  completely  on  the 
right.  This  would  apply  to  the  exchangeable  model  based  on  the  Weibull 
distribution  with  possible  change  in  the  scale  parameter  9.  Let  us 
consider  the  other  possibility  -  the  exchangeable  model  involving  the 
Weibull  distribution  with  possible  change  in  the  shape  parameter  T). 

We  shall  show  that  this  model  is  outlier-resistant  completely  on  the 
right. 

Consider  the  situation  where  n-1  observations  are  from 
WEI(l,p,0)  and  one  observation  is  from  WEI(1 ,k*p,0) ,  0  <  k*  1  and, 
a  priori,  each  observation  is  equally  likely  to  be  the  spurious  one. 
Then  f(x;_9)  is  WEI(l,p,0)  and  f(x;£)  is  WEI (1  ,k*p, 0)  and  the 
likelihood  may  be  written  as 

1  n 

L (_x;  p,k*)  =  -  l  n  f(x  ;9)f(x  ;£),  x  >  0,  0  <  k*  <  1 
n  r=l  i*r 

(k*  <  1  since  we  are  considering  as  a  possible  outlier).  Then 

the  joint  density  of  the  order  statistics  may  be  written  as 

f 


(x 


(1)’ 


n 


n! 


r=l 


f  (xr;  pk*) 

f(xr;  p) 


n 

n 

i=l 


f  (x1;  p) 


and 


- 
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8(*(l)’*(n-iyx(iO) 


X(n-1) 

X 

(n-2) 

X(3) 

=  J 

J 

•  •  • 

j  f(x(1)>X{2), 

X(l) 

X 

(1) 

(1) 

n 

-  (n-1)! 

I 

r=l 

hr(x 

(l)’x(n-l)’x(n)) 

(n)'(2)*  *  (n-  2) 


where 


hr(x(l)’x(n-l)>x(n))  =  •> 


f (xr;  \ )  n 

- r —  n  f  (x. ;  ri)dx  .  .dx, 

f(x  ;  ri)  i  (2)  (n-2 ) 

S  ^  r  l-l 

n-3 


and  S  0  is  the  region  x..  .  <  x,ON  <  ...  <  x,  <  x,  ..  . 

n-3  s  (1)  (2)  (n-2)  (n-1) 

For  r  =  2 , . . . ,n-2 


,  ,  N  f(x(l);,1>f(x(n-l);T’:,f<x(n);,l) 

hr(X(i),X(n_i),X(n))  (r-2) ! (n-r-2) ! 


x . 


J  (F(x^x;  ti)-F(x/1  ri)  }r  2  (F(x/__1  ti)-F(x,^  ;  n)  }n  r  2 


x 


(r)’  ^v~(l) 


(n-1)’  iV^(r) 


(1) 


f(x(r);n1)dx(r) 


while 
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hl(x(l)’x(n-l)’x(n)) 


=  f(x(1);Tll)f(X(n_1) 


;  Tl)  f  (X(ny  Tl) 


[F(x 


(n-1) 


;  ti)-F(x^;  tO] 


n-3 


(n-3 ) ! 


and 

hn-l (x(l ) ,x(n-l ) ,X(n) ^ 


=  £(x(1);9)f(x(n.1);11l)f(x(n);T1) 


[F(X(n_i  );'n)~F(X(i);r|)] 

(n-3Ti 


n-3 


and 


hn(x(l)’x(n-l),x(n)) 


=  f(x(1);r1)f(x(n_1);T,)£(x(n);^1) 


[F(x(n_i);9)-F(x(1) 

(n-3)! 


Letting  t(y)  =  f(y;ri)  and  s(y)  =  fCyjr^)  we  obtain  again  expression 
(4.2.2)  and,  setting 


p  ”  P^X(r)  >  (k+1)X(n-l)  kX(l)^ 


(n— 1)  00 

’ll  I 

0  0  (k+1 )x 


(n-1)  kX(l) 


g(X(i)>X(n-1),x(n))dx(n)dx(l)dx(n-l) 
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and  x 


n-3 


t {(k+1 )y-kx}dxdy 


where  S  and  s  are,  respectively,  the  distribution  function  and 
probability  density  function  of  the  spurious  observation  and  T  and  t 
are,  respectively,  the  distribution  function  and  probability  density 
function  of  the  non-spurious  observations. 

Theorem  5.2.2:  The  exchangeable  model  with  (at  most)  one  spurious 
observation  based  on  the  Weibull  family  of  distributions  indexed  by 
shape  parameter  p  is  outlier-resistant  completely  on  the  right. 

Proof:  From  Theorem  4.2.1  we  know  that  if 


p  =  P{X(r)  >  (k+l)X^n_1}  -kX(1)} 


then,  for  non-negative  random  variables 


where  x  =  x^^,  y 


distribution  function  and  p.d.f.  of  a  spurious  observation  and  T 
and  t  are,  respectively,  the  distribution  function  and  p.d.f.  of  a 


non-spurious  observation. 


' 
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If,  for  all  k  >  0,  n  >  2  sup  p  <  1  then  the  family  is  outlier- 
resistant  completely  on  the  right.  It  is  then  sufficient  to  show  that 


go  y 

(5.2.1)  J  =  k(n-l ) (n-2)  j  t(y) J  S (x) {T(y)-T(x) }n  3t { (k+1 )y-kx}dxdy 

0  0 


>  0. 


Now 


t(y)  - 


y  >  0 


0  elsewhere 


s(x)  = 


p-1  "1 

1  -x 

p^x  e  ,  x  >  0 

0  elsewhere 


T  (b)  = 


1-e 

0 


-b' 


,  b  >  0 
elsewhere 


S(b)  = 


( 


V1 

1-e  ,  b  >  0 


0  elsewhere 


Then  the  left  hand  side  of  expression  (5.2.1)  becomes 


p  y  T  p  p  n-  3 

-x  N  r ,,  ~Y  \  -x 


(5.2.2)  J  =  k(n-l ) (n-2 )  J  py^e  7  j  (1-e  x  )  {(1-e  7  )-(l-e  X  )} 

0  0 


(continued) 
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X  nUk+Dy-kx^'V1  (k+1)y'kx)  Vdy. 


J 


5 

>  k(n-l)(n-2)J 

3 


_i  _  Tl  2 

Tiy11  e  7  J  (1-e 
1 


^1  r|  r\  n-  3 

~X  )  {(l-e_7  )-(l-e_X  ) } 


Ti{(k+l)y-kx}^~h~^  (k+1  )y-kx]  ^jxdy . 

Now  for  r)  >  0,  k  >  0,  r^>0,  3  £  y  _<  5,  1  x  2 


e 


>  0 


Tli 


(1-e  X  )  _>  (1-e  ^  )  =  l-  ^->0  independently  of  >  0, 


T]  P  0  h  Tl  o  0  Tl  o  T3  o 

r/i  _y  n  /i  _x  N-,n-3  (  -x  -y  ,n-3  v  r  -2  -3  ,n-3  v  n 

{(1-e  3  )  -  (1-e  )}  ={e  -  e  '  }  >{e  -e  }  >0 


{(k+l)y  -kx}11"1 


{ (k+1  )y-kx  }T*  ((k+1  )3-2k  }r|  =  (k+3)11 

{(k+1 )y-kx }  -  {(k+l)5-k  }  {4k+5  } 


e~[  (k+1  )y-kx] 11  >  e-[(k+l)5-k]T1  =  e-[4k+5]T1  >  Q 


Therefore 


J  >  kp^(n-l ) (n-2) 


11  -3p,n-3 

-e  } 


x 


(k+3)"1 

4k+5 

e-[4k+5] 


Tl 

2 


>  0  for  all  k  >  0,  n  >  2  ,  p,  p^  >  0. 
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Thus  the  exchangeable  model  with  (at  most)  one  spurious  observation 
based  on  the  Weibull  family  of  distributions  indexed  by  shape  parameter 
p  is  outlier-resistant  completely  on  the  right. 


5.3  Detection  of  Outliers 


Consider  a  situation  where  n-1  observations  are  distributed  as 
WEI(l,r),0)  and  one  is  distributed  as  WEI(1  ,k*ri,0) ,  0  <  k  <_  1  and  we 

have  the  exchangeable  model.  If  k*  <  1  we  have  exactly  one  spurious 

dG  k*n-ri  -(x^^-x^) 

observation.  Then  Y(x)  =  —  =  k*x  e  is  not  monotone. 

dr 

Letting  u(r;n,p,k*)  =  P(X^^  is  the  spurious  observation),  we  have 


.  00  hi  T)  14-1  k*p 

,  t  jl  \  rn- 1\(  ,,  -x  \r-l,  -x  \  n-r.  .  k*p-l  -x  , 

u(r;n,p,k*)  =  (  ,  JJ  (1-e  )  (e  )  k*rix  e  dx 

0 


-x 


k*p 


and,  setting  y  =  e  ,  we  obtain 


.  ,  m-1  >  i  r-i  ~(ln  1/y)  r  1/y)  >n-r, 

u(r;n,k*)  =  [  .  JJ  (l-e  J  (e  J  dy 

0 


which  is  independent  of  p.  Therefore,  without  loss  of  generality 
(w.l.o.g.),  we  may  take  p  =  1/k*,  giving 


,  °°  1/k*  .  1/k* 

,  *  j* \  rn-lw  -x  'vr-l/  -x  \n-r  -x, 

u(r ;n,k*)  =  (  ,  JJ  (l-e  J  (e  J  e  dx 

0 


and  in  particular 


' 
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1  -v1/k*  -v/(n-l)k*  . 

- t*  J  e  e  dv 

(n-ir  0 


^L(u  =  1/k*,  p  =  - “T"*) 


(n-1) 


(n-1 ) 


k* 


oo  k 

where  L(u,p)  =  J  e  V  e  ^Vdv  is  the  Laplace  transform  of  e 

0 


-v 


u 


(Re  u  >  0). 


Now  u(l;n,k*)  =  J  ( 

0 


00  1/k*  n 

-x  ^n-1  -x 
e 


)n  e  Xdx  ,  n  >  1 


=  j  tlf.  f  xme-(r,‘1)xl/k*dx 

L  m!  J 


m=0 


0 


1  /k* 

and  setting  t  =  (n-l)x  we  obtain 


u(l;n,k*) 


OO 


I 

m=0 


(~l)m 

m! 


00 


J 

t=0 


e"tk*tk*(m+1)"1 

(n_i)k*(m+l) 


y  (-l)\*r{k*(iiH-l)} 

n  •  r  i  \k*( m+1 ) 

m=0  m!(n-l) 


(5.3.1) 


oo  m 

y  (-1)  r{k*(m+l)+l} 

n  /  ,  i  m /  i  >k*(m+l ) 

m=0  (m+1)! (n-1) 


For  computational  accuracy  to  the  third  decimal  place,  we  would  require 
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q  terms  where  q  is  chosen  such  that 


r{k*(q+l)+l } 


k*(q+l)  — 


<  0.0005. 


(q+l) ! (n-l) 

(Apostol  1964).  For  fixed  n,  as  k*  ->  °°  (heterogeneity  increases)  q 
increases;  as  k*  ->  1  (heterogeneity  decreases)  q  decreases. 

Table  5.3.1  shows  q  values  for  some  selected  (n,k*)  values. 


Table  5.3.1  Number  of  terms  required  to  calculate  u(l;n,k*) 
accurate  to  third  decimal 


5 

Sample 

10 

size 

20 

n 

50 

Coefficient 

1/4 

5 

5 

4 

4 

of  spuriosity 

1/2 

5 

4 

3 

3 

k* 

3/4 

5 

3 

3 

2 

1 

5 

3 

2 

1 

We  can  also  develop  a  recursive  formula  for  u(r;n,k*)  in  terms  of 
u(l;j,k*),  j  <  n.  From  equation  5.3.1  we  have 


u(l ;n,k*) 


y  (-1 )mr{k*(m+l )+l } 

n  /  ,n,/  .  v  k*(m+l  ) 

m=0  (m+1 ) ! (n-l ) 


n  >  1 


and,  for  the  case  of  n  =  1,  we  may  define  u(l;l,k*)  =  1.  For  the 
special  case  of  r  =  n  we  may  write 

°°  _  1/k*  n-l 

u(n;n,k*)  =  J  (l-e  X  )  e  Xdx 
0 


• 
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00  n-1  1  .  1/k* 

j  l  C')(- 1)V*  )VXdx 

0  j=0  J 


n-1  1  .  00  1/k* 

1  +  I  (  l1  K-l) JJ  e-JX  e'Xdx 


j=l 


0 


n-1  .  °°  ,  .  Nm  .  1/k* 

I  l  *V^  dx 

j=l  J  m=0  m-  0 


,  "y1  rn-lv  j  y  (-l)m  k*r{k*(m+l ) } 

Jl  j  Jo  m!  jk*(n+1) 


,  .  "J  fn-lv  .j  y  (-l)m  £  {k*(m+l  )+l ) 

jti  3  Jo  (m+1)!  .k*(m+l ) 


n-1 

=1+1  (  "  )(-l)Ju(l;j+l,k*) 

j=l  J 

n-1 

=  1  )(-l)Ju(l; j+l,k*)  where  u(l;l,k*)  =  1  by 

j=0  J 

definition. 


dx 


£ln  (r:x)(-i)3[  i  k*(m+1) 

j=0  3  n=0  (m+l)!(n-r+j)  {  3 


(-l)mr{k*(m+l)+l} 


] 


(t-})  l  (ri1)(-l)^u(l;n-r+j+l,k*) 
j=0  J 


Thus  we  may  write 


l  for  r  =  1,  n  >  2 


m=0  (m+l)!(n-l) 


k*(m+l ) 


u(r;n,k*)  =  -«  1  if  r  =  1,  n  =  1 

i  r_1  i 

'n-1 ^  v  ( r-1 


l  (r"  )(-l)Ju(l;n-r+j+l,k*)  for  r  =  2 
j=0  J 


0 


otherwise 


(See  Appendix  IV).  Consider  now 


lim  u(r;n,k*) 
k*+l 


CD 

1/k*  . 

j  lim 

(  -x  -\n-l  -x ,  . 

(e  J  e  dx  ,  r  =  1 

0  k*-»l 

CO 

1/k*  1 

J  lim 

f.  -x  %n-l  -x, 

(1-e  J  e  dx  ,  r  =  n 

0  k*->l 

CO 

1/k*  ,  1/k* 

("I  )J 

lim  (l-e  X  )r  (e  X 

r  1  0 

k*+l 

.n-r  -x 
)  e  dx 


r  2 ,  . . . ,  n- 1 


otherwise 


Thus  we  find  that 
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-xn,  , 

e  dx  ,  r  =  1 


0 


11m  u(r;n,k*)  = 
k*+l 


I  -xn-1  -x 

J  (1-e  )  e  dx  ,  r  =  n 

0 


("  i)  J  a-e'x) 

0 


-xNr-l  -x(n-r+l )  0  , 

'  e  dx  ,  r  =  2 ,  . . . ,n-l 


VP 


,  otherwise 


1/n  , 

1/n  , 


r  =  1 
r  =  n 


fn-l >  r(r)r(n-r+l ) 
lr-l  J  r(n+l 

0  ,  otherwise 


=  1/n  ,  r  =  2,3,...,n-l 


The  interpretation  of  this  result  is  that,  as  k*  -*  1  (k*  =  1  is  the 

case  of  homogeneous  data),  the  heterogeneity  of  the  data  is  decreasing 

and  the  probability  of  any  order  statistic  being  spurious  approaches 

1/n,  i.e.  all  equally  likely. 

Now  consider  lim  u(r;n,k*).  We  have 
k*-K) 


°°  1/k*  .. 

J  lim  (e  X  )n  e  Xdx  ,  r 
0  k*->0 


=  1 


lim  u(r;n,k*) 
k*+0 


<»  _  1/k*  1  _ 

]  lim  (l-e  X  )n  e  Xdx  ,  r  =  n 
0  k*->0 

.  ®  1/k*  .  1/k* 

fn-l'ii  ,  .  (.  -x  >r-l,  -x  >n-r  -x 

L  ,  J  J  lim  U“e  J  (e  )  e  dx  , 

0  k*+0 

r  2,..., n- 1 


otherwise 


Now 
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00  _  1/k* 

lim  u(l;n,k*)  =  lim  J  (e  X  )n  e  Xdx 
k*->0  k*-*0  0 


=  lim  £ 
k*->0  m=0 


(-Dm 

m! 


co  1/k* 

J  xm(e'x  )n_1dx 
0 


llm  l  (-l)mr{k*(m+l)) 

k*+0  111=0  m!  |p-|  (n-l)k  (m 


n  >  1 


v  (-i)V(i+i)r{k*(®fi)} 
lim  ^  - k*Vm+ll - 

k*-K)  m=0  (m+l)!(n-l)  ^  ' 

-i  .  v  (-I)mr{k*(m+1)+1} 

^"lm  ^  k*Cm+l') 

k*+0  m=0  (m+l)!(n-l)  V  ; 

°°  in  0  13 

V  (-D  -  (-D  .  (-D  ,  (-1)  . 

1  (m+1 ) !  1!  2!  3! 

m=U 


=  1  -  -  =  .63  . 
e 


For  r  =  n 


lim  u(n;n,k*)  =  lim  J  (l-e 
k*->0  k*-K)  0 


co  1/k*  n-1 

-x  -x 


)  e  dx 


co  n-1  .  .  1/k*  j 

=  lim  j  l  ("-1) (-l)J(e-X  )  e'xdx 

k*-K)  0  j=0  J 
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=  llm  Y  (n‘1)(-Dj 

k*+0  j=0  J 


00  .  1/k* 

f  e’JX  e‘Xdx 
0 


11m 

k*-K) 


n-1 


1 

j=i 


(n?)(-i)3 


00  .  1/k*  )  co 

J  e’JX  e"Xdx  +  j  e"Xdx 
0  )  0 


lim 

k*-K) 


n-1 

I 


j=i 


J  m=0 


(-l)mr{k*(m+l)+l}-|j 

(m+l)!jk*(m+1)  i 


+  1 


n-1 


)<-Dj[  l 

m=0 


.m. 


(-D  r(i) 

(m+1 ) ! 


]  +  1 


.37 


For  r  =  2,3,..., n-1 
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lim  u(r;n,k*)  =  lim  (n_j)  J  (l-e 
k*+0  k*+0  r  1  0 


-x 


1/k*  r-1  1/k*  n-r 


)  (e  X  )  e  Xdx 


n-1  00  r-1  .  1/k*  j+n-r 

llm  ( r_! )  i  1  (r4  X  )  e  Xdx 

1  0  j=0  J 


k*+0 


li® 

k*+0  j=0  J 


“  /  i  Nm  00  /  • ,  \  1/k* 

l  I  xme-0-Hi-r)x  dx, 

m=0  m'  0 


i  r"l  i 
n-1  >  v  ( r-1  ,o 


=  iim  (,: )  i  (7)(-d 

k*-K)  j=0  J 


y  (-1)  i{k*(m+l) } 
n  i j 1  | ,  . ,  N k*(m+l ) 

1=0  m!  |^*|  ( J+n_r) 


-  r]n  r1) c-n3  uJ  l 

r  1  j=0  3  k*-»0  (m=0  (m+l)!(  j+n-r )k  1 


r^n  (r_1i(-D j)  ?  <-i>m 
lr-iy0  1  j  K  15  j miQ  (-fl)l 


( j=o  J 


-  ^Ic-i+D 
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=  0  . 


Thus,  as  k*  -*■  0  (i.e.  one  observation  with  Weibull  shape  parameter 
much  smaller  than  the  rest),  we  have 


,  otherwise 

In  this  situation,  heterogeneity  is  more  evident  and  the  probability 
that  the  first  order  statistic  is  spurious  approaches  1  “  ~  =  *63,  the 
probability  that  the  last  order  statistic  is  spurious  approaches 
=  .37,  and  the  probability  that  any  other  statistic  is  spurious 
approaches  zero. 


lim  u(r;n,k*)  = 
k*->0 


i- 1 
e 

1_ 

e 

0 
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5.4  Estimation  for  Weibull  parameters 
5.4.1  Standard  Estimators 

Complete  sufficient  statistics  do  not  exist  for  the  Weibull 
distribution.  Much  of  the  estimation  for  the  Weibull  distribution 
falls  into  one  of  three  categories: 

i)  For  M.L.E’s,  the  distributional  results  are  not 
mathematically  tractable  hence  Monte  Carlo  simulation  is  needed  to 
tabulate  percentage  points.  While  the  M.L.E.’s  have  optimality 
properties  as  n  increases,  they  are  not  in  closed  form  and  require 
computer  solution  to  obtain 


t)  = 


n 


£  X.  Jin  x 
i=l 


n 


i  n 

- I  Jin  X. 

n  .  L.  l 


l  < 


i=l 


i=l 


n 


-1 


and  0^  =  — 


l 


n 


/:  no 

For  n  >_  100,  f)  is  asymptotically  N(r|,  - — — - )  [see  Cohen  (1965)  and 

Harter  and  Moore  (1965)].  Inference  procedures  based  on  M.L.E.’s 


depend  on  pivotal  quantities 


—  and 

r\ 


[see  Thoman,  Bain  and 


Antle  (1970) ] . 

ii)  Best  Linear  Estimators  have  properties  similar  to  M.L.E.’s. 

Mann  and  Fertig  (1973)  considered  BLIE  for  the  parameter  b  of  EV  . 

r 

The  BLIE  b  T  =  T  c.  X,.N  is  a  linear  function  of  the  BLUE  and 

BLIE  ^  i,r,n  (l) 

r 


a  maximal  invariant  since 


=  0.  (r  is  the  size  of  the 
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BLIE 

censored  sample).  The  distribution  of  — - —  is  independent  of  b 

b 

and  £.  Mann  and  Fertig  (1977)  also  took  Hassanein’s  (1972) 

k 

asymptotically  unbiased  optimum  estimators  b,  =  Y  b.  ,  X,  N  , 

k  i^1  i,k  (n1) 

n.  =  [y  n]  +  1  and  obtained,  for  small  samples,  the  unbiased 

X  1 ,  tt 


b*  = 


k,n 


where  b 


.  =  E[  T  b.  .Z,  J,  where  Z.  ~EV_(0,1)  and 

k,n  i_=\  1,k  1  I  cb  K2 / 


df 


and  the  BLIE  b* 


b* 


BLIE  1+c 


k,n 


where  c.  =  Varf-^— )  . 
k,n  vb  1 


iii)  Simple  estimators,  such  as  good  linear  unbiased  estimators 
(GLUE’s)  and  modified  GLUE’s,  are  identical  to  BLUE’s  if  r  =  2  and 
similar  but  not  identical  if  r  >  2,  where  r  is  the  size  of  the 
censored  sample.  They  are  essentially  equivalent  to  M.L.E.’s  for 
censored  sampling  but  are  slightly  less  efficient  for  the  complete 
sample  case.  Bain  (1973)  suggested  using  the  approximate  BLUE 


r-1 


^B  nk 


ly(i)-y(r)l  (r_1)y^r)  -  y(1) 


r,n 


nk 


r,n 


1  n_1 

where  y.  ~EVT(£,b)  and  k  =  —  —  E  {  Y  (w.-w  )}  where  w.  are 
J  i  I  ’  r,n  n  1  “  l  r  J  l 


■ 


Ill 


the  order  statistics  of  EV_(0,1).  Then  2nk  ~  Xo  (if 

I  r,n  b  2nk  df 

r,n 

nk  -1 

r/n  <  1/2)  and  r|  =  — — y- — -r—  is  approximately  unbiased  for  t).  For 

d  nk  b„ 
r,n  B 

A 

complete  samples  (i.e.  r  =  n)  ,  b  has  zero  asymptotic  relative 

B 


efficiency. 


„  £  |y(i)'y(s) 

Englehardt  and  Bain  (1973)  suggested  b^  =  - — -  as  an 


improvement  where 


s  >  r  ,n 


s  = 


n 

n-1 

[.892  n]  +  1  , 
r 


2  £  n  £  15 
16  £  n  £  24  . 

n  >_  25 
r  <  . 9n 


and  r  >  .9n 


and  h  -r—  ~  X?  where  h  = 

b  h  df 


b  * 

Var(^) 


For  —  ->  0  (very  heavy  censoring)  as  n  ->■  ®  ,  b_  and  M.L.E. 
n  B 

A  A 

appear  to  agree.  Compared  to  b  ,  b  has  relative  efficiency  ranging 

B  s 

from  1  (n=2)  to  0.82  (n=25). 


Mann  and  Fertig  (1975)  converted  Bain’s  estimator  to  an 

A 

b 


approximate  BLIE.  For  n  >_  20, 


r,n 


C2/X  df 


where 


r,n 


JJy<s>  7(1)1  ~  bs 

b  =  - : -  and  approximate  BLIE  b._  =  7-7- 

s  nk  MF  1+x 

s,r,n 


where 


r,n 


c r> 
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>  =  Var  (t-— )  =  -r - - 

r,n  kb  1  nk 


s,r,n 


X  b 

.  Then  MSE(b  )  =  ■  —  .  For  heavy 

Mr  IT  X 


r,n 


censoring  (r  «  n) ,  b,„  is  closer  to  b,„ _  then  b  . 

MF  MLE  s 

In  terms  of  the  Weibull  shape  parameter  p,  af)  =  a  [■■J' ■  )  is 

^MF 

2  1 

unbiased  for  p  if  a  =  ■.  ■_  and  dp  =  df - 1  has  minimum  MSE  if 

h+2  v  ~  1 


MF 


j  h“4  u  h+2  i  ,  „  , .  i_  2 

d  =  i where  — —  =  1  +  X  (i.e.  h  =  — - r z — nr 

h-2  h  r,n  Var [b  /bj 


)• 


Englehardt  and  Bain  (1977a)  proposed,  as  a  new  simple  unbiased 

A 

estimator  for  complete  samples  to  replace  b  , 

s 


n 


+ « i-L7^ 


EB 


nk 


n 


where  s  =  [qn],  0  <  q  <  1  is  chosen  to  minimize  the  asymptotic 
variance  of  .  This  results  in 

DD 


and 


s  =  [ . 84n ] 


n 


kn  =  "(i)  +  ™  ±J+1  w(i)^/n  >  “i  ~  EV(0>1) 


=  -[e  l  W  +  sy]/(n-s) 
i=l  K  J 


' 
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and  y  is  Euler's  constant.  Then 


b  V  +  (— )u 

n  Var(-&2.  -  -2- 


sir 


n-s *  1  s+1 


(n-s)k 


n 


n 


n 


nit 


where  =  Var  w(±)),  Ug+1  =  Var(^  w(1)),  VarfJ^  w(1))  =  6 


To  obtain  an  approximate  BLIE  (or  MLE)  we  may  use  b*  = 

EB 


EB 


1+Var  (— r— ) 


Three  other  simple  estimators  have  been  suggested.  Menon  (1963) 

X  T) 

used  the  transformation  Z  =  An(-^-)  where  X  ~WEI(9,p,0)  and 

6s^  1/2 

/n  s  y  r-|  y  x 

In  X  ~EVI(C,b).  He  suggested  b^  =  ( - ^ - ) 

% 


/\ 

and  ^MN  =  ^  9 


n 


n 


9mn  e 


I  fa  x  .  2 

i=l  ,  ~  *  1  f  l.lri  \  , 

—  +  .  Then  =  g-  ~  N(r|,  — — )  and 

M 

2  2 

1  •  2  b  0 

~  N(  9,  — ^ - ).  Dubey  (1967)  suggested  estimating  t)  by  using 

percentiles.  Based  upon  two  percentiles, 


Anf-Antl-p^)  }-Jbi{-ln(l-p2)  } 
^  In  y  -  An  y 


where  p^  =  17%  and  p£  =  97%  minimize  Var(fi^)  and 


yP  = 


x,  N  if  np  is  an  integer 
(np) 

x r  if  np  is  not  an  integer. 

[np]+l 


114 


*  n  -  „r  ri2(. 91627479)^ 

Asymptotically  ~  N(p,  — — — - — J  . 

Murthy  and  Swartz  (1975)  used  an  approach  similar  to  Dubey  but  used  the 
EV^.  distribution  and  obtained  an  estimator  of  b  based  upon  two  order 
statistics 


bMS  “  T(j)  '  *nT(«}  B(N’X’j) 


where  B(N,l,j)  = 


and  Y  = 


jn  T  —  Sn  T 

(j)  (JO 


This  is  MVUE  and 


2E(Y)  2b 

has  relative  efficiency  w.r.t.  the  Cramer-Rao  lower  bound  (CRLB) 


approaching  70%.  The  following  Table  5.4.1  gives  the  optimum  values 


for  j  and  i. 


Table  5.4.1 


Optimal  Order  Statistics  for  Murthy-Swartz  Estimator  of  b  =  — 


n 

j 

1 

2  <_  n  <  5 

n 

1 

6  £  n  £  10 

n 

2 

11  1  n  1  16 

n 

3 

17  £  n  <  23 

n 

4 

24  i  n  £  26 

n 

5 

5.4.2  Estimators  suggested  for  use. 


Under  the  exchangeable  model  with  a  shape  change,  we  assume  n-1 


115 


rpl  r\ 

observations  are  governed  by  f(x;9,r])  =  — -  exp  {-  },  x  >  0, 


t|,  9  >  0  and  one  observation  is  governed  by  f(x;9,k*r|),  0  <  k*  <_  1. 

If  0  <  k*  <  1,  we  have  exactly  one  spurious  observation.  We  are 

interested  in  estimating  r|  (or  in  the  EV^.  case,  b  =  -^) . 

From  Table  5.4.1,  for  2  <_  n  <_  26,  the  optimal  form  of  the 

Murthy-Swartz  estimator  uses  the  largest  order  statistic.  Dubey’s 

estimator  involves  one  or  both  of  the  smallest  and  largest  order 

statistics  for  2  <  n  <  33.  Also  b  ,  2  <  n  <  15,  b  and  b  all 

—  —  S  —  —  B  EB 

involve  these  two  highly  suspect  values.  Menon's  estimator 


MN 


f6sL  x 


has 


MSE<bMN>  ’  Var<b,V  +  {Blas(bMN)}' 


1,1  b  ■  +  b20  (It)  +  {b0  (— -)  }2 


n 


n 


n 


under  the  homogeneous  model.  However,  with  the  above  exchangeable 
model , 


E  (b  )  =  —  \JVax —  +  tsee  Cramer  (1946)>  27.7.1] 

het^  MN'  n  *■  In  x 


=  ±  {Jj'i  b2( s=k)  +  4(?t)2  -  +  boC— H 

x  *  6  v  n  7  6  ^k*  7  n  vn  7  J 


' 
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+  bo(-) 


and  as  k*  ■>  0,  bias  ■*  a3. 

Also  MSE  (^)  =  v.r  (fi^)  +  (BlasCb^)}2  and 


6  rV^ 


’"wCO  -  \  l4SIT^+  tsee  Cramer  (1946),  27.7.2] 

it  ^2  n 


it 


,  L  n-l+  -h-  4(2.404)Y(n-l+  -tr)(n-l+  if) 

3t  fe(  .  k*  )  + - - 


2  it' 


2  L15  ^  n 


n 


x2Y2(n-1+  1  )(n.1+  1_)2  3YVl+lf)4  ,4  12 

l  k*  _  k*  1 

+  5  A  0  J 


n 


n 


+b2o  (i») 


n 


36n 


(see  Appendix  V). 

Taking  lim  MSE,  (h._T)  ,  we  find  MSE(b  )  ->  00  as  k*  >  0.  Thus 
°  ,  .  -  bet  MN  MN 

k*-K) 

this  estimation  technique  seems  poor  if  an  outlier  is  present.  Also 
2 

s  would  be  strongly  affected  by  the  presence  of  an  outlier 

aJI  X 

occurring  at  or  at  X(n)  an<^  these  are  precisely  the  values 

where  u(r;n,k*)  is  largest  as  k*  ■*  0.  We  might  obtain  a  more  robust 


estimator  by  using  trimming,  winsorization  or  semi-winsorization.  Thus 


we  might  consider 
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i) 


RA 


6  a  1/2 

-4) 


and 


It 


3RA 


^  +  ^RA 


ii) 


6°W  1/2  - 

'  2''  and  ’RW  4  +  -bRW 


6a\  1/2 

1115  bRS  =  (“2 )  and  =  +  YbRS 

% 

1  a2  a2  a2  a  a  a 

where  =  In  ,  b  =  -  ,  ^  =  In  0  and  o^,  o^,  ,  p.^,  p^  and  p,. 


are  as  defined  in  Appendix  I  and  y  is  Euler’s  constant. 


Table  5.4.2  Computation  of  means  and  variances  and  winsorized  means  and  variances  for  25  samples  of  size 
five  based  on  EV  (0,1 /n)  with  one  spurious  observation  present  (see  Appendix  1) 
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We  could  seek  M.L.E's  of  9  and  v\  under  the  exchangeable  model 


where 
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The  above  equations  appear  mathematically  intractable. 
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Appendix  I:  The  A-,  W-,  and  S-Rules  and  Premium/Protection. 


We  adopt  the  following  notation: 


zi  =  7i  -  y 


y  = 


n 

£  7± 

i=l 


n 


^  —  9 

I  (y,-y> 

2  i=l 

S  =  - = - 

n-1 


n 

l  y(1) 
1=2  •* 
'(1)  n-1 


n-1 

I  y 

1=1 


n 


(i) 


(  I  ym^nry 


Y(n)  n-1  y(j,k) 


i=l 


(i)  J(j)  J(k) 


n 


n 


2  i=2 


I  (y (1)-?(i)) 


2 


(1)  n-2 


(n)  n-2 


2 

s  /  . 


ZTT  (  l  vJ2  +  vJ2  -  ( 


y^-y/.,  ..0  }• 


(j,k)  (n-1)  w(l)  '  w(j)  7  ( j  »k)  1  w  (k)  (  j,k) 


Anscombe’s  Rule  (A-Rule): 


1  y 

if 

|z(n) 

|  <  Cs 

and 

^d)1 

<  Cs 

II 

?(D 

if 

lZ(l) 

1  >  cs 

and 

|zci)! 

>  |z(n) 

1%) 

if 

|2(n) 

1  >  Cs 

and 

|2(n)l 

>  lZ(l) 

Ds2 

if 

2 

Z(D 

<  Ks2 

and 

Z(n>  < 

Ks2 

2  i 

C A. 

I  DS(1) 

if 

2 

z(l) 

>  Ks2 

and 

Z(D  > 

2 

Z(n) 

1  Ds2  . 
(n) 

if 

2 

Z(n) 

>  Ks2 

and 

2  s 

z  > 

(n)  - 

2 

Z(D 
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Winsorization  (W-Rule): 


if  z,, N  <  Cs  and  z,  s  <  Cs 
1  (1)'  1  (n) 1 


7  (2 , 1 )  if  I  Z  (1 )  >  -  Cs  and  i Z  (1 )  I  >  I  2  (n) 


y,  ,  x  if  z  ,  N  >  Cs  and  z,  N  >  z,1N 
(n-1 ,n)  1  (n) '  -  1  (n) 1  1  (1) 


2 

°w  = 


Ds 


if  z}..  v  <  Ks2  and  z2  N  <  Ks' 
(1)  (n) 


D  max [ s ^ , *,  s2  .  if  z2  >  Ks2  and  z2  N  >  z2 

L  (2,1)’  (n, 1 ) J  (1)  —  (1)  (n) 

2 


2 

!o 

2 


D  maxfs/.  N,  st  1  N  1  if  z2  N  >  Ks2  and  z2  N  >  z'),  N 
L  (1 ,n)  (n-1 ,n) J  (n)  -  (n)  (1) 


Semi-Winsorization  (S-Rule): 


f  _ 

y 


= 


if  z,-N  <  Cs  and  \z,  N  <  Cs 
1  (l)1  1  (n)  1 


(n-l)y^+(y-Cs) 


n 


(n-l)y^+(y+Cs) 


n 


if 

lz(i)l 

>  Cs 

and 

lzd)l 

> 

|z(n)i 

if 

|z(n)l 

>  Cs 

and 

|z(n)l 

> 

|z(1)!  • 

= 


r  ds' 

D 


2  2  2  2 
if  z/1N  <  Ks  and  z,  .  <  Ks 
(1)  (n) 


2 

(1) 

+  Ks2} 

if 

2 

Z(D 

> 

Ks2 

and 

2 

(n) 

+  Ks2} 

if 

2 

Z(n) 

> 

Ks2 

and 

2  n  2 
'(1)  >  Z(n) 


(n)  ‘(1) 
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If  we  adopt  the  premium-protection  approach  suggested  by  Anscombe 
(1960)  we  define 


Premium 


MSE(new  estimator)-MSE(old  estimator) 
MSE(old  estimator) 


and 


Protection  = 


assuming  homogeneous  data 


MSE(old  estimator)-MSE(new  estimator) 
MSE(old  estimator) 


assuming  spurious  values (s)  are  present. 


Appendix  II:  Asymptotic  approximation  for  f(v+l) 


Consider  T(v+1)  =  J  e  u  du.  Setting  u  =  vt,  du  =  vdt  and 

0 

00  00 

r(  v+l)  =  J  e  Vt(vt)Vvdt  =  V^1]  e t+'^n  t-^dt.  This  last  integral  is 
0  0 


of  the  form  j  #(t)e 
0 


vh(t) 


dt  where  v  is  a  large  positive  constant 


and  i?(t)  and  h(t)  are  real  and  continuous  in  [0,®).  Now 
h(t)=-t  +  Jhit  which  has  a  single  maximum  at  t  =  1  with 
h’(l)  =  0,  h"(l)  =  -1.  Applying  the  Laplace  approximation  to  the  two 
intervals  0  <  t  <  1  and  1  <  t  <  00  we  obtain  the  following: 


j"  ev(-t+!n  t)dt  _  j1  ev(-t+in  t)dt  +  j”  e  v(-t+Jln  t)dt  _ 
0  0  1 


Since  h(t)  has  a  maximum  at  t  =  1,  the  main  contributions  to  the  two 

integrals  on  the  right  occur  in  the  neighborhood  of  t  =  1.  We  may 

.  }  v(-t+Jln  t),  £  v(-t+jfrit)Jfc  ,  }  v(-t+£n  t).  , 

rewrite  J  e  v  dt  =  J  e  dt  +  J  e  'dt  and 


0 

v(-t+ln  t) 


dt  =  J 


1+6 


0 

v(-t+ln  t) 


dt  +  J 


1-e 

v(-t+Jln  t) 


dt. 


Setting 


1  1  1+6 

2 

x  =h(l)-h(t),  2xdx  =  -h’(t)dt  and  expanding  h(t)  in  a  Taylor’s 

2 

Series,  h(t)  =  h(l)  +  h’(l)(t-l)  +  h"(5)  (t^y)  ,  1  <  l  <  1+6,  we 

2 

obtain  x^  =  h(l)  -  h(t)  =  -h"(£)  ^ —  and 
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h’(t)  =  h"(5)(t-l)  . 


Now 


2x 

h7^) 


h"(5) 


(t-1) 


h’*(0(t-l) 


2 


-1 

J-  \  h"<5) 

-1 

y-  j-  h-ci> 


Thus 


1+6 


J 

t=l 


v(-t+£n  t) 


dt  =  J 


x=0 


v(-l-x2) f  ~2x  1 

V(t)J 


dx 


dx 


i  2 

-V  sr-  ,  ~VX 

=  e  /2  J  e 
0 


dx  . 


i  2  •  °°  2 

But  J  e  dx  =  J  e 
0  0 


dx  since  the  major  contribution  to  this 
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latter  integral  occurs  in  the  neighborhood  of  x  =  0.  Thus 


1+6  /  ,  n  x  •  00  2 

,  v(-t+ln  t).  ,  - vx 

J  e  dt  =  e  /2  J  e  dx 

1  0 


Similarly  we  can  show 


1 


J 

1-e 


v(-t+ln 

e 


and  thus 


r(v+l)  = 


v+1  r  f  v(-t+!n  t) 

v  IJ  e 

l-£ 


1+6 


dt  +  J 

1 


v(-t+£n 

e 


v+1 

=  V 


V  V 
V 


/2-jtv 


CO 


Appendix  III:  Evaluation  of  j  <})(x)<I>( vx)dx . 

— oo 


00 

Theorem:  J  <j)(x)<I>(  vx)dx  =  1/2. 

—  OO 


Proof:  Let  f(v)  =  ]  <t>(x)<3K  vx)dx . 

—CD 

CD 

Then  f'(v)  =  J  <j)(x)  <})(  vx)xdx 


i  00  i  ,  /0,  2,  2  2. 

1  I  1  -l/2(x  +v  x  ) 

—  nzz  J  I  ^ 

/2n  -<»  /2x 


xdx 


4- j'S  e“  2  (1^  )xd* 

/2n^l+v2  /2tc 


and  setting 


u  =  xV’  '  2 


1+v  we  obtain 


.  1  00  -z2/2 

f’(v)  =  -3-  - - J  z  g-==—  dz 

/271  Vl+v2  ~°°  /27t 

=  0. 


Then  f(v)  =  c  for  all  v. 

In  particular,  if  v  =  1, 

=  1/2. 


00 

f(l)  =  c  =  J  <Kx)3>(x)dx  =  — 


Thus  f(v)  =  1/2  for  all  v. 

OO 

i.e.  j  4)(x)<?E>(  vx)dx  =  1/2. 
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Appendix  IV 

Tables  of  u(r;n,k*)  =  P(X^^  is  spurious  in  a  sample  of  n|k*) 
where  the  exchangeable  model  of  Weibulls  with  at  most  one  outlier 
(involving  shape  change)  is  assumed. 
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* 

K  s  .0 1050 
N  s  50 


T1 


1 

1 . OOOOO 

2 

.62989 

3701  1 

3 

.62721 

00535  .36743 

4 

.  62565 

00470  . 00334  .36632 

5 

.62454 

00444  .00273  .00263 

36566 

6 

.62368 

00431  .00249  .00205 

00226  .36521 

7 

.62297  . 

00422  00237  .00182 

00171  .00202  .36487 

6 

.62238 

00417  .00229  .00171 

00146  .00151  .00186 

3646  1 

9 

.62186 

00412  .00223  .00164 

00137  .00127  .00138 

00  173 

36439 

1  O 

.62141 

00409  .00219  .00159 

00131  .00120  OOI 10 

00  128 

00163  .36421 

1  1 

.62100 

00407  .00216  .00154 

00127  OOI 10  OOIOO 

OOIIO 

00120  .00154 

. 36406 

1  2 

.  6  2063 

00405  . 002 13  00151 

00123  OOI 10  . 00095 

0009  1 

00099  OOIIO 

00  146 

.36392 

1  3 

. 62030  . 

00403  . 002 11  .00146 

00120  OOIOO  00092 

00084 

00084  00095 

OOIIO 

00  140 

3638  1 

1  4 

.61999 

00402  .00210  00146 

00120  OOIOO  . 0009 1 

00081 

00075  .00079 

.  0009  2 

OOIOO 

.00134 

36370 

1  5 

.61970  . 

36  36  1 

00401  .00209  .00145 

OOI 10  .00098  .00089 

00080 

00072  .00069 

.  00076 

00089 

. 00099 

00  129 

1  6 

.61943 

.00125 

00400  .00208  .00144 

.  36353 

OOI 10  .00095  .00087 

00079 

0007 1  . 00064 

. 00065 

00075 

.00086 

0009  5 

1  7 

.61918 
. 00090 

00399  .00207  .00143 

.00120  .36345 

OOIIO  .00093  .00084 

00078 

0007 1  . 00063 

. 00059 

0006  3 

.00073 

.  00064 

1  8 

.61895 

0008  1 

00398  .00206  .00142 

. 00086  .00120  . 36338 

OO 1 1 0  . 0009 1  . 0008 1 

00076 

00070  . 00063 

. 00057 

00055 

. 0006 1 

. 00072 

1  9 

. 6 1 873 
0007  1 

00397  .00205  .00142 

.00078  .00082  .00120 

OOIIO  .00089  .00079 

36332 

00074 

.00070  .00063 

. 00056 

.  00052 

. 00053 

0006  1 

20 

.61852 
. 00060 

00396  .00204  .00141 

.00070  .00075  .00078 

OOIIO  .00088  00077 

OOIIO  .36326 

00071 

00068  .00063 

. 00057 

. 0005 1 

.00048 

. 0005 1 

2  1 

.  6  1 832 
. 0005 1 

00396  . 00204  .00140 

00060  .00069  .00072 

OOIIO  . 00087  . 00075 

.  00075  OOIIO  .36320 

0006  9 

00066  .  0006  3 

. 00058 

. 0005 1 

00046 

. 00046 

22 

.61814 

. 00044 

.00395  00203  .00140 

00050  .00060  .00067 

OOIIO  .00087  .00074 

. 00069  . 00072  OOIIO 

00067 

.36315 

00064  . 00062 

00058 

00052 

. 00046 

. 00043 

23 

.61796 
. 00040 

.00395  .00202  .00139 

.00043  .00050  .00060 

OOIIO  .00087  .00073 

.00066  .00066  .00069 

. 00065 
.00110 

00062  .00061 
363  10 

. 00058 

. 00053 

. 00047 

. 00042 

24 

.61778 

00039 

00394  .00202  .00139 

.00038  .00042  00051 

OOIIO  .00087  .00072 

.00060  .00064  .00063 

0006  4 
. 00067 

00060  . 00059 

OOIIO  . 36306 

. 00057 

. 00053 

. 00046 

00042 

25 

.61762 
. 00039 

00394  .00201  .00138 

.00036  00037  .00042 

OOIIO  .00087  .00072 

.00051  .00059  .00062 

. 00062 
0006  1 

00058  .  00067 

. 00066  OOIIO 

00056 

3630  1 

00054 

.  00049 

00043 

26 

.  6  1746 
00039 

.00394  .00201  .00138 

.00036  .00034  .00037 

OOIIO  .00087  .00072 

.00043  .00051  .00058 

.  0006  2 
. 00060 

00056  .00055 

.00058  .00063 

.  00054 
OOIOO 

00053 

.36297 

. 00049 

. 00044 

27 

6  1731 
00040 

.00394  .00201  .00137 

.00036  .00033  .00033 

OOIIO  .00087  .00073 

.00036  .00043  .00052 

0006  1 

00058 

.00055  .00053 

.  00058  .  00056 

0005  3 
0006  1 

00052 

OOIOO 

.00050 

36293 

00046 

28 

.61717 
. 0004 1 

.00393  00200  .00136 

.00036  .00033  .00031 

OOIIO  . 00087  . 00073 

.00032  .00036  00044 

.  0006  1 

. 00052 

00054  . 0005 1 

.00056  .00055 

. 0005 1 
. 00053 

0005  1 

. 00060 

.00050 

OOIOO 

. 00046 
36289 

29 

6  1703 
. 00042 
.36266 

.00393  .00200  .00136 

.00037  .00033  .00031 

OOIIO  .00087  .00073 

.00030  .00032  .00037 

. 00062 
. 00044 

00053  . 00049 

00052  . 00055 

00049 

00053 

00049 
. 0005 1 

00049 

. 00058 

. 00047 

OOIOO 

30 

.61689 
. 00043 
OOIOO 

00393  .00200  .00136 

.00039  00034  00031 

36282 

OOIOO  .00087  .00073 

.00029  .00029  .00031 

. 00062 
.  00037 

.00053  .00048 

00045  00052 

00047 
. 00054 

0004  7 

. 0006 1 

00048 

00049 

. 00047 
.  00057 

3  1 

.61676 
. 00044 
. 00056 

.00393  .00200  00135 

00040  00035  .00031 

OOIOO  .36279 

OOIOO  .00087  .00074 

.00029  .00028  .00028 

. 00062 
. 00031 

00053  .00047 

. 00038  00046 

. 00045 

. 00052 

. 00045 
. 00052 

00047 

. 00049 

. 00046 

00046 

32 

.61663 

. 00044 
00046 

.00392  .00199  .00135 

.00041  .00036  .00032 

00055  OOIOO  .36275 

OOIOO  .00086  .00074 

.00029  .00027  .00027 

. 00063 
. 00028 

.00053  00047 

.00032  .00038 

. 00044 

00046 

. 00044 
0005  1 

00045 
. 00060 

00045 

00047 

33 

6  1651 
. 00044 

.  00045 

. 00392  .00199  .00134 

.00041  .00037  .00033 

.00045  .00055  OOIOO 

OOIOO  .00086  .00074 

.00029  .00027  .00026 

.  36272 

. 00063 
. 00026 

00054  . 00046 

00027  .00032 

. 00043 

. 00039 

. 00042 
. 00047 

00043 
. 00050 

00044 

00046 

34 

6  16  39 
. 00043 
. 00047 

. 00392  .00199  00134 

.00042  .00038  .00034 

.00043  .00043  .00054 

OOIOO  .00085  .00074 

.00030  .00027  .00026 

OOIOO  .36269 

. 00064 

. 00025 

00054  .00046 

.00025  .00027 

00042 

00033 

00040 

. 00040 

.00041 

. 00047 

.  00043 
00049 

35 

.61628 
. 00042 

0004  8 

. 00392  .00199  .00134 

.00042  .00039  .00035 

.00045  .00041  .00042 

OOIOO  .00085  .00074 

.00031  .00028  .00026 

.00053  OOIOO  .36266 

. 00064 
. 00024 

00054  . 00047 

.00024  00025 

. 0004 1 
00028 

00039 

00033 

. 00039 
0004  1 

. 0004 1 

00047 

36 

37 

38 

39 

40 

4  1 

42 

43 

44 

45 

46 

47 

46 

49 

50 
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6  16  16 

00392 

00199 

00134 

OOIOO 

00084 

00073 

00064 

00055 

00047 

0004  1 

00038 

00038 

. 0004 1 

. 00047 

. 0004 1 

. 00047 

.  00040 

. 00043 

. 00036 
. 00039 

. 00032 
.  0004  1 

. 00028 
. 00053 

. 00026 

.OOIOO 

. 00024  . 00023 

36263 

. 00023 

.  00024 

00028 

. 00034 

6  1606 

0039  1 

00199 

00134 

OOIOO 

00084 

00073 

00064 

00055 

00047 

0004  1 

00037 

00036 

. 00039 
00042 

.  00040 
. 00046 

00040 

. 00045 

. 00037 

. 0004 1 

.  00033 

.  00038 

. 00029 

. 0004 1 

.  00026 

. 00052 

. 00024 

OOIOO 

. 00023 

.  36260 

. 00022 

.  00022 

00024 

. 00028 

6  1595 

0039  1 

00  199 

00134 

OOIOO 

00083 

00073 

00064 

00056 

00048 

0004  1 

00037 

00035 

. 00038 
. 00036 

00039 
. 00043 

.  00040 

. 00046 

. 00038 
00044 

. 00035 
.  00039 

. 00030 
. 00037 

.  00027 
.  00040 

. 00025 
.  00052 

. 00023 
OOIOO 

. 00022 

.  36  258 

00022 

00022 

. 00024 

6  1  585 

0039  1 

00198 

00  134 

OOIOO 

00083 

00072 

0006  4 

00056 

00048 

0004  1 

00037 

00034 

. 00036 
. 00030 

.  00038 
.  00037 

. 00039 
. 00043 

. 00038 
. 00045 

.  00035 
. 00042 

. 00032 

. 00037 

. 00028 

. 00035 

.  00025 
. 00039 

. 00023 
. 0005 1 

. 00022 
OOIOO 

. 0002 1 
.36255 

. 0002 1 

. 00022 

6 1 575  . 

0039  1 

00198 

00  133 

OOIOO 

00082 

00072 

00064 

00057 

.00049 

00042 

. 00037 

00034 

. 00034 

.  00036 

.  00038 

. 00038 

.  00036 

. 00033 

. 00029 

. 00026 

. 00023 

. 00022 

. 0002 1 

.  0002 1 

. 00020 

. 00025 

. 0003 1 

.  00038 

. 00043 

. 00044 

. 00040 

. 00036 

. 00034 

. 00039 

. 0005 1 

OOIOO 

36252 

6  1565 

0039  1 

00198 

00133  . 

OOIOO 

00082 

00071 

00064 

00057 

00049 

00042 

00037 

00033 

. 00033 

. 00035 

.  00037 

. 00037 

00037 

. 00034 

.  00030 

00026 

00024 

00022 

0002  1 

. 00020 

. 00020 

. 0002 1 

00025 

. 00032 

. 00038 

. 00043 

. 00043 

. 00039 

00034 

. 00033 

. 00038 

. 0005 1 

OOIOO 

.  36250 

6  1555 

0039  1 

00198 

00  133 

OOIOO 

0008  1 

0007  1 

00064 

00057 

.00050 

00043 

00037 

00033 

. 0003 1 

.  00033 

00035 

. 00037 

. 00037 

. 00035 

. 0003 1 

00027 

. 00024 

. 00022 

. 0002 1 

00020 

. 00020 

. OOO 1 9 

.  00022 

. 00026 

. 00033 

00039 

00043 

. 00042 

. 00037 

. 00033 

. 00033 

00038 

. 0005 1 

OOIOO 

6  1546 

0039  1 

00198 

00  133 

OOIOO 

0008  1 

00070 

00063 

00057 

00050 

00043 

, 00037 

00033 

. 00030 

. 0003 1 

. 00033 

. 00035 

. 00036 

. 00035 

00032 

. 00029 

. 00025 

. 00022 

. 0002 1 

. 00020 

. 00020 

. OOO 1 9 
. 36245 

.  OOO  1  9 

00022 

00027 

00034 

00040 

0004  2 

00040 

. 00035 

0003  1 

. 00032 
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486  92 
00349 

00384 

0346  1 
00336 
. 00427 

,01816 

00325 

004  9  8 

0  1247 
003  16 

006  30 

009  5  9 
003  1  0 

009  73 

00785 

00305 
. 3 1789 

006  6  9 
00303 

00586 

00302 

00525 

00303 

.00477 

00306 

004  40 

003  1  2 

004  10 

00322 

00386 
0033  5 

0036  6 

0035  5 

35 

4859  1 

00345 

00349 

03456 
0033  1 

. 00378 

.01813 

00320 
004  2  1 

0  1245 

003  1  1 
0049  1 

00956 
. 00304 

006  23 

00783 

0029  9 
009  6  5 

006  6  6 
0029  5 

.31761 

00584 

0029  3 

00522 
. 00293 

004  7  4 

. 00295 

00437 

0029  9 

00407 

00305 

00382 
003  1  5 

0036  2 

00329 

36 

37 

38 

39 

40 

4  1 

42 

43 

44 

45 

46 

47 

4  8 

49 
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48492  .03451  01810  .01242  .00954  .00780  .00664  00581  .00519  .00472  .00434  00404  .00379  00358 

.00341  00327  .00316  .00306  .00299  00293  00288  00286  00285  .00285  .00287  00292  .00299  .00309 

.00323  .00343  00372  .00415  .00485  00617  .00957  .31734 


48396  .03447  01807  .01240  00952  .00778  .00662  .00579  .00517  00469  .00431  .00401  .00376  .00355 

.00338  .00324  .00312  .00302  .00294  00287  .00283  00279  00277  00277  .00278  00281  .00285  .00293 

.00303  .00318  .00338  00367  .00410  .00480  .00611  00950  31707 


48303 
. 00335 
. 00287 

0344  2 
. 00320 

.  0029  8 

. O 1 804 

. 00308 
.00312 

.01237 

.00298 

.00333 

009  50 
. 00290 
00362 

00776 

. 00283 
. 00405 

006  60 

. 00277 
00475 

. 00577 
.00273 

.00605 

. 005 1 5 
. 0027 1 

009  4  2 

0046  7 
. 00269 
.31682 

. 00429 

. 00269 

0039  8 
. 0027 1 

. 00373 
00274 

00352 
. 00279 

482  1  3 

. 00332 
. 00274 

03438 

.00317 

. 00282 

0  1802 
00305 

.00292 

0  1  235 
. 00294 
.  00307 

00947 

. 00286 

. 00328 

00774 
.00278 
. 00357 

.  006  5  8 

00273 
. 00400 

00575 

. 00268 
.  00470 

. 005 1 3 
0026  5 

. 00599 

0046  5 
. 00263 

009  35 

004  27 
. 00262 
.31657 

00396 

0026  3 

00370 

0026  5 

.  00350 

.  0026  8 

48  125 

. 00329 
. 00263 

03434 
003  1  4 

.  0026  9 

.01799 
. 00302 
.00277 

0  1  233 
. 0029 1 

. 00288 

00945 

00282 
. 00303 

00772 
. 00275 
. 00323 

006  56 

.  0026  8 

. 00352 

00573 
0026  4 

.  0039  5 

005  1  1 

. 00260 
. 00465 

00463 
00257 
. 00594 

. 00424 
. 00256 

. 00929 

0039  3 
.  00255 
.31633 

0036  8 
. 00256 

.  00347 
. 00259 

48039 

00327 

. 00253 

03429 

. 003  1  2 
. 00258 

0  1796 

. 00299 
. 00264 

,01231 
. 00288 

00272 

00944 

. 00279 

. 00283 

00770 
0027  1 

. 00298 

006  54 

. 00265 
. 003 1 9 

005  7  1 

00259 

. 00348 

. 00509 

00255 
0039  1 

.  0046  1 
00252 
00460 

. 00422 
.00250 

.  00589 

. 0039 1 
00249 

009  22 

0036  6 
0024  9 

.31610 

00345 
0025  1 

47955 

00324 

. 00245 

03425 
. 00309 

.  002  4  8 

0  1794 

.  0029  6 
. 00253 

.01229 
. 00285 

. 00259 

009  4  2 
00276 
. 00267 

0076  8 
. 00268 
. 00279 

006  52 
.00261 
0029  4 

0056  9 
00255 

. 003 1 5 

.  00507 

0025  1 
. 00344 

00459 
. 00248 

.  00387 

00420 

00245 

. 00456 

00389 

00244 

00584 

0036  4 
. 00243 

. 009 1 6 

.  00342 
.  00244 

.31588 

47874 

.  00322 
00238 
.31566 

03422 
. 00307 

. 00240 

0  1792 
0029  4 
0024  3 

0  1227 

. 00283 
. 00248 

009  40 

00273 
. 00255 

00767 

.00265 
0026  3 

006  50 

0025  8 
. 00275 

00567 

00252 
.  00290 

00505 

00247 

. 003 1 1 

004  5  7 
. 00243 

. 00340 

004  1  9 
00240 

00383 

0038  7 
00239 

. 00452 

0036  2 
00237 
. 00579 

.  00340 
.  00237 

.00910 

47794 
. 00320 

. 00232 
.00905 

034  1  8 
. 00305 

.00233 

.31545 

0  1790 
0029  1 

.  00235 

0  1225 

00280 
. 00239 

009  38 

.00270 
. 00244 

0076  5 

.00262 
. 00250 

006  4  9 

0025  5 

.  0025  9 

0056  6 
002  4  9 

.00271 

00504 
. 00244 

. 00286 

00455 
. 00240 

00307 

004  1  7 

.00236 
. 00336 

00386 

00234 
. 00379 

00360 

. 00232 
. 00448 

0033  8 
00232 
. 00575 

477  1  7 

.  003  1  8 
. 00227 

.  0057  1 

034  1  4 
. 00303 
. 00227 

. 00899 

0  1787 
00289 

00228 

.31525 

,01223 
. 00278 

. 0023 1 

00937 

.00268 

00235 

00763 

0  0  2  5  9 
00240 

006  47 

. 00252 
00247 

00564 

00246 

00255 

00502 
0024  1 

0026  7 

004  5  4 
.00236 

00283 

004  15 

. 00233 
00303 

00384 

00230 

00333 

00358 

. 00228 
. 00375 

00337 
.00227 
004  4  4 

4764  1 

.00316 
00222 
.  004  40 

034  1  1 

.00301 
00222 
. 00566 

0  1785 

00287 

00222 

.00894 

0  1222 
. 00276 
00224 

.31505 

00935 

.00266 

00227 

00762 

00257 
0023  1 

006  4  6 

.00250 

00236 

00563 
00243 
. 00243 

0050  1 

00238 

00252 

00452 

. 00233 
. 00264 

004  14 

.  00229 
00279 

00382 
.00226 
.  00300 

00357 

00224 

00329 

00335 
.00222 
.  00372 

47567 

003  1  5 
.00217 
.  0036  8 

03407 

0029  9 
.  002  1  7 
. 00436 

0  1783 

00285 
.00217 
. 00562 

0  1  220 
. 00274 
002  1  8 
. 00889 

009  34 
.00264 

.  00220 
.31486 

0076  1 

00255 
. 00223 

006  4  5 

00247 

00227 

005  6  2 
. 0024 1 

.  00232 

004  9  9 

00235 

00239 

00451 

. 00230 
00248 

004  1  2 

00226 

00260 

00381 

00223 

00276 

00355 

00220 
0029  7 

00333 

.00218 

00326 

47495 
. 003 1 3 
002  13 

00323 

03404 

. 00297 
002  1  2 
.  0036  5 

0  17  8  1 
.00284 
.00212 
00433 

0  12  18 
00272 
002  1  3 
OOS5  9 

009  32 
. 00262 
. 002 1 4 
. 00884 

00759 

00253 
002  1  6 
.31467 

006  43 
00245 
.00219 

005  60 
.00238 
00223 

004  9  8 

00232 

00229 

00450 

00227 
. 00236 

004  1  1 

. 00223 
00245 

00379 

00220 

00257 

00353 
002  1  7 
00273 

00332 
002  1  5 
.  0029  4 

47424 

.00312 
002  10 
.  0029  1 

0340  1 
.00296 

. 00208 
.00320 

0  1779 
.  00282 
00208 

. 00362 

0  12  17 
00270 
00208 
. 00430 

00931 

00260 

.00209 

. 00555 

00758 
0025  1 
002  10 
. 00879 

00642 
. 00243 
002  12 
.  3  1449 

005  5  9 
. 00236 
002  16 

00497 
00230 
. 00220 

004  4  8 
00225 

. 00226 

004  10 

0022  1 
. 00233 

00378 
002  1  7 

. 00242 

00352 
002  14 

00254 

00330 
002  11 
00270 

47354 

.00310 
00206 
0026  7 

0339  7 
. 00294 

.00205 
. 00288 

0  1777 
. 00280 

00204 
003  17 

0  12  15 
. 00268 

. 00203 
. 00359 

00930 
. 00258 

.00204 
004  26 

00757 

00249 

00205 

0055  1 

006  4  1 

0024  1 

. 00206 
008  7  4 

OOS58 

00234 

00209 

.31431 

004  9  5 

00228 

.00212 

00447 

00223 
002  17 

00408 
002  18 

00222 

00377 
002  14 

.  00230 

00351 

002  1  1 
00239 

00329 

00208 

.00251 
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K  =  .  1  5000 

N  =  50 


n 

1  1 . OOOOO 


2 

.60119 

3988  1 

3 

56384 

07470 

.36146 

4 

.54224 

.06479 

. 04726 

.34571 

5 

. 52707 

.06067 

.03858 

.03730 

.  336  38 

6 

.51541 

05830 

. 03507 

.  02923 

.  03200 

329  9  8 

7 

50596 

056  7  1 

.  033  1  1 

. 02599 

.  024  36 

0286  6 

.  3252  1 

8 

. 49803 

.05555 

03  184 

024  1  9 

.02128 

02  133 

.  02633 

.32145 

9 

.49119 

. 05465 

0309  4 

. 02303 

.01959 

.01839 

.01925 

024  59 

.31837 

1  O 

48520 

.05391 

03026 

. 02222 

.01850 

0  1676 

.01640 

. 0 1772 

.02324 

.31579 

1  1 

4  79  87 

.  05330 

. 02972 

.02161 

.  0 1 773 

.01572 

.01484 

0  1496 

.01654 

.02214 

.  3  1358 

1  2 

. 47508 

.  05277 

02928 

.02113 

.01716 

.01498 

.01383 

.01344 

0  1385 

0  1559 

.02124 

3  1  16  4 

1  3 

47072 

.05231 

.02891 

. 02075 

.01672 

.01444 

.01312 

.01246 

O  1  237 

.01297 

0  1482 

02047 

. 30994 

1  4 

. 46672 

.05191 

.02860 

. 02043 

.01636 

.01402 

.01260 

.01177 

.01  140 

0  1150 

O  1225 

.01417 

.01982 

. 30841 

1  5 

46304 
. 30704 

.05154 

.02832 

020  1  6 

0  1607 

.01368 

.01219 

.01130 

. 01070 

.01060 

.01080 

.01170 

0  1362 

.01925 

1  6 

45963 

.01874 

.05121 
. 30579 

02808 

.01993 

0  1582 

0  1340 

.01186 

.01090 

0  1030 

009  9  4 

. 009  9  2 

0  1030 

oiiio 

.01315 

1  7 

456  45 
.01273 

0509  1 

.01830 

.02787 

.  3046  5 

.01972 

.01560 

.01316 

.01 160 

0  1060 

009  87 

. 00945 

.  009  28 

.00937 

. 00978 

0  1070 

1  8 

.  45347 

.01030 

0506  3 
.01236 

.02767 
. 01790 

.01954 

. 30359 

0  1542 

0  1296 

.01  140 

.01030 

. 00956 

.  00908 

0088  1 

.  008  7  4 

00890 

. 00937 

1  9 

. 45067 

. 00900 

05038 

.  009  9  9 

.  02750 
. 0 1 203 

0  19  38 
.01753 

0  1526 

.  3026  2 

0  1279 

.01 120 

.  0  1  O  1  0 

. 00930 

00878 

00844 

00828 

. 00829 

.  00850 

20 

. 44803 
.  008  1  5 

050  14 
.00869 

02734 
. 00970 

0  19  24 
.  O 1  174 

.01511 

.01720 

0  1264 
.30171 

. O 1 1 00 

009  89 

00909 

. 00853 

008  15 

. 00792 

00783 

00790 

2  1 

. 44554 
. 00756 

049  9  1 
00785 

027  1  9 
. 0084 1 

.01911 
. 00943 

0  1499 

0  1150 

0  1250 

0  16  90 

.01090 
. 30087 

.  009  73 

0089  1 

00832 

00790 

.00763 

. 00747 

. 00744 

22 

.44317 

.007  1  1 

049  70 
. 00726 

02706 
.  00758 

.01899 
. 008 1 6 

0  14  87 
009  1  9 

0  1238 

0  1120 

. 0 1 070 

.01662 

009  59 
30008 

00876 

. 008 1 5 

00770 

.00739 

. 007 1 9 

. 00709 

23 

44092 

. 00677 

049  50 

006  82 

026  9  3 
00700 

0  1888 
. 00734 

.01476 

. 00793 

0  1228 
0089  7 

0  1060 

O  1  1  OO 

.  009  47 

0  16  37 

0086  2 

.29933 

00800 

00753 

. 0071 9 

006  9  5 

006  8  1 

24 

43877 

. 00649 

049  32 
006  4  8 

026  8  1 
006  5  7 

0  1878 
006  77 

0  1467 
007  1  2 

0  12  18 
00773 

. O 1 050 
. 00877 

009  36 
. 01080 

00850 

.01613 

007  8  6 

.  29  86  3 

00738 

00702 

00676 

.  006  5  8 

25 

43673 
. 00627 

049  1  4 

.  006  2 1 

.02670 
006  23 

0  1869 

006  34 

0  1458 

006  5  6 

0  1  209 
006  9  3 

0  1040 

00754 

009  26 
0085  9 

00840 

.01060 

.  00775 

0  15  9  1 

00725 

29  79  7 

006  8  7 

. 00659 

.  006  39 

26 

43477 

00608 

0489  7 
. OOS99 

026  S  0 
00596 

0  1860 

. 0060 1 

0  1450 

006  1  3 

0  1201 
. 00637 

. 0 1030 

006  75 

009  1  7 

. 00737 

00830 

. 00842 

.  0076  4 

0  1040 

007  1  4 

0  1570 

006  75 

. 29734 

006  45 

006  23 

27 

43289 
. 0059 1 

04  88  1 

005  80 

026  50 
00574 

0  1852 
0057  4 

0  1442 

0058  1 

0  119  3 

00595 

0  1  030 
006  1  9 

00909 

. 00659 

0082  1 

007  2  1 

00755 

00826 

00704 

0  1030 

006  6  4 

.01551 

006  33 
.  29  6  74 

. 00609 

28 

4  3  109 
. 00578 

0486  6 

. 00564 

026  40 

OOS56 

.01844 

00552 

0  1435 
00554 

0  1186 

0056  2 

.01020 

005  7  8 

. 0090 1 

00604 

008  14 

00644 

.  00747 

00707 

.  006  9  5 
.00812 

006  54 
.01010 

006  2  2 

0  1533 

.  005  9  7 
296  1  8 

29 

4  29  36 
005  6  6 

.29564 

0485  1 

. 00550 

026  32 
. 00540 

.01837 

00534 

0  14  28 

00533 

0  1180 

005  36 

0  10  10 

00546 

.  008  9  5 
005  6  3 

00807 

. 00589 

00739 

006  30 

006  8  6 

006  9  3 

006  4  5 
0079  8 

006  1  2 
009  9  7 

00586 

0  15  16 

30 

42769 

. 00555 
0  1499 

.04837 

00538 

.29512 

026  23 

005  26 

.01830 
. 005 1 9 

0  1422 
005  1  5 

0  1  17  4 

005  1  5 

.  O  1  0  1  0 

00520 

008  88 

0053  1 

00800 

0054  9 

.  00732 
OOS76 

.00679 

006  1  7 

006  37 

006  8  1 

00603 

00786 

00576 

009  8  4 

3  1 

4  2608 
0054  5 

. 0097 1 

04823 

005  2  8 

0  1484 

026  1  5 
.  005  1  5 

29462 

0  1824 
00505 

0  14  16 
004  9  9 

0  1  170 

004  9  7 

. 0 1 OOO 

0049  9 

008  83 

. 00506 

0079  4 

005  1  7 

00726 

00536 

006  72 

005  6  4 

006  30 
00605 

005  9  5 
006  6  9 

005  6  8 
0077  4 

32 

42453 
00537 
. 00763 

048  10 

005  1  9 

.  009  5  9 

02607 

00504 

.01470 

.01818 

004  9  4 

.29415 

0  14  11 

00486 

O  1  160 
.  00482 

009  9  6 

.  004  8  2 

008  77 

00485 

007  8  8 

004  9  2 

007  20 
00505 

006  6  6 

005  24 

006  23 
005  5  2 

005  8  8 
0059  4 

005  60 

006  5  8 

33 

42303 

00529 

006  4  8 

04798 

005  1  0 
00752 

02600 
004  9  5 
009  48 

0  18  12 
0048  3 

0  1456 

0  1  405 

00475 

29  36  9 

0  1  160 

0046  9 

.00991 

0046  7 

008  72 
00467 

00783 
0047  1 

007  1  5 

004  80 

00660 

00493 

006  1  7 

005  1  3 

005  82 
00542 

0055  3 
. 00584 

34 

42  158 

. 00522 
005  7  4 

04785 

. 00503 
. 00638 

025  9  3 
00487 
00742 

0  1806 
00474 

00937 

0  1400 

004  6  5 
.01443 

0  1150 
00458 

.29326 

009  87 
004  54 

008  6  8 

. 00452 

00778 
004  5  4 

007  10 

004  5  9 

00655 
0046  8 

006  1  1 
0048  2 

00576 

. 00502 

0054  6 

00532 

35 

4  2018 
005  1  6 

00522 

04  774 
.  004  9  6 

.  005  6  5 

02586 
004  79 

006  29 

0  1801 

0046  6 
00733 

0  1396 
0045  6 

009  27 

0  1150 
004  4  8 

0  1430 

00982 

00442 

29  284 

0086  3 

00440 

00774 

00439 

00705 

00442 

006  50 

004  4  8 

00606 

004  5  8 

00570 
004  72 

0054  1 

0049  3 

' 


36 

37 

36 

39 

40 

4  1 

42 

43 

44 

45 

46 

4  7 

48 

4  9 

SO 
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4  1882 
005  1  O 

. 00484 

04762 

. 00490 

. 005 1 4 

.  02580 
.  00473 
.  00556 

.01796 
00459 
. 00620 

.01391 
00447 
. 00724 

.01140 
. 00439 

.  009  1  8 

.  009  7  8 

004  32 
.01418 

00859 

. 00428 

29  243 

00770 

. 00427 

0070  1 

.  004  27 

006  46 
. 0043 1 

0060  1 

00437 

0056  5 

.  00448 

. 00535 
0046  3 

4  1750 

00505 

.  00454 

.04751 

0048  4 

. 00475 

02574 
0046  7 

00505 

0  1791 
. 00452 
00548 

.01387 

00440 

.  006  1  2 

0  1140 

. 0043 1 

.  007 1 6 

00974 

00424 

.  00909 

00855 
004  1  9 

.01407 

0076  6 
. 004 1 6 

. 29204 

.  006  9  7 
. 004 1 5 

00642 

. 004 1 6 

. 00597 

004  2  1 

.  0056  1 

. 00428 

00530 

. 00439 

4  1621 

00500 
. 00430 

0474  1 
00479 

. 00446 

02567 
.00461 
. 00468 

0  1786 

. 00446 
. 00498 

.01383 
. 00434 

. 0054 1 

.01  140 
.00424 

. 00605 

. 0097 1 

004  1  6 

.  00708 

00852 

004  1  0 

.  00900 

0076  2 

. 00406 

.01396 

.  006  9  3 
. 00404 

.29166 

006  38 
00404 
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.01491 
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K  =  .  5  5000 

N  =  50 

n 

1  1  OOOOO 
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.02462 

.  026  54 

.02976 

.  03648 

0660  1 

05843 

04276 

.  036  74 

03320 

03075 

.02891 

.02747 
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.01948 

0  1930 

.01916 

.01905 

.01898 

.01893 

.01893 

.01896 

.01904 

.01918 

.01937 

.01964 

. 02000 

. 02048 

.02112 

. 02200 

. 02325 

025  1  2 

.02824 
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.01955 

.01931 

0  1909 
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.01855 

0286  5 

.01894 

.01883 

.02695 

.01873 

.01920 

02560 
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026  5  1 
.01838 
0  1846 

025  1  8 

.01819 

.01883 

024  10 

0  1802 

.01932 

02320 

.01789 

.01996 

02243 
. O 1 777 

.02083 

.02176 

.01769 
. 02206 

.02118 
.01763 
.  02388 

0206  6 

0  1760 

.  026  9  1 

0202  1 

.01760 

033  19 

05246 

.01948 

.01720 

.03272 

03846 

0  19  12 

.01724 
. 06000 

03309 

0  1881 

.01732 

029  9  2 

.01853 

.01744 

02773 

.01827 

.01760 

.  02608 

.01805 

.01782 

024  78 
.01786 

.01811 

02372 

.01769 

0  1848 

.  02283 

.01754 

.01897 

02207 

.01742 

.01961 

02141 

.01733 

. 02048 

02083 

.01726 

.02169 

02033 

.01721 

02350 

0  19  88 

0  17  19 

.  026  50 

05  160 

0  19  16 

.01681 

026  1  1 

03784 

.01881 

0  1683 
. 03226 

03256 

.01850 

.01688 

.05926 

02945 

.01821 

0  16  9  6 

.02729 

.01796 

.01708 

0256  7 

0  1774 
.01725 

02439 

.01754 

.01748 

. 02335 
.01737 
.01777 

02247 

. 0 1 722 

.01814 

02  172 
. 0 1 709 

.01863 

02  107 

0  16  9  9 

0  19  27 

. 02050 

.01691 
.  020 1 4 

02000 

0  1685 
02  134 

0  19  5  6 

0  16  8  1 

.02314 

05077 

.01886 

.01644 

. 02279 

03724 

.01851 

0  16  44 

. 02573 

03205 

.01820 
0  16  47 

03  183 

02899 

.01792 

.01653 

. 05856 

026  87 

.01766 

.01662 

02528 

.01744 

.01675 

02402 

.01724 

.01692 

.02299 
.  O 1 706 
.01715 

022  1  2 

.01691 

.01745 

02  139 

0  16  7  8 
0  1782 

02075 

0  16  6  7 
.01831 

020  1  9 
.01658 

.01895 

0  19  6  9 
.01651 

.01981 

0  19  25 

.01646 

.02101 

04998 

.01857 

.01610 

. 02069 

03667 

0  1822 

.01609 

02245 

03  156 

.01791 

.01610 

.02537 

02855 

.01763 
0  16  13 
03141 

02647 

.01738 

.01620 

05787 

02490 

0  17  15 
.01630 

0236  6 
.01695 
0  16  43 

0226  5 

0  1677 
.01661 

02  179 

0  16  6  2 
.01684 

02  107 

.01648 

.01714 

0204  3 

0  16  36 

0  1752 

0  19  88 
.01627 

0  1801 

0  19  39 

.01619 

0  1864 

0  1896 

0  16  13 

0  19  50 

0492  1 

0  1829 

0  15  7  8 
.01920 

036  1  1 

0  1795 
.01575 

.02038 

03  109 

0  1764 

.01575 

022  1  3 

028  13 

0  1736 
.01577 

. 02503 

026  08 
.01711 

.01581 
03  100 

02453 

0  16  88 
0  1588 
.05721 

0233  1 
.01668 
.01599 

0223  1 

0  16  50 
0  16  13 

02  147 

.01634 

.01631 

02076 

0  16  20 

.01654 

020  13 

0  1608 

0  16  8  4 

0  19  5  9 
.01597 

. O 1 722 

0  19  11 

.01589 

.01771 

0  1868 

0  1582 
0  1835 

04847 

0  1802 

0  15  4  8 

.01807 

03558 

.01768 

0  1544 

0  18  9  1 

0306  4 
.01737 

.01542 
. 02009 

02772  . 

0  17  10 
0  1543 
.02182 

02570 

0  16  85 

0  1545 
02  46  9 

024  1  8 

0  16  6  2 
.01550 
.  0306  1 

0229  8 
.01642 

.01558 
056  58 

02199 

.01623 

.01569 

02117 

.01607 

0  1584 

0204  6 

.01593 

.01602 

0  19  84 

0  1580 

0  16  26 

0  19  3  1 

0  1569 

0  16  5  6 

0  1883 

0  1560 

0  16  9  4 

0  1840 
.01553 
0  1743 

04775 

0  1777 

0  15  19 

.01716 

03506 

.01743 

.01515 

0  1779 

03020 

0  17  12 

.01512 
0  1863 

02732 

.01684 

.0151 1 
. 01980 

02533 

.01660 

0  15  12 

.02152 

02384 

0  16  37 

.01516 

.02437 

02265 

0  16  17 

0  1521 

. 03023 

02168 

0  15  9  8 
.01530 

.05596 

02087 

0  1582 
.01541 

020  17 

0  1567 

.01556 

0  19  5  6 
.01554 

0  1575 

0  1903 
.01543 

0  15  9  9 

0  1856 

0  1533 

0  16  29 

0  18  14 

.01526 

.01667 
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K  =  .  80000 

N  =  50 

n 

1  1 . ooooo 


2 

. 5 1 505 

4849  5 

3 

.36741 

.29529 

. 33731 

4 

.29144 

. 22789 

. 2 1 504 

.26563 

5 

.24412 

. 18929 

.  17  184 

. 17216 

.  22259 

6 

.21142 

. 1 6348 

.  1  4628 

. 14012 

.  14514 

.  1  9356 

7 

.  1  8730 

. 14472 

.  1  2863 

. 12100 

. 11940 

.  1  264  1 

.  1  7249 

8 

.  1  6  868 

. 1 3036 

.  11550 

. 10770 

. 10410 

. 10470 

.  1  1260 

.  1  564  1 

9 

. 15381 

. 11890 

. 10510 

09760 

.09340 

.09189 

09362 

. 10190 

1  4367 

1  0 

.14163 

.  109  60 

.09680 

.  08958 

.08522 

.08290 

08257 

. 08499 

09  34  1 

. 13329 

1  1 

. 13145 

.  10  180 

.08987 

0830  1 

.07868 

. 07603 

.07482 

. 07520 

. 07803 

.08645 

.  1  2460 

1  2 

. 12280 

.09519 

.08401 

07750 

.  07328 

0705  1 

.06888 

06  837 

06922 

.07230 

. 08064 

.  11730 

1  3 

.11530 

08949 

.07897 

.  07280 

.06871 

06  59  2 

.06410 

063  1  4 

06  309 

.06425 

06  748 

.  07570 

. 1 1 1 OO 

1  4 

. 10880 

.  0845  1 

.07459 

.06872 

.06479 

06  203 

060  1  3 

. 05892 

0584  1 

. 05868 

. 06006 

06337 

.07144 

105  50 

1  5 

.  10310 
10070 

.08013 

. 07073 

.06515 

06  137 

.05867 

.05674 

.05541 

054  6  3 

. 05443 

.  05493 

056  46 

.05962 

.06773 

1  6 

.  09803 
.06447 

.07624 

.09637 

.06731 

06198 

.05835 

.05573 

.05380 

.05241 

.05148 

05101 

.05104 

.05171 

.05334 

.05672 

1  7 

. 09349 
.05398 

07276 

.06157 

06  425 
09253 

.05916 

.05567 

.05312 

05  122 

.04980 

04879 

048  1  6 

.04791 

048  1  1 

. 04890 

05060 

1  8 

.08939 

.04818 

.06962 

.05154 

06  150 
05898 

056  6  2 

. 08906 

.05326 

. 05079 

04892 

04750 

.  04645 

. 04571 

.04530 

04522 

.04555 

.04643 

1  9 

.08568 
. 04424 

.06677 

04602 

05  900 

.04935 

.05432 

05664 

.05109 

.08591 

.04870 

046  87 

.04546 

04438 

04358 

. 04306 

.04280 

.04286 

04329 

20 

.08230 

04  127 

06418 

04228 

.05672 
. 04407 

05223 

04738 

049  1  1 

.05452 

046  80 

. 08304 

.  04  50  1 

.04362 

04  253 

.04170 

.04111 

.  04074 

. 0406 1 

.04077 

2  1 

079  2  1 

03890 

06  180 

.03947 

05464 
. 0405 1 

.05031 

. 04232 

.04730 

.04559 

04506 

.05259 

04333 

.08041 

04196 

.04087 

. 04002 

.03938 

03894 

038  70 

03867 

22 

.07637 

.03693 

.05962 

03723 

05273 

03784 

04855 
.  038?  1 

04565 

. 04073 

.04348 

.04396 

04  179 

.  05083 

04044 

0779  9 

039  36 

. 03851 

. 03784 

.03735 

03702 

036  88 

23 

07376 
. 03525 

0576  1 
.03537 

05096 
. 03572 

046  9  3 
.03637 

.04412 

. 03746 

.04201 

.03927 

04037 

.04247 

03905 
. 04920 

03799 

.07575 

037  1  3 

036  45 

.  03592 

03554 

.  0353  1 

24 

07  133 
03378 

05  57  4 

.03378 

04  9  32 
0339  5 

.04543 

03434 

0427  1 

03502 

0406  6 

036  1  3 

03906 

.03793 

.03777 

04  109 

03673 

04770 

.03587 
0736  8 

.  035  1  8 

.03463 

.03421 

0339  3 

25 

.06908 

03248 

.05401 

03239 

04780 

.03245 

04403 

.03266 

. 04140 

. 03309 

.03941 

03379 

03785 

.03491 

.03659 
0367  1 

03556 

. 03982 

.  03472 
04632 

.  03  402 
07  175 

.03346 

.03301 

0326  9 

26 

06699 

.03132 

05239 

.03117 

04  6  38 
.03114 

04273 

.03123 

040  1  8 

.03148 

03825 

.03193 

.  03673 
.03266 

.03550 

03378 

.03449 

03557 

0336  5 
0386  4 

0329  5 
04502 

03238 

06995 

03192 

.03157 

27 

06  503 

. 03026 

05089 

.  03007 

04  506 
.02998 

04  152 

. 02999 

03904 

.  030 1  2 

037  16 
03040 

0356  8 

.  03087 

03448 

.03161 

03348 
. 03274 

0326  6 
0345  1 

03  197 
0375  5 

03  139 

04382 

0309  2 

06826 

0305  5 

28 

06  320 
029  30 

04947 

.02908 

04382 
0289  4 

04038 

02889 

.03797 

.02893 

036  1  4 
029  10 

03470 

02940 

03352 
. 02989 

.03255 

. 03064 

03  174 

03  177 

03  105 
03353 

.  03048 

036  53 

. 03000 
04  270 

029  6  1 

06668 

29 

06148 

02842 
065  1  9 

.04814 

.02817 

04265 
.  02800 

.03931 

.02790 

036  9  7 
.02788 

035  19 
.02796 

03378 

.02815 

03263 

02848 

.03168 

02899 

. 03087 
02974 

. 03020 
03087 

029  6  2 
03262 

.02914 

0  3  5  5  8 

02874 

.04164 

30 

.05986 

.02761 

.04065 

046  90 

.02734 
.  06379 

04  156 

. 027 1 4 

0383  1 
.02701 

03602 

.02694 

03429 

.02696 

.  0329  2 
02707 

03  179 

02728 

03086 
0276  2 

03007 

.02814 

029  40 
02890 

02883 

03003 

. 02835 
.03176 

.02794 

03468 

3  1 

05834 
026  85 

03384 

.04572 

026  57 
039  72 

04052 
026  35 

.06247 

03736 
026  1  9 

035  1  3 
02609 

03344 

02606 

032  10 

026  1  1 

03  100 

026  23 

. 03009 
02647 

029  3  1 
026  82 

0286  5 

02735 

02809 

028  1  2 

02760 
029  24 

027  19 
0309  6 

32 

056  90 
026  1  5 
0302  1 

0446  1 

02586 

03305 

03954 

.02562 

03884 

036  46 

02544 

.06122 

.03429 

.02532 

0326  4 

02525 

03  133 
02525 

03026 

02531 

.02936 

02546 

02860 
0257  1 

0279  5 

02608 

02739 
0266  1 

026  9  1 
02738 

.02650 

02850 

33 

05554 

.02549 

.02781 

04355 

.02519 

029  49 

0386  2 

02495 

03230 

.  0356  1 
.02475 

. 03801 

03350 

.02461 

06003 

03189 

.02451 

0306  1 
02447 

029  5  6 
.  02449 

0286  8 
02457 

0279  3 
. 02474 

.02729 

02  500 

026  73 
02538 

026  26 
02592 

025  8  4 
.02669 

34 

05425 

. 02488 
02604 

04256 

02457 

027  15 

03774 

.02432 

.02882 

0348  1 

.024  1  1 
03  160 

03274 

0239  5 
03722 

03117 

.02383 

.05890 

029  9  2 
.02376 

02889 

02374 

02  803 

02378 

02729 

02389 

026  66 
02406 

026  1  1 

02434 

0256  4 

02473 

025  23 
02527 

35 

05302 

02430 

02466 

04161 

02399 

02543 

036  9  1 
02373 

026  53 

03405 

.02351 

028  1  9 

03203 

02334 
0309  3 

0304  9 

02320 

03648 

029  26 
023  1  1 
05783 

02826 

02306 

0274  1 
02307 

026  6  9 

023  12 

02607 
0232  4 

02553 

02343 

02506 
0237  1 

024  6  5 
024  11 

36 

37 

38 

39 

40 

4  1 

42 

43 

44 

45 

4  6 

47 

4  8 

4  9 

50 
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05  186 

.02375 

.02353 

. 0407 1 

02344 
.  02408 

.  036  1  2 

.  023  1  8 
.  02485 

, 03332 
02295 

. 02595 

.03135 
. 02276 

.  02759 

.  02984 

02262 

.03030 

. 02864 

.02251 

.03577 

.  02766 

. 02244 

.05681 

.02683 
0224  1 

.02612 

02243 

0255  1 

02250 

.02497 

.02263 

.02451 

.  02284 

.02410 

023  13 

05075 

.  02323 

.02257 

. 03985 
.02292 
.02298 

.  03536 

.02266 

.02354 

03263 

.02242 

.02431 

.  03070 
.  02223 
.02539 

02923 

. 02207 

.  02702 

. 02805 

.02195 

.02970 

. 02709 
.02186 

.  03509 

.02627 

.02181 

.05583 

.  02557 

02  180 

02497 

.02184 

02445 

.02192 

0239  9 

.  02207 

02359 

.02228 

04  9  70 

.02274 

.02175 

03903 
.02243 
.  02205 

. 03464 

.  022  1  6 

.  02246 

03  197 
.02193 
. 02302 

03008 
.02173 
.  02379 

0286  4 

.02156 

.02487 

02749 

.02143 

.02648 

.  026  54 

.02133 
029  1  2 

02574 

02  126 
.03445 

.02506 

.02123 

.05490 

024  4  7 

02  123 

02395 

.02128 

.02349 

02  138 

02309 
02  153 

04  86  9 

.02228 

.02102 

.03825 

.02197 

02  125 

03396 

.02170 
02  156 

03  134 

02  146 

.02197 

,  02949 
.02126 

.  02253 

02808 

.02108 

. 02329 

02695 
0209  4 

02437 

.  02602 

. 02083 

.02596 

02524 
. 02074 
.02858 

02456 

. 02069 
.03384 

.  0239  8 
.  0206  7 

.05401 

02347 
. 02069 

02302 

02075 

0226  3 
02086 

04773 
.02184 
. 02037 

.03751 
.02153 
.  02054 

.  03330 
.02126 

.  02078 

.  03073 
.02102 
.02109 

02892 

.  0208  1 
.02151 

02754 

02063 

. 02207 

026  44 
. 02048 
02283 

.  02552 
. 02036 

. 02389 

02475 
. 02026 
.02547 

02409 

. 02020 

02806 

02352 
.  020 1 6 
03325 

02302 
. 0201 6 

.05316 

.02258 
020  1  9 

022  1  8 

.  02026 

046  8  1 

02  142 
.01979 

03679 

.02111 

.01991 

03267 

. 02084 
. 02009 

030  1  6 

02060 

. 02033 

.02838 
. 02039 

.02065 

02703 

. 02020 
.02107 

02595 

.  02005 

.02162 

02505 

.01992 

. 02238 

02429 

0  19  8  1 
.02344 

0236  4 

.01973 

.02500 

02308 

.01968 

.02756 

02258 

.01966 

0326  9 

022  1  5 

.01967 

.05234 

02  176 

0  19  7  1 

0459  3 
.02101 
.01926 

.03611 

. 02071 

.01935 

.  03207 
.02044 

0  19  4  8 

029  6  1 
.  02020 

.01966 

.  02787 

.01999 

.01990 

.02654 

0  1980 

.  02022 

.  02547 

0  1964 

.02064 

.02459 

0  19  50 

.02120 

.  02385 
.01939 

.02195 

0232  1 
.01930 

.  02300 

02266 
.01924 
. 02456 

022  17 

.01920 
. 02709 

.02174 
.01919 
.  032  1  5 

02  136 
.01921 

.05156 

04508 

.  0206  3 
.01877 
05080 

03545 
. 02033 
.01883 

.03149 
. 02006 

.01392 

02908 

. 01982 
.01906 

.02737 

.01961 

.01925 

02607 

.01942 

0  19  50 

02502 

.01925 

.01982 

024  1  6 

.01911 

. 02024 

.02343 

.01899 

02080 

.02280 

0  1889 

.02155 

02225 

.01882 

.02259 

02  177 
.01877 

.02413 

.02135 

0  18  7  4 

026  6  3 

.02097 

0  1874 

.03164 

04427 
. 02026 

.01832 

.03114 

03482 

0  1997 

.01836 
. 05008 

.03094 

.01970 

0  1842 

0285  7 
.01946 

.01852 

.02689 

.01924 

.01867 

0256  1 

0  1905 
.01886 

02459 

.01888 

.0191  1 

02374 

.01874 

.01944 

02302 

.01861 

.01986 

02240 

.01851 

.02042 

02  187 

01843 
02  1  16 

02  139 

.01836 

.02220 

02097 

.01833 

.02372 

,  02060 

0  1831 
026  1  9 

04350 

.01991 

. 01790 
.02578 

03422 

.01962 
.01792 
.  0306  7 

03040 

.01935 

0  1796 

.04938 

02808 
.01911 
. O 1 804 

.02643 

.01890 

0  18  14 

025  18 

.01870 

.01829 

024  17 

.01853 

.01849 

02334 

0  18  38 
.01874 

02263 

0  1826 
.  0 1  907 

02202 

.01815 

.01949 

02149 
. 0 1 806 

.02005 

02  103 
.01799 
. 02079 

0206  1 

.01794 

.02182 

0202  4 

.01791 

.02332 

04275 

0  1958 
0  1751 

0229  5 

0336  4 
.01928 

.01751 

.02537 

02989 

.01902 
.01754 
. 03021 

0276  1 
.01878 

0  1759 

.04871 

02599 

.01856 

.01767 

02476 

.01837 

.01778 

02377 

.01820 

.01793 

0229  5 

.01805 

.01813 

02226 

0  1792 
.01839 

02166 

0  17  80 
.01872 

02114 

.01771 

.01914 

02068 

.01763 

.01970 

02027 
.01757 
. 02043 

0  19  9  0 

.01753 

.02145 

04203 

0  19  25 
.01714 

02  1  10 

03308 

.01896 

.01713 

.02258 

029  40 
.01870 

.01714 

.02499 

027  1  6 
.01846 

.01717 
029  77 

02557 

.01825 

0  1  723 
.04806 

02436 

0  1805 
. O 1 732 

02338 

.01788 

.01744 
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.03197 

03339 

03254 
.  03196 

03244 

03195 

03236 

.03195 

03230 

03196 

03224 

03196 

032  19 
.03198 

.03215 

03200 

032  1  1 
03204 

03208 

. 03209 

03205 

032  1  6 

03202 

.03226 

32 

032  19 
03101 

.03127 

03  184 
. 03099 

.03142 

03166 
0309  8 
03  170 

03  153 

.  0309  7 
. 03237 

0  3  14  4 

0309  6 

03  136 

0309  5 

03  130 
0309  5 

03  124 

0309  6 

03  120 

.  0309  7 

03116 

. 03099 

.03  1  12 

.03101 

03  109 

.03105 

03  1  06 

.03109 

03  103 

03116 

33 

.03122 
. 03008 
.  03023 

03088 

. 03006 
. 03033 

0307  1 
03005 
03048 

03059 

03003 
.  03075 

03050 

03002 

.03140 

0304  2 
03002 

03036 

03002 

0303  1 

.  03002 

03026 

. 03002 

0302  2 
.  03003 

.03019 

03005 

.03015 
. 03008 

030  1  3 
030  11 

030  10 

030  16 

3  A 

0303  1 

.  029  20 
02928 

029  9  9 
029  1  8 

. 02935 

029  82 
029  17 
029  45 

029  70 

.02916 
029  6  0 

029  6  1 

029  1  5 
029  86 

029  54 

029  1  4 
03050 

029  48 
029  1  3 

029  43 

029  1  3 

02938 

.02913 

029  3  4 
029  1  4 

.02931 

.02915 

02928 

029  1  7 

02925 

02  9  20 

029  22 
. 02923 

35 

. 02945 

.02838 

0284  1 

029  1  4 

.02836 

.02845 

02897 

02834 

02852 

02886 

.02833 

.02862 

02878 
.02832 
.  02  8  77 

0287  1 

02  83  1 
02902 

0286  5 

.02830 

029  6  5 

02860 

02830 

02855 

. 02830 

0285  1 
02830 

02848 

.02831 

02845 

02832 

02842 

02834 

02840 

02837 

■ 

36 

37 

38 

39 

40 

4  1 

42 

43 

4  4 

45 

46 

47 

48 

49 

50 
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0286  5 

. 02760 

. 02759 

028  34 

.02758 

.02762 

028  18 
.02756 
. 02767 

02807 
. 02755 

.02774 

. 02799 
. 02754 
. 02784 

.02792 
. 02753 
0279  8 

.02786 
. 02752 
.  02823 

.02781 

.02751 

.02884 

.  02777 
.  0275  1 

. 02773 

.  0275  1 

02770 

.02752 

. 02767 
. 02753 

. 02764 
.02754 

0276  2 

02756 

02788 
. 02686 
. 02682 

.  02758 
.02684 

.02685 

02742 

.02682 

.02688 

02732 
. 0268 1 
.02693 

. 02724 
. 02680 

.  02700 

027  1  7 
. 02679 
. 02709 

027  1  2 
.02678 
. 02724 

02707 
. 02677 
.02748 

. 02703 

.  02677 

. 02808 

026  9  9 
02677 

026  9  6 

.02677 

. 02693 

.02678 

.02690 

.02679 

02688 

.02680 

027  15 
.02616 
026  10 

026  86 

.02614 

.02612 

0267  1 
026  1  2 

.  026  1  5 

0266  1 
.  026  1  1 

.02618 

02653 

.02610 

.02623 

02647 

.02609 

.02630 

026  4  1 
.02608 

.02639 

.02637 

.02607 

02653 

.02633 

.02607 

.02677 

026  29 

. 02606 

02735 

026  26 
.02606 

. 02623 
.02607 

026  20 
.02607 

.02618 

. 02609 

02646 

. 02550 
. 02542 

026  1  8 
.02548 
.  02543 

02603 
02546 
. 02546 

.02593 

.02545 

.02548 

025  86 
.  02543 
02552 

02580 

. 02542 

. 02557 

02574 

025  4  1 

.  02563 

02570 
.02541 
. 02572 

.  02566 
.02540 
.02586 

.02562 

.02540 
026  10 

.  02559 
.02540 
.  02667 

02556 

.02540 

. 02554 

.02540 

02552 

025  4  1 

0258  1 

.02486 

02478 

025  5  3 

.02485 

.02479 

02539 

.02483 

.02480 

.  02529 
. 02482 
.  02482 

025  22 
.02480 
02485 

025  1  6 

.02479 

.02489 

025  1  1 

.02478 
.  02493 

02506 
.02477 
. 02500 

02502 

.02477 

.  02509 

. 02499 

02476 
.  02522 

02496 

.02476 

02546 

.  02493 
.02476 

. 02601 

0249  1 
.02476 

02489 
.  02477 

025  1  9 
.02427 

024  1  7 

02492 

.02425 

.02417 

02478 
.  02423 
.02419 

0246  8 
. 02422 
.02420 

.02461 

.02420 

.02422 

02455 

.02419 

.02425 

02450 

.02418 

. 02429 

02446 
. 024 1 7 
.  02433 

.  02442 
.02417 
. 02440 

02439 

.02416 

.  02449 

02436 

024  16 
02462 

02433 

.02416 

.02485 

0243  1 

.  024  1  6 
02539 

02429 

.02416 

02459 
0236  9 
.02359 

02433 
0236  8 

.  02359 

.02419 
0236  6 
.02360 

024  10 
02364 
.  0236  1 

02403 

.  02363 

.  02363 

02397 

.02362 
. 02365 

0239  3 
.02361 
. 02368 

02388 
.02360 
. 0237 1 

02385 

.02359 

. 02376 

0238  1 

02359 

.02382 

.02379 

02358 

.02391 

02376 

02358 

02404 

. 02374 
.02358 

.  02427 

0237  1 
02358 

02480 

02403 

023  1  5 

.  02303 
02424 

02377 

.02313 

.  02  304 

0236  4 

. 0231 2 
. 02305 

02355 
.02310 
. 02305 

.02348 
. 02309 
.  02307 

02342 

. 02308 
02308 

02338 

.  02307 
.02310 

02333 
. 02306 

.02313 

02330 

02305 

.  023  1  7 

02327 
. 02304 

.  0232  1 

02324 

02304 

.  02328 

0232  1 
.  02303 

.02336 

.02319 

02303 

.02349 

023  1  7 
. 02303 
0237  1 

02349 

. 02263 

0225  1 
.02318 

02324 

0226  1 
.02251 

.  02370 

023  11 
02260 

.02252 

02302 

02258 

02252 

02295 

.02257 

.  02253 

02290 

.02256 

. 02255 

02285 
.  02255 
. 02256 

0228  1 

.02254 

.02258 

02278 

.02253 
.  0226  1 

.02274 
02252 
0226  5 

02272 

02252 
0226  9 

0226  9 

. 022S 1 
.02275 

0226  7 

.02251 
.  02284 

0226  5 

0225  1 
0229  7 

02297 
. 022 1 3 
.02201 

.02246 

02273 
.022  1  1 

0220  1 
.02268 

02260 

.02210 

.02201 

.02318 

0225  1 
02208 

.02202 

02245 
.  02207 

.  02203 

02239 
. 02206 

.  02204 

.02235 
. 02205 

.  02205 

0223  1 
. 02204 
.  02207 

02228 
02203 
.  02209 

02225 

02202 

.0221 1 

02222 

02202 

.02215 

022  1  9 
.  0220  1 
022  1  9 

022  1  7 
0220  1 
02226 

022  15 
0220  1 
02234 

02248 

.02166 

02153 

.02186 

02224 

02164 

.02153 

.02198 

022  1  1 

02  162 

02  153 
022  1  9 

02203 
.02161 
02  154 

.02269 

02  197 
.02160 
02  154 

02191 

02  158 
.02155 

02187 

.02157 

.02156 

02  183 
.  02  1  S6 
.02157 

02  180 

.02156 

.02159 

02  177 

.02155 

02161 

02  174 

02154 

02164 

02  1  72 
.02154 
.02167 

02169 

.02153 

.02172 

02167 

02  153 
02  178 

02200 

.02120 
02  107 

.02132 

02  177 

.02118 

.02107 

.02140 

02165 

.02  1  17 
.02107 
02  152 

02  157 

.02115 
.02107 
02  173 

02  150 

.02  1  14 
.02108 

.02222 

02145 

.02113 

.02108 

02141 
02112 
02  109 

02  137 
.02111 
02  1  10 

02  134 
.021 10 

.02112 

02131 

.02109 
.02  1  13 

02  128 

02  109 

.02116 

02  126 

.  02  108 
02118 

02  124 

02  108 

.02122 

02  122 
02  107 

.02126 

02  155 
02076 

0206  3 
02082 

02  132 
.02075 

.  0206  3 
02088 

02  120 
.  02073 

.02063 
. 02096 

02112 
. 02072 

. 02063 
.02108 

02  106 
0207  1 
.  0206  4 

.02128 

02101 
.  0206  9 
. 02064 

.02176 

0209  7 
.  0206  8 

.  0206  5 

02093 
0206  7 
0206  5 

02090 

.  0206  7 
.  0206  7 

02087 
0206  6 
. 02068 

02085 
0206  5 
.02070 

02082 

.02064 

02072 

02080 
.  0206  4 

.  02075 

0207  8 

0206  4 
.02078 

02111 
02035 
0202  1 

0203  6 

02089 

.02033 
0202  1 
02040 

02077 
. 0203 1 
. 02021 

02046 

02070 

02030 
0202  1 
02054 

0206  4 

. 02029 
. 0202 1 
. 02066 

02059 
. 02028 
. 02022 
02086 

02055 

.  02027 
.  02022 
.02133 

0205  1 

02026 

02023 

0204  8 

. 02025 
. 02024 

0  2  0  4  5 
02024 

.02025 

02043 

02023 

02026 

02040 

02023 

.02028 

02038 
02022 
.  02030 

.02036 

02022 

02033 

02070 

0  19  9  4 

0  19  8  1 

0  19  9  2 

02048 

0  19  9  3 

0  19  8  1 

0  19  9  6 

02036 

0  19  9  1 

0  1981 
02000 

02029 

0  19  90 

0  19  8  1 
02006 

02023 

.01989 

.01981 

.02014 

020  1  8 

0  19  87 

.01981 

02025 

020  14 

.01986 

.01981 

.  02045 

020  10 

0  19  85 

0  1982 
0209  1 

02007 

.01985 

.01982 

02005 

.01984 

.01983 

02002 

0  19  83 

0  19  84 

02000 

0  19  8  2 
.01986 

0  19  9  8 

.01982 

0  19  88 

0  19  9  6 

.01981 
0  19  9  0 

’ 

Appendix  V:  Evaluation  of  moments  of  EV^.(£,b) 
By  definition,  for  p  >  0 


<Kp)  = 


r*(p) 

r(P) 


But,  for  any  positive  integer  k 


r(k)(p)  =  J  xp  1  exp(-x)(ln  x)kdx 
0 


and 


f(k)(i)  =  \  =  J  (£n  x)e  X  dx  . 

0 


i  -x 


In  general,  define  A  =  J  ( In  x)  e  dx.  In  particular  A.  =  T 

0 

=  -y  =  -.5772.  From  Jahnke  (1960)  (p.  12) 


=  rtt) 


-  l 


k=0  (z+k) 


and,  by  differentiating  successively  and  setting  z  =  1, 


k-1 

4>  '(z)  =  (-1)  k!C(k+l),  where  C(x)  is  Riemann's  zeta  function 


Using  values  of  C(k)  [Jahnke  (1960,  p.  37)],  we  have 


(1) 
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4> 


(1) 


(1) 


2 

n 

6 


<1* 


(2) 


(1) 


-2.404 


4> 


(3) 


(1) 


4 

n 

15 


Under  the  homogeneous  model,  we  obtain 

\i2  =  E  [  In  X  -  E(  An  X)]2  =  b2(X2~X2) 

H3  =  b3(X3-3X2X1+2X2)  =  -2.4036b3 

^4  =  b4(X4-4X3X1+6X2X2-3X^)  =  14.6119b4 


(see  [Menon  (1963)].  Here 


X2  =  y2  +  —  ,  X3  =  -5.4445,  X4  =  235601. 


However,  under  the  exchangeable  model,  using  Y  =  An  X, 


E  (Y)  =  ^{n  -  1  + 
bet  n 1  k* J 


2,  2 


Ehet(,2>  ■  -  1  +  7^  +  (“  -  1  + 

k* 


where 
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Ehet(Y3)  =  V-  <-»)V2)(l>  +  \  (fel)3-t-<2)(l) 


=  ♦(2>(1)('l,})  {„  -  i  +  _!_} 

n  k*3 


Ehet(Y4)  =  +  4(c*)V3)(l) 


n 


n  ^k* 


^(3W4){n.1+^} 

n  k* 


E  (Yr)  =  -br4>(r)(l)  [see  Patel  et  al  (1976)] 

homogeneous 


Appendix  VI  Computer  Programs 

Program  A  was  used  to  compute  u(r,n,k*)  for  the  Weibull 
distribution  using  the  recursive  formula  from  Chapter  V.3,  p.  103. 
Extended  precision  Fortran  was  used  to  perform  the  calculations. 


-188- 


189 


Program  A 

IMPLICIT  LOGICAL*l  (A-Z) 

DIMENSION  U  (50,50) ,  EFUN(501),  UBAR(50) 

REAL* 16  U,  EFUN,  UBAR ,  SUM,  K,  C,  Y,  YINCR 
INTEGER  I,  J,  N,  R 
C 

WRITE (6, 200) 

200  FORMAT ('  Enter  K  (with  decimal  point)  and  N 

1  1  (integer)'/  each  followed  by  a  comma.') 

READ (5, 100)  K,  N 
100  FORMAT (F7. 5, 12) 

WRITE (7,201)  K,  N 

201  FORMAT (' IK  =  ' , F7 . 5 , / '  N  =  ',12, /’O') 

C 

C  INITIALIZATIONS 
C 

DO  300  1=1, N 
DO  301  J= 1 , N 
U(J,I)  =  0.0Q0 
301  CONTINUE 
300  CONTINUE 
C 

YINCR  =  (1.Q0/500.Q0) 

DO  302  1=1,501 

Y  =  QFLOAT (l~l) *YINCR 
IF  (Y  .GT.  l.Q-30)  GO  TO  3021 
EFUN (I)  =  0.0Q0 
GO  TO  302 

1  EFUN (I)  =  QEXP (- (QABS (QLOG (1 . QO/Y) )  )  **  (1 . QO/K) ) 

CONTINUE 

CALCULATE  VALUES  OF  U-BAR  (FIRST  COLUMN  OF  U) . 

USE  SIMPSON'S  (1/3)  RULE  FOR  THE  INTEGRATION. 

DO  310  R=1 , N 

IF  (EFUN (1)  .LT.  l.Q-70)  EFUN(l)  =  l.Q-70 
SUM  =  EFUN ( 1 )  **  (R- 1 ) 

DO  311  J=1 , 249 

IF  (EFUN (2* J)  .LT.  l.Q-70)  EFUN (2* J)  =  l.Q-70 
IF  (EFUN (2*J+1)  .LT.  l.Q-70)  EFUN(2*J+1)  =  l.Q-70 
SUM  =  SUM  +  4.Q0*EFUN(2*J)**(R-1) 

1  +  2.Q0*EFUN(2*J+1)**(R-1) 

311  CONTINUE 

SUM  =  SUM  +  4.Q0*EFUN(500)**(R-1)  +  EFUN  (501)  **  (R-l) 
UBAR (R)  =  SUM* (. 002Q0/3 . QO) 

310  CONTINUE 
C 

C  SET  1ST  COLUMN  OF  U  EQUAL  TO  UBAR 
C 

DO  315  1=1, N 

U(I,1)  =  UBAR  (i) 

315  CONTINUE 


' 


non 
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WRITE  (7 , 250)  U  (1 , 1) 

250  FORMAT(14(’0'  ,  F8 . 5)  / ,  3  (3X,  14  (IX,  F8 . 5)  /)  /) 

CALCULATE  THE  REST  OF  U  AND  WRITE  IT  OUT. 

DO  320  1=2, N 
DO  321  R=2 , I 
SUM  =  0.Q0 
DO  322  J=1 , R 

CALL  COMB (R— 1 , J~1 , C) 

SUM  =  SUM  +  C*((-1)**(J"1))*UBAR(I-R+J) 
322  CONTINUE 

CALL  C0MB(I-1,R-1,C) 

U (I , R)  =  C*SUM 
321  CONTINUE 

WRITE  (7 , 250)  (U(I,R)  ,R=1, 1) 

320  CONTINUE 
STOP 
END 

SUBROUTINE  C0MB(l,J,C) 

IMPLICIT  REAL* 16  (A~H,K,0-Z) 

REAL* 16  K (15) 

DO  300  IK-1,15 
K (IK)  =  0.Q0 
300  CONTINUE 

LSTOP  =  MIN0(I-J,J) 

L=1 

PN=QFLOAT (L) 

PD=QFLOAT (L) 

IF  (LSTOP  .LE.  0)  GO  TO  3999 
DO  3000  L=l, LSTOP 
A=QFL0AT(L) 

PD=A*PD 
M=I~L+1 
A=QFL0AT (M) 

PN=A*PN 

DO  3001  IK-1,15 
K  (IK)  =  K  (IK)  + 1 
3001  CONTINUE 

IF  (K (1)  .NE.  2.Q0)  GO  TO  3101 
FACT=2.Q0 
PN=PN/FACT 
PD=PD/FACT 
K  (1)  =  0.Q0 
3101  CONTINUE 

IF  (K (2)  .NE.  3.Q0)  GO  TO  3102 
FACT=3.Q0 
PN=PN/FACT 
PD=PD/FACT 
K (2)  =  0.Q0 


' 


3102  CONTINUE 

IF  (K (3)  .NE.  5.Q0)  GO  TO  3103 
FACT=5 . QO 
PN=PN/FACT 
PD=PD/FACT 
K  (3)  =  O.QO 

3103  CONTINUE 

IF  (K (4)  .NE.  7 . QO)  GO  TO  3104 
FACT=7 . QO 
PN=PN/FACT 
PD=PD/FACT 
K  (4)  =  O.QO 

3104  CONTINUE 

IF  (K (5)  .NE.  11. QO)  GO  TO  3105 
FACT=1 1 . QO 
PN=PN/FACT 

pd=pd/fact 

K  (5)  =  O.QO 

3105  CONTINUE 

IF  (K (6)  .NE.  13. QO)  GO  TO  3106 
FACT=13 . QO 
PN=PN/FACT 
PD=PD/FACT 
K  (6)  =  O.QO 

3106  CONTINUE 

IF  (K(7)  .NE.  17. QO)  GO  TO  3107 
FACT=17 . QO 
PN=PN/FACT 
PD=PD/FACT 
K (7)  =  O.QO 

3107  CONTINUE 

IF  (K (8)  .NE.  19. QO)  GO  TO  3108 
FACT=19.Q0 

pn=pn/fact 
pd=pd/fact 
K(8)  =  O.QO 

3108  CONTINUE 

IF  (K (9)  .NE.  23. QO)  GO  TO  3109 
FACT=23 . QO 
PN=PN/FACT 

pd=pd/fact 
K(9)  =  O.QO 

3109  CONTINUE 

IF  (K (10)  .NE.  29. QO)  GO  TO  3110 
FACT=29.Q0 

pn=pn/fact 

PD=PD/FACT 
K  (10)  =  O.QO 


3110  CONTINUE 

IF  (K (1 1)  .NE. 
FACT=31 . QO 
PN=PN/FACT 
PD=PD/FACT 
K (1 1)  =  O.QO 

3111  CONTINUE 

IF  (K(12)  .NE. 
FACT=37.Q0 

pn=pn/fact 
pd=pd/fact 
K  (12)  =  O.QO 

3112  CONTINUE 

IF  (K (13)  .NE. 
FACT=41.Q0 
PN=PN/FACT 
PD=PD/FACT 
K  (13)  =  O.QO 

3113  CONTINUE 

IF  (K (14)  .NE. 
FACT=43 . QO 
PN=PN/FACT 
PD=PD/FACT 
K  (14)  =  O.QO 

3114  CONTINUE 

IF  (K (15)  .NE. 
FACT=47.Q0 

pn=pn/fact 

PD=PD/FACT 
K  (15)  =  O.QO 

3115  CONTINUE 
3000  CONTINUE 
3999  CONTINUE 

C=PN/PD 

RETURN 

END 


31. QO)  GO  TO  3111 


37. QO)  GO  TO  3112 


41. QO)  GO  TO  3113 


43. QO)  GO  TO  3114 


47. QO)  GO  TO  3115 


. 


Program  B  was  used  to  generate  25  random  samples  each  of  size  5 
with  one  outlier,  starting  from  a  uniform  (0,1)  distribution.  The 
samples  of  Weibulls  are  printed  out  as  X’s  and  transformed  to  W’s, 


samples  of  EV^.  ,  by  a  In  transformation.  For  each  sample  w  and 
2 

s^  are  computed.  Also  for  each  sample  the  smallest  observation  is 

replaced  by  the  second  smallest  and  by  the  largest  and  the  new  means 
-  2  2 

w,~  ,  N ,  w,  and  variances  s,»  ,  N  ,  s.  are  computed.  Also  the 

(2,1)’  (n, 1)  (2,1)  (n, 1 ) 

largest  observation  is  replaced  by  the  second  largest  observation  and 

by  the  smallest  observation  and  again  means  w,  -  N  ,  w,,  N  and 

(n-1 ,n)  (1 ,n) 

2  2 

variances  s,  ,  N  ,  s,,  N  are  computed  for  each  sample. 

(n-1 ,n)  (l,n) 
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Program  B 

C  This  program  computes  Weibull  and  EV  variables, 

C  every  fifth  one  spurious. 

C  It  prints  out  samples  of  size  5 

C  along  with  mean  and  variance. 

C  For  each  sample  it  also  computes 

C  the  Winsorized  mean  and  variance 

C  by  replacing  the  smallest 

C  observation  or  the  largest  observation. 

IMPLICIT  REAL *8  (A-H,0~Z) 

DIMENSION  X(125),  W(125),  W21 (125) ,  W51(125), 
W45  (125)  ,  W15C125) 

DIMENSION  WB(25),  WV (25) ,  W21B(25),  W21V(25), 
W51B  (25) ,  W51V  (25) 

DIMENSION  W45B (25) ,  W45V(25),  W15B(25),  W15V(25) 
REAL  RAN (125) 

INTEGER  ISEED (2) 

INTEGER* 2  HSEED (3) 

REAL *8  DSEED 

EQUIVALENCE  (ISEED (2) , HSEED (1) ) 

C 

WRITE (6, 200) 

200  FORMAT (’  Enter  BETA  and  K') 

READ (5 , 100)  BETA,  C 

100  FORMAT (2F12. 5) 

WRITE (7 , 201)  BETA,  C 

201  FORMAT  ('1  BETA  =  \F12.5,'  K=  '  ,F7.5) 

BETAC  =  BETA*C 

NSIZE  =  5 
NS AMP  =  25 

NTOTAL  =  NSIZE*NSAMP 


Get  random  number  depending  on  time  of  day  into  DSEED 


CALL  TIME (4,0, ISEED) 

HSEED (3)  =  HSEED (1) 

HSEED (1)  =  HSEED (2) 

HSEED (2)  =  HSEED (3) 

DSEED  =  DABS (DFL0AT (ISEED (2))) 

C 

C  Call  *IMSL  routine  GGUBS  to  get  NTOTAL  random 
C  numbers  into  RAN 
C 

NR  =  NTOTAL 

CALL  GGUBS (DSEED,  NR,  RAN) 


' 


noon  non  nn 
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Compute  Weibull  r.v.'s  X  and  W 

DO  301  K=1 , NTOTAL 
FACT  =  BETA 

IF  (MOD (K, 5)  .EQ.  0)  FACT  =  BETAC 
X  (K)  =  DBLE  (-AL0G  (RAN  (K)  )  )  **  (1 .  DO/FACT) 
W(K)  =  DLOG  (X  (K)  ) 

301  CONTINUE 


Print  out  X  and  W 
WRITE  (7, 202) 

202  FORMAT('-  X'/) 

WRITE (7, 203)  X 

203  FORMAT (25 (3X, 5 (IX, E13 . 6) /) ) 

WRITE  (7, 204) 

204  FORMAT (’0  W') 

WRITE  (7, 203)  W 

Compute  Means  and  Variances  for  each  sample  of  size 
NSIZE  in  X  and  W 

K  =  1 

DO  310  KK=1 ,NSAMP 

CALL  STAT (W ,  K,  NSIZE,  WB (KK) ,  WV (KK) ) 

K  =  K+NSIZE 
310  CONTINUE 
C 

C  Find  the  smallest,  second  smallest,  largest,  and  second 
C  largest  entry  in  each  of  the  samples  in  W. 

C 

K  =  1 

DO  320  KK= 1 , NS AMP 
WMINl  =  1.D75 
WMAXl  =  -1.D75 
DO  321  Kl=l, NSIZE 

IF  (W(K)  .GE.  WMINl)  GO  TO  3211 
WMIN2  =  WMINl 
WMINl  =  W(K) 

KMIN  =  K 
GO  TO  3212 

3211  CONTINUE 

IF  (W(K)  .GE.  WMIN2)  GO  TO  3212 
WMIN2  =  W(K) 

3212  CONTINUE 
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C 

IF  (W(K)  .LE.  WMAXl)  GO  TO  3216 
WMAX2  =  WMAX1 
WMAXl  =  W(K) 

KMAX  =  K 
GO  TO  3217 

3216  CONTINUE 

IF  (W(K)  .LE.  WMAX2)  GO  TO  3217 
WMAX2  =  W(K) 

3217  CONTINUE 
K  =  K  +  1 

321  CONTINUE 
C 

C  Replace  the  smallest  and  largest  entries  with 
C  the  second  smallest 

C  and  largest  and  with  the  second  largest  and 
C  smallest  entries  to 

C  form  the  altered  samples  W(21),  W(51),  W(45),  and  W(15). 
C  Recompute  the  mean  and  variance  for  each  of 
C  the  altered  samples. 

C 


K  =  K  -  NSIZE 
DO  322  Kl=l, NSIZE 
W21  (K)  =  W(K) 

IF  (K  .EQ.  KMIN)  W21 (K) 
W51  (K)  =  W(K) 

IF  (K  .EQ.  KMIN)  W51 (K) 
W45(K)  =  W(K) 

IF  (K  .EQ.  KMAX)  W45 (K) 
W15  (K)  =  W(K) 

IF  (K  .EQ.  KMAX)  W15 (K) 
K  =  K  +  1 
322  CONTINUE 
320  CONTINUE 


WMIN2 

WMAXl 

WMAX2 

WMINl 


C 


330 


K  =  1 

DO  330  KK=1 , NSAMP 

CALL  STAT (W21 ,  K,  NSIZE, 
CALL  STAT (W5 1 ,  K,  NSIZE, 
CALL  STAT (W45 ,  K,  NSIZE, 
CALL  STAT (W15 ,  K,  NSIZE, 
K  =  K  +  NSIZE 
CONTINUE 
WRITE  (7, 225) 


W21B (KK) ,  W21 V (KK) ) 
W51B (KK) ,  W51V (KK) ) 
W45B(KK),  W45V (KK) ) 
W15B(KK),  W15V (KK) ) 


225  FORMAT (' $**$FORMAT=FMTLl  FONTNEXTIMAGE= 

1200. MEDIUM. 9. FIXED. LANDSCAPE. 1 ') 

WRITE (7, 250) 

250  FORMAT ('1  Sample  WBAR 

WVAR  W  (21)  BAR  W  (21)  VAR  W(51) 

BAR  W (51) VAR  W (45) BAR  W (45) VAR 
W  (15)  BAR  W(15)VAR 
DO  350  J=1 , NSAMP 
C 

C  Reassign  these  so  the  WRITE  statement  will 
C  fit  on  one  line. 

C 

W1  =  WB ( J) 

W2  =  WV(J) 

W3  =  W21B ( J) 

W4  =  W21 V ( J) 

W5  =  W51B (J) 

W6  =  W51V  (J) 

W7  =  W45B  ( J) 

W8  =  W45V ( J) 

W9  =  W15B  ( J) 

W10  =  W15V (J) 

WRITE (7 , 255)  J,  Wl ,  W2,  W3,  W4,  W5,  W6, 

W7,  W8,  W9,  W10 

255  FORMAT (6X , I 2 , 2X , 5 ( '  ' , 2 (F8 . 5 , IX) ) , ' ' ) 

350  CONTINUE 
STOP 
END 

SUBROUTINE  STAT(X,  KK,  NSIZE,  XBAR ,  XVAR) 
IMPLICIT  REAL *8  (A-H.O-Z) 

DIMENSION  X (125) 

KSTOP  =  KK  +  NSIZE  -  1 
SUM  =  0.0D0 
SSQ  =  0.0D0 
DO  300  K=KK, KSTOP 
SUM  =  SUM  +  X  (K) 

SSQ  =  SSQ  +  (X(K)**2) 

300  CONTINUE 

XBAR  =  SUM/DFLOAT (NSIZE) 

XVAR  =  (SSQ  -  ( (SUM**2) /DFLOAT (NSIZE)))/ 
(DFLOAT (NSIZE)  -  1 

RETURN 

END 


